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Abstract 

The  saturation  of  raicroinstabi 1 i t ies  has  been  theorized  to 
cause  anomalous  transport  in  MPD  thrusters.  However,  calculation 
of  anomalous  transport  properties  1.'  an  area  which  is  still  in 
its  infancy.  It  is  the  object  of  this  paper  to  discuss  the 
importance  of  anomalous  transport  in  MPD  thruster  operation,  how 
microinstabilities  are  related  to  anomalous  transport,  the  types 
of  microinstabilities  most  likely  to  occur,  and  at  what  locations 
in  the  MPD  thruster.  First,  preliminary  measurements  were  made 
of  mean  plasma  properties  and  fluctuations  using  imbedded  probe 
techniques  in  the  plume  and  internal  locations  of  a  MW  level  MPD 
thruster  operating  near,  but  below,  global  onset  conditions  using 
probe  techniques.  Analysis  of  these  results  as  well  as 
significant  incidental  results  are  discussed.  Next,  a  kinetic 
theory  method  for  calculating  plasma  transport  properties  in  a 
two- temperature,  multiple  species,  non-equilibrium  plasma  has 
been  modified  to  account  for  anomalous  transport.  The 
differences  between  transport  coefficients  calculated  using  that 
method  and  a  method  based  on  the  mean  free  path  formula  for 
conductivity  have  been  compared.  In  addition,  as  an  illustrative 
example  of  application  of  the  kinetic  method,  anomalous  and 
classical  properties  have  been  calculated  at  different  spatial 
locations  in  the  MW  level  MPD  thruster  plasma  using  the  mean 
plasma  properties  acquired  through  measurements  in  the 
preliminary  experiments. 

The  experimental  results  did  not  show  oscillations  to 
support  the  presence  of  microinstabilities  with  characteristic 
frequencies  near  the  lower  hybrid  frequency,  which, 
theoretically,  should  have  been  present.  However,  much  lower 
frequency  oscillations  were  observed,  which  have  been  theorized 
to  be  associated  with  the  onset  phenomenon.  In  addition,  the 
insertion  of  the  Langmuir  probe  into  the  thruster  current  path 
appeared  to  induce  onset  at  operating  conditions  where  onset 
should  not  normally  occur.  Based  on  the  measured  mean  plasma 
properties  in  the  MW  thruster,  the  computational  work  shows  that 
anomalous  conductivity  does  not  appear  to  be  as  significant  as 
previously  suspected.  This  difference  appears  to  stem  from  the 
inclusion  of  more  particle  interactions  in  the  kinetic  model, 
specifically  electron-electron  interactions. 
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I.  INTRODUCTION 

In  the  past,  the  movement  of  objects  in  the  vacuum  of  space 
has  relied  on  the  use  of  chemical  rocket  engines.  Unfortunately, 
the  majority  of  the  mass  in  a  chemically  propelled  spacecraft 
consists  of  the  chemical  propellants  themselves.  In  order  to 
reduce  the  mass  of  the  spacecraft  without  reducing  thrust,  the 
exhaust  velocity  of  the  exiting  gases  must  be  increased. 

Currently,  chemical  rocket  engines  are  limited  to  exhaust 
velocities  between  4  and  5  km/sec  due  to  temperature  and  pressure 
limitations  on  the  combustion  chambers,  as  well  as  the 
limitations  on  energy  which  propellants  can  add  to  the  flow 
through  exothermic  chemical  reactions.  Thrusters  using  an 
external  power  supply  to  increase  the  thermal  energy  of  the 
propellant  are  generally  limited  to  exit  velocities  of  10  km/sec 
before  temperatures  threaten  the  structure  of  the  engine  [1]. 

However,  if  electromagnetic  body  forces  are  used  to 
accelerate  the  propellant  instead  of  the  conventional  means  of 
heating  and  expanding  through  a  nozzle,  exhaust  velocities  an 
order  of  magnitude  greater  than  chemical  or  electrothermal 
devices  are  possible.  One  type  of  device  which  uses 
electromagnetic  body  forces  to  produce  thrust  is  the  self-field 
magnetopl  ismadynamic  (MPD)  thruster. 

Figure  1.1  shows  a  schematic  depiction  of  a  self-field  MPD 
thruster.  The  self-field  thruster  is  very  simple  in  its  basic 
design.  Propellant  is  injected  at  the  back  of  the  thruster,  and 
a  high  current  electric  arc  is  produced  between  a  cathode  and  a 


Figure  1.1  -  Schematic  depiction  of  self-field  MPD  thruster 

operation. 

coaxial  anode.  When  the  anode  and  cathode  are  coaxial,  the 
current  produces  a  magnetic  field,  B,  primarily  in  the  azimuthal 
direction  while  ionizing  the  propellant.  As  the  current,  I, 
travels  radially,  the  I  x  B  interaction  provides  an 
electromagnetic  body  force  which  accelerates  the  ionized  plasma 
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in  an  axial  direction  while  simultaneously  pushing  the  pLasma 
radially  inward.  This  type  of  acceleration  may  also  be  referred 
to  as  Lorentz  force  acceleration  [2].  The  resultant  exhaust 
velocities  are  in  the  range  of  10-50  km/sec,  depending  on  the 
current  level  and  propellant  used  [3].  For  a  steady-state 
device,  a  separate  power  supply  provides  the  current,  while  in 
the  case  of  a  quasi-steady  thruster,  a  large  capacitor  bank 
provides  the  current.  A  quasi-steady,  MW  level  self-field  MPD 
thruster  is  the  device  examined  in  this  research. 

Currently,  an  MPD  thruster  has  very  low  efficiencies,  on 
the  order  of  30  per  cent  at  most,  where  a  minimum  60  per  cent 
efficiency  is  required  to  make  the  MPD  thruster  competitive  with 
propulsion  systems  currently  in  use  [3].  What  this  means  is  that 
only  30  per  cent  of  the  input  energy  is  becoming  converted  into 
directed  kinetic  energy,  which  is  the  energy  used  to  produce 
thrust.  The  question  is,  what  happens  to  the  other  70  per  cent 
of  the  energy?  Recently,  Choueiri,  et.  al .  has  suggested  that  a 
portion  of  the  energy  which  is  not  converted  into  thrust  provides 
the  energy  source  for  microscopic  plasma  instabilities,  or 
microinstabilities  [4]. 

The  primary  concern  about  microinstabilities  is  that  they 
have  been  theorized  to  dictate  anomalous  transport.  Anomalous 
transport  may  be  defined  as  additional  effects  not  accounted  for 
in  the  classical  kinetic  theory  of  gases  that  increase  such 
properties  as  viscosity,  while  decreasing  others  such  as 
electrical  conductivity.  This  anomalous  transport  leads  to 
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electrical  conductivity.  This  anomalous  transport  leads  to 
undesirable  plasma  heating  due  to  additional  viscous  dissipation, 
more  ohmic  heating,  etc.  As  the  MPD  thruster  uses  primarily 
electromagnetic,  rather  than  electrothermal  energy,  to  provide 
thrust,  the  energy  going  into  heating  the  plasma  is  essentially 
wasted.  This  wasted  energy  results  in  the  very  low  efficiencies 
associated  with  the  MPD  thruster.  It  is  hoped  that  a  better 
understanding  of  microinstabilities  will  enable  the  scientific 
community  to  determine  the  significance  of  the  impact  which 
microinstabilities  have  on  these  energy  losses. 

1.1  Definition  of  Microinstability 

An  instability  may  be  described  as  an  oscillation  or  other 
disturbance  in  one  or  more  parameters  of  a  system  that  changes  in 
time  and  space  [2].  The  progression  of  an  instability  is  easy  to 
understand  in  a  continuum  flow.  For  example,  turbulent  flow  may 
be  considered  as  resulting  from  an  instability;  eddies  can  be 
said  to  form  and  grow  with  time  and  distance  downstream.  Such 
descriptions  of  fluid  turbulence  are  based  on  the  existence  of  a 
continuum  of  particles.  In  contrast,  a  plasma  can  be  considered 
as  either  a  continuum  when  the  plasma  state  is  highly 
collisional,  or  it  may  be  considered  as  a  group  of  independent 
particles,  for  example  when  the  plasma  state  indicates  that  it  is 
collisionless.  In  the  latter  case,  internal  electric  fields  can 
play  a  role  in  creating  and  propagating  disturbances. 

A  microinstability  can  be  described  simply  as  a  plasma 
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instability  that  occurs  on  a  very  small  scale.  For  continuum 
state  instabilities,  the  smallest  eddies  are  still  much  larger 
than  a  mean  free  path  between  molecular  or  atomic  collisions; 
the  eddies  and  their  interactions  can  be  considered  to  occur  on  a 
macroscopic  scale.  In  a  plasma,  electric  fields  can  play  a  role 
in  instabilities.  A  plasma  typically  observes  the  quasi¬ 
neutrality  condition,  that  is  there  are  an  equal  number  of 
negative  and  positive  charges  in  a  plasma  and  macroscopically  the 
plasma  appears  to  be  neutral.  However,  because  the  charge  in  a 
plasma  is  associated  with  different  species,  there  may  be 
specific  locations  where  effects  due  to  a  net  electric  field 
appear.  These  effects  do  not  extend  very  far;  they  most  likely 
occur  only  in  small  zones  where  there  is  a  local  charge 
separation.  As  a  result,  these  effects  occur  on  a  length  scale 
that  is  often  less  than  the  mean  free  path,  A.  [4],  which  is  on 
the  order  of  a  millimeter  for  a  plasma  with  a  number  density,  ng, 
of  about  1015  m'^  and  electron  temperature,  Te,  about  3  eV.  This 
diminutive  length  scale  gives  microinstabilities  their  name. 

One  general  class  of  microinstabilities  is  the  category  of 
streaming  instabilities.  A  streaming  instability  may  exist 
whenever  a  relative  drift  between  different  species  occurs  [6]. 

As  the  current  which  provides  the  electromagnetic  acceleration 
also  produces  a  relative  drift  between  electrons  and  ions, 
streaming  instabilities  are  expected  to  occur  in  the  MPD 
thruster.  A  further  discussion  of  streaming  instabilities  is 
presented  in  Chapter  II. 
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1.2  Previous  Work 

Microinstabilities  in  I  x  B  plasmas  have  been  studied 
extensively  by  the  plasma  physics  community  [7].  Unfortunately, 
many  of  these  works  are  not  directly  applicable  to  MPD  devices 
because  the  plasma  properties  in  other  plasma  devices  differ 
significantly  from  the  MPD  thruster  [1]. 

Initial  work  concerning  the  effects  of  instabilities  in  the 
MPD  plasma  concentrated  primarily  of  the  prediction  of  the 
phenomenon  known  as  "onset".  "Onset"  is  the  operating  point 
where  various  processes  detrimental  to  the  MPD  thruster,  such  as 
cathode  erosion,  begin  to  occur.  This  point  is  defined  as  the 
global  current  and  global  voltage  level  above  which  high 
frequency  terminal  voltage  fluctuations  are  observed,  and  this 
point  scales  with  I2/ ft  ,  where  ft  =  mass  flow  of  propellant [8, 9] . 
Recently,  macroscopic  instabilities  have  been  examined  as  the 
cause  of  onset  [10-14],  where  a  macroscopic  electron  acoustic 
instability  [10]  or  a  Pierce  instability  [13,14]  appears  to  be 
responsible  for  onset.  Experimentally,  inconclusive  attempts 
have  been  made  to  identify  the  instability  associated  with  onset 
phenomenon  by  examining  the  power  spectra  of  both  terminal 
voltage  and  light  intensity  f luctuations[ 15] . 

In  recent  years,  however,  some  emphasis  has  begun  to  be 
placed  on  the  investigation  of  microinstabilities  in  MPD 
thrusters  operating  below  the  onset  current  level. 
Microinstabilities  have  been  studied  extensively  for  other  plasma 


7 


devices[ 16-25  ] ,  but  recent  work  has  focused  on  the  theoretical 
identification  of  the  various  instabilities  which  may  occur  in 
the  self-field  MPD  thruster  [4,26,27]  and  the  identification  of 
the  instability  which  dictates  the  anomalous  transport  properties 
(e.g.  electrical  conductivity,  viscosity,  thermal  conductivity, 
etc.)  in  the  plasma  of  the  MPD  thruster.  At  the  time  of  this 
study,  the  generalized  lower  hybrid  drift  instability  ( GLHDI )  has 
been  linked  to  anomalous  transport  [4,28,29].  A  more  detailed 
discussion  of  the  GLHDI,  as  well  as  other  microinstabilities,  is 
presented  in  Chapter  II. 

Experimental  verification  of  both  the  presence  of 
microinstabilities  in  the  MPD  thruster  as  well  as  a  measurement 
of  the  linear  dispersion  relation  of  the  GLHDI  has  also  been 
attempted,  most  notably  at  Princeton  University  [1,4,30,31]. 

1.3  Thruster  Efficiency 

Before  discussing  the  effects  of  microinstabilities  on 
efficiency,  it  is  first  necessary  to  look  at  what  defines 
efficiency.  There  are  two  modes  of  input  energy  conversion  which 
operate  in  an  electromagnetic  thruster:  plasma  heating  and 
electromagnetic  power.  Most  of  the  electromagnetic  power  is 
converted  into  thrust  with  the  exception  of  some  undirected  flow 
[2].  However,  the  plasma  heating  mode  has  very  little  of  its 
portion  of  total  energy  go  into  the  form  of  directed  kinetic 
energy.  The  majority  of  the  power  input  for  heating  is 
dissipated  in  electrode,  plasma  thermal,  and  internal  mode 
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(vibration,  ionization)  losses,  and  also,  in  undirected  kinetic 
energy  [2].  Because  the  MPD  thruster  has  a  low  operating 
pressure,  the  stabilizing  effect  of  collisions  is  not  readily 
available.  Thus,  the  flow  in  a  MPD  thruster  is  in  a  highly  non¬ 
equilibrium  state,  and  this  fact  makes  it  much  more  difficult  to 
convert  heat  into  propulsive  kinetic  power  [2].  As  a  result  the 
MPD  thruster  efficiency  becomes  reduced.  Meanwhile,  the  thermal 
energy  which  does  not  get  converted  into  propulsive  kinetic  power 
serves  as  the  energy  source  for  microinstabilities.  If  the 
microinstabilities  can  be  re<  ced,  in  principle,  more  energy  will 
be  available  for  conversion  into  directed  kinetic  energy. 

Two  methods  can  be  used  to  increase  the  efficiency,  assuming 
that  electromagnetic  power  losses  are  negligible  in  comparison. 
One  method  is  to  try  to  increase  the  percentage  of  electrothermal 
thrust  by  incorporating  such  components  as  a  thrust  nozzle.  This 
is  attempted  in  hybrid  arcjet-MPD  thrusters.  In  the  second 
method,  the  overall  plasma  heating  is  reduced.  One  means  of 
accomplishing  this  is  to  control  the  occurrence  of 
microinstabilities.  By  controlling  the  microinstabilities,  the 
additional  plasma  heating  associated  with  anomalous  transport  is 
reduced,  thus  decreasing  the  wasted  input  power. 
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Assuming' an  azimuthally  symmetric,  radial  discharge,  the 
thruster  efficiency  can  be  written  as  follows  [2]: 


n 


( i .  n 


In  the  preceding  equations,  r|  =  efficiency,  o  =  electrical 
conductivity,  j  =  current  density,  V  =  volume  (a  volume 
integral),  th  =  mass  flow  and  the  electromagnetic  thrust 
coefficient,  b,  is 


(1.2) 


where  r,  =  radius  of  anode,  r.  =  radius  of  cathode,  and  = 
permeability  of  free  space.  Considering  the  expression  for 
efficiency,  one  can  see  that  as  o  decreases,  the  second  term  in 
the  denominator,  representing  the  increase  in  thermal  energy  due 
to  ohmic  heating,  will  increase,  decreasing  the  efficiency.  The 
classical  value  of  electrical  conductivity  can  be  expressed  by 
the  mean  free  path  formula: 


where  or-:  h  classical  electrical  conductivity,  e  a  charge  of  an 
electron,  3  mass  of  an  electron,  and  voi  a  electron-ion 
collision  frequency. 

The  microinstabilities  can  be  expected  to  alter  the 
effective  collision  frequency,  thereby  causing  changes  in  the 
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electrical  conductivity.  However,  these  instabilities  will  not 
grow  indefinitely.  Eventually,  mechanisms  appear  which  partially 
or  totally  drain  the  energy  source  fueling  the  instability.  When 
this  occurs,  the  plasma  will  reach  a  new  state  of  equilibrium, 
and  often  this  steady-state  occurs  so  as  to  allow  both  the  energy 
source  and  the  sink  (the  instability)  to  exist  simultaneously 
[2],  This  new  steady-state  is  known  as  saturation  and  is  of  a 
turbulent  nature  [2].  As  linearized  dispersion  relations  predict 
unlimited  growth  of  the  instability  as  time  increases,  to  produce 
saturation  a  non-linear  damping  mechanism  must  be  present  which 
is  not  accounted  for  by  linearized  theory.  One  type  of  non¬ 
linear  damping  mechanism  is  a  collision,  and  its  effect  is 
expressed  in  terms  of  the  particle  collision  frequency.  However, 
collisions  are  not  always  enough  to  stop  the  runaway  growtn  of  an 
instability.  In  that  case,  an  interaction  between  the  electric 
field  produced  by  the  unstable  wave  and  a  charged  particle  must 
provide  the  mechanism  for  saturation  of  the  instability. 

1.4  Saturating  Instabilities 

There  are  many  possible  mechanisms  by  which  such  saturation 
can  occur.  Three  of  the  most  common  are  quasi linear 
stabilization,  saturation  due  to  resonance  broadening,  and 
saturation  due  to  electrostatic  trapping.  The  firnt  two  are 
supposed  to  have  a  negligible  effect  on  microinstabilities  in  an 
MPD  thruster  [4].  Electrostatic  trapping  can  be  described  as  the 
situation  where  a  charged  particle  is  travelling  at  a  velocity 


near  the  wave  velocity  and  is  trapped  in  a  potential  "well"  of 
the  wave.  Any  charged  particle  which  is  trapped  is  forced  to 
slow  down  when  it  is  travelling  upwards  in  a  potential  "well", 
and  increases  in  speed  as  it  travels  down  the  well.  Some  energy 
transfer  occurs  between  the  wave  and  particle,  tending  to 
equalize  the  two  velocities  [4].  One  can  think  of  a  ball  in  a 
bowl.  If  the  ball  is  let  go  near  the  top,  it  will  roll  almost  as 
high  up  on  the  other  side,  and  then  roll  back,  and  so  on. 
Eventually  the  initial  potential  energy  is  dissipated  through 
friction  and  the  ball  settles  to  the  bottom  of  the  bowl.  In  a 
similar  manner,  some  of  the  particle  excess  energy  is  transferred 
to  the  wave,  and  vice  versa  in  an  oscillating  manner,  until 
eventually  the  particle  velocity  matches  the  wave  phase  velocity. 
In  addition,  the  wave  of  one  species  acts  on  the  other  species  in 
the  plasma,  effectively  transferring  energy  from  one  particle  to 
another.  This  may  explain  plasma  heating:  energetic  electrons 
transfer  their  energy  to  ions  through  wave  interactions  as  the 
instability  saturates  in  an  attempt  to  produce  a  plasma  in 
thermal  equilibrium. 

In  general,  the  effective  collision  frequency  depends  on  the 
method  of  non-linear  damping  assumed  to  offset  the  linear  growth. 
Choueiri  has  professed  that  if.  ion  trapping  is  the  method  of 
damping,  collision  frequencies  can  be  as  much  as  two  orders  of 
magnitude  higher  than  the  classical  collision  frequency  [4]. 

This  is  the  mode  of  non-linear  damping  which  he  expects  to 
produce  anomalous  transport  in  the  MPD  thruster  [4].  On  the 
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other  hand,  Choueiri  has  also  shown  that  if  resonance  scattering 
or  electron  trapping  is  the  damping  method,  the  anomalous 
collision  frequency  can  become  much  less  than  the  classical 
value.  Thus,  the  change  in  transport  with  these  types  of  damping 
mechanisms  becomes  fairly  insignificant  [4], 

Whichever  type  of  non-linear  damping  mechanism  is  used  to 
account  for  anomalous  transport,  the  primary  conditions  under 
which  the  microinstability  is  expected  to  occur  must  be  known. 

The  ratio  of  electron  drift  velocity  to  ion  thermal  velocity  is  a 
key  parameter  in  the  conditions  for  the  occurrence  of 
microinstabilities  [4,15].  As  current  increases  above  a  critical 
value  of  electron  drift  to  ion  thermal  velocity  (about 
1.5)[4/15],  the  collision  frequency  stays  fairly  constant.  The 
mechanism  of  ion  trapping  is  characterized  more  by  the  ratio  of 
electron  collision  frequency  over  the  lower  hybrid  frequency,  or 
as  the  inverse  of  the  Hall  parameter  scaled  by  the  square  root  of 
the  mass  ratio.  Thus,  the  high  Hall  parameter  regions  of  the 
discharge  are  the  regions  where  anomalous  dissipation  is  more 
likely  to  occur  [4].  One  such  region  would  be  near  the  anode, 
where  the  number  density,  and  consequently,  collision  frequency 
is  much  lower.  The  electron  Hall  parameter  is  defined  as  [19]: 

v  eH 

where  a>ces  electron  cyclotron  frequency  (hz),  and  v 9H  s  electron- 
heavy  species  collision  frequency.  Thus  lower  collision 
frequencies  lead  to  higher  hall  parameters,  which  would  result  in 
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stronger  anomalous  transport  effects  due  to  ion  trapping.  It  has 
also  been  found  that  collisional  processes  tend  to  disrupt  these 
non-linear  collisionless  damping  mechanisms,  allowing  ion¬ 
trapping  also  to  produce  effective  collision  frequencies  below 
the  classical  value  even  in  collision  cominated  regions  [ 4 ] . 

The  saturation  of  microinstabilities  seems  to  be  the  major 
player  in  producing  anomalous  transport,  and  thus  decreasing  the 
plasma  electrical  conductivity  below  the  classically  predicted 
values.  This  reduced  electrical  conductivity  produces  a  larger 
voltage  drop  for  a  given  current,  which  means  more  power  must  be 
added  to  the  system  to  keep  it  operating  at  its  current  level. 

If  the  plasma  is  fully  ionized,  the  ratio  of  electron  drift 
velocity , U.-e ,  to  ion  thermal  velocity,  v.  , ,  usually  exceeds  the 
threshold  value  of  1.5  for  the  occurrence  of  microinstabilities. 
Thus,  another  condition  which  for  all  practical  purposes  may  be 
considered  as  a  necessary,  but  not  sufficient  condition  for  the 
presence  of  microinstabilities,  is  that  the  plasma  is  fully 
ionized.  Introducing  the  dimensionless  total  current[20] 

5  s  (1.5) 

^  f  2 

where  Ifi  *  current  required  for  full  ionization  of  the  plasma, 
microinstabilities  may  occur  if  £>  1  [20].  The  utility  of  using 
this  parameter  is  that  the  occurrence  of  microinstabilities  can 
be  predicted  to  some  degree  based  on  global  parameters,  expressed 
in  an  alternate  form  of  (I. 5) [4, 20]: 

In  the  above  equation,  =  ionization  potential,  M  =  atomic 
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mass,  and  p  is  the  magnetic  permeability  constant,  which  for  all 
practical  purposes  is  the  permeability  of  free  space  [4].  An 
additional  use  of  ;  is  to  determine  how  much  voltage  is  required 
to  produce  electromagnetic  acceleration.  This  relationship  is 
expressed  as: 

for  £  <  1 ,  <$>«£ 

for  £>1,  (1.7) 

4>  = 

where  V'3*  =  total  input  voltage  and  V:;glf  s  the  portion  of  total 
voltage  which  produces  electromagnetic  acceleration  at  a  current 
of  I.-;.  As  thrust  scales  linearly  with  current  (F  a  I  x  B),  the 
more  voltage  which  is  required  to  produce  that  thrust  causes  the 
input  power  required  to  increase  as  well.  Thus,  an  increase  in 
input  power  without  a  corresponding  increase  in  thrust  will 
result  in  a  lower  thruster  efficiency.  Consequently,  efficiency 
decreases  as  £  increases  above  1.  Since  microinstabilities  are 
likely  to  occur  when  5  is  above  1,  and  efficiency  drops  at  that 
condition  as  well,  one  may  infer  that  microinstabilities  may  be  a 
cause  for  the  loss  of  engine  efficiency  once  the  plasma  becomes 
fully  ionized. 


1.5  Outline  of  Report 


The  current  investigation  is  designed  to  accomplish  two 
major  tasks:  1)  To  measure  mean  plasma  properties  and  plasma 
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fluctuations  to  determine  if  and  where  microinstabilities  appear 
in  the  thruster,  and  2)  to  develop  an  improved  method  for 
calculation  of  anomalous  transport  properties  based  on  mean 
plasma  properties  measured  inside  the  thruster. 

To  accomplish  these  two  tasks,  experiments  were  performed  on 
a  MW  level  MPD  device  to  gain  spatial  maps  of  the  magnetic  field, 
the  electron  number  density,  the  electron  temperature,  and  the 
relative  magnitudes  of  electron  number  density  oscillations. 

These  measurements  were  obtained  using  magnetic  and  Langmuir 
prooes  to  examine  the  possible  occurrence  of  microinstabilities. 
The  data  also  serves  as  inputs  to  a  kinetic  theory  computer  code 
for  determining  plasma  transport  properties  taking  into  account 
microinstabilities,  for  example  those  that  occur  at  the  lower 
hybrid  frequency. 

Chapter  II  will  discuss  the  types  of  microinstabilities 
which  are  expected  to  occur  in  an  MPD  thruster  as  well  as  discuss 
the  relationships  between  them.  Chapter  III  presents  the  results 
of  experiments  performed  in  a  quasi-steady  MW  level  self-field 
MPD  device.  Chapter  IV  examines  the  modification  of  a  plasma 
transport  code  to  account  for  anomalous  transport,  and  conditions 
from  two  spatial  locations  determined  through  the  experimental 
work  discussed  in  chapter  III  are  examined.  Finally,  chapter  V 
discusses  ways  to  improve  the  experiments  and  computer  program, 
points  out  avenues  of  possible  future  research  and  draws  some 
tentative  conclusions  from  this  work. 
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II.  CLASSIFICATIONS  OF  MICROINSTABILITIES 
There  are  several  type  of  microinstabilities  which  may  be 
present  in  an  MPD  plasma.  The  ones  which  are  most  likely  to 
occur  are  those  classified  as  cross-field  microinstabilities  [4]. 
The  name  cross-field  is  derived  from  the  energy  source  for  the 
instabilities,  namely  the  current  arc,  which  is  generally 
perpendicular  to  the  magnetic  field  in  a  self-field  MPD  thruster. 
The  cross-field  microinstabilities  fall  under  the  broader 
category  of  streaming  instabilities.  Streaming  instabilities  are 
caused  either  by  a  beam  of  energetic  particles  travelling  through 
the  plasma,  or  by  a  current  driven  through  the  plasma;  both 
result  in  different  species  suffering  drifts  relative  to  each 
other.  The  drift  energy  (or  current)  provides  the  source  for  the 
instability.  Experimental  measurements  have  shown  that  the 
electron  drift  velocity,  Ujg  is  typically  1-100  times  larger  than 
the  ion  thermal  velocity,  v.,  [31].  On  this  basis,  it  was 
assumed  that  current  driven  instabilities  are  the  dominant 
instability  in  MPD  thrusters  [31].  The  cross-field 
microinstabilities  which  are  most  likely  to  be  observed  are  the 
electron-acoustic  instability  (EAI)  (not  to  be  confused  with  the 
macroscopic  electron-acoustic  instability  which  has  been 
postulated  by  Rerapfer,  et.  al .  to  play  a  role  in  onset 
phenomena [ 10] ) ,  ion-acoustic  instability  ( IAI ) ,  electron- 
cyclotron  drift  instability  (ECDI),  drift-cyclotron  instability 
( DCI ) ,  generalized  lower  hybrid  drift  instability  (GLHDI ) ,  and 
the  two-stream  instability  (TSI)  (otherwise  known  as  Bunemann  or 
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hydrodynamic  instability). 


II. 1  DESCRIPTION  OF  RELEVANT  MICRCINS7ASILITIES 
II. 1.1  Electron  acoustic  instability  ( EAI ) 

The  electron  acoustic  instability  propagates  longitudinally. 
Its  direction  of  propagation  and  temperature  operation  regime  is 
given  by  [4]: 


r. 


k2  cos-0  v.J 

1  4.  _ 


(II. 1) 


where  0  s  angle  of  propagation  with  respect  to  the  magnetic 
field,  <j).y  2  lower  hybrid  frequency,  k  s  wave  number,  and  T.  m 
ion  temperature.  In  the  limit  of  6  =  n/2  (propagation 
perpendicular  to  the  magnetic  field)  this  equation  simplifies  to 
T; /Ta  >>  1.  If  the  wave  is  propagating  at  an  angle  oblique  to 
the  magnetic  field,  the  T./Tg  value  required  for  the  existence  of 
the  instability  must  be  even  larger.  Although  it  is  possible 
that  this  mode  could  be  present  in  extreme  cases,  it  is  not  very 
likely  to  occur  because  in  the  MPD  thruster  Tg  is  on  the  order  of 
T  .  In  any  event,  this  instability  will  not  exist  unless  the  dot 
product  of  the  wave  number  and  the  electron  drift  velocity  is 
less  than  zero[4,17]. 

k  •  0m  <  0  (II. 2) 

In  other  words,  the  instability  will  only  occur  if  the  electrons 
flow  in  the  opposite  direction  of  the  wave. 

Examining  figure  II. 1  to  see  the  likely  flow  pattern  of 
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current  in  the  MPD  thruster,  one  can  see  that  the  for  the 
instability  to  exist,  the  wave  must  have  a  component  in  the 
direction  of  the  contour  lines  going  from  anode  to  cathode.  As 
the  instability  propagates  in  a  longitudinal,  or  acoustic,  manner 
it  can  only  propagate  upstream  at  subsonic  flow  velocities. 
Because  the  ion  drift  velocity  is  hypersonic  in  the  MPD  thruster, 
the  EAI  is  only  expected  in  the  inter-electrode  space  and  near 
the  anode.  In  any  event,  this  instability  is  not  likely  to  be 
observed  because  it  is  damped  in  most  experimental  situations 
C17]  . 


II. 1.2  Ion  Acoustic  Instability  ( IAI ) 

The  ion-acoustic  instability  (IAI)  occurs  if  [17] 


KT.: 

~ktI 


<  l+k2rL2 


(II. 3) 


where  r s  =  electron  Larmor  radius,  and 


> 


a 


(II. 4) 


where  v.,  =  ion  sonic  velocity,  k,/k  =  cosine  of  angle  of  wave 
propagation  with  respect  to  magnetic  field,  Mj  a  mass  of  ion,  vtg 
2  electron  thermal  velocity,  and  u>pg  a  electron  plasma  frequency. 
In  the  limit  of  an  infinitely  strong  magnetic  field  (r^  =  0)  the 
thermal  criteria  for  existence  of  the  instability  reduces  to 
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T./T3  <<  1.  If  k./k  is  as  small  as  )•  (nearly  perpendicular 

to  the  magnetic  field)  the  condition  for  existence  of  the 
microinstability  if  lhg  £  (1-2)  v;.  in  the  whole  range  T.  >  T 
[17].  This  instability  may  occur  in  the  MPD  thruster  although 
the  thruster  is  considered  largely  isothermal  because  of  an 
extreme  sensitivity  to  the  actual  electron  and  ion  temperature 
distributions.  Due  to  lack  of  dependence  on  the  direction  of  the 
drift  velocities,  theoretically  it  could  occur  anywhere  in  the 
thruster.  Higher  electron  velocities  are  expected  near  the 
cathode,  though,  which  makes  the  cathode  the  most  likely  region 
of  existence  of  the  IAI .  However,  like  the  electron-acoustic 
mode,  the  ion-acoustic  mode  is  damped  for  most  experimental 
situations,  so  it  is  unlikely  to  be  seen  in  MPD  plasmas. 

11.1.3  Electron  Cyclotron  Drift  Instability  (ECDI) 

The  electron-cyclotron  drift  instability  (ECDI)  operates  at 
discrete  frequencies  near  the  harmonics  of  the  electron  cyclotron 
frequency.  The  criteria  for  the  instability  to  exist  is  [18] 

k  <  Inl-j^r  (II. 5) 

I  UJe  i 

where  |n|  5  the  number  of  the  harmonic  of  the  electron  cyclotron 
frequency.  This  is  the  range  where  Bernstein  harmonics  have  a 
negative  energy.  A  Bernstein  wave  is  an  electrostatic  wave  which 
propagates  perpendicular  to  the  magnetic  field  at  harmonics  of 
the  cyclotron  frequencies,  and  there  are  both  ion  and  electron 
modes  of  this  wave.  When  this  negative  energy  wave  comes  into 
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contact  with  a  frequency  resonance  or  sink  to  which  it  can  lose 
its  energy,  an  instability  will  result  [18].  What  happens, 
looking  at  the  dispersion  relation,  is  that  the  branch 
corresponding  to  Bernstein  waves  intersects  the  ion  acoustic 
branch  when  the  energy  of  the  wave  is  negative.  Thus,  the 
instability  will  probably  appear  in  the  same  regions  as  the  IAI, 
namely  in  the  cathode  region.  Lashmore-Davies  gives  a  solution 
for  the  ECDI  which  is  valid  for  [18] 
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It  has  been  theorized  that  these  waves  may  be  connected  with 
anomalous  resistivity  [20].  Part  of  the  problem  associated  with 
this  instability  is  the  non-linear  mechanism  which  saturates  the 
instability.  It  turns  out  that  the  collective  friction  force 
produced  by  the  instabilities  does  not  always  stabilize  them 
(e.g.  for  >>  vti).  Vortical  instabilities  in  an  electron  gas 
and  anomalous  viscosity  could  be  some  means  of  saturation[20] . 


I I. 1.4  Drift  Cyclotron  Instability  (DCI) 

The  drift-cyclotron  instability  (DCI),  first  considered  by 
Mikhailovsky  and  Timofeev  [16],  is  caused  by  gradients  in  the 
plasma  properties,  and  operates  at  discrete  frequencies  near  the 
ion  cyclotron  frequency.  Because  it  is  gradient  driven,  one 
might  be  tempted  to  classify  this  microinstability  as  a  Rayleigh- 
Taylor  instability.  However,  the  instability  energy  source  is 
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the  cross-field  current,  not  a  non-electromagnet ic  body  force; 
therefore,  it  still  falls  under  the  category  of  streaming 
instabilities.  The  instability  exists  if  [23]; 

>  ;  (II. 7) 

V..  ~  \j  M.- 

where  (h.  s  ion  drift  velocity.  Thus,  this  instability  would 
most  likely  exist  where  the  ion  drift  velocity  is  relatively 
high,  such  as  near  the  thruster  exit  plane  and  out  into  the 
plume.  Finite  3  (ratio  of  plasma  pressure  to  magnetic  pressure) 
increases  the  frequency  and  decreases  the  growth  rate  but  does 
not  completely  stabilize  the  instability.  When 


where  r-.  =  ion  larmor  radius  and  2  length  scale  of  density 
gradient,  the  DCI  transitions  to  become  the  lower  hybrid  drift 
instability  ( LHDI ) ,  or  the  lower  hybrid  gradient  driven 
instability  ( LHGDI ) ,  discussed  in  the  next  section. 

II. 1.5  Generalized  Lower  Hybrid  Drift  Instability  (GLHDI) 

The  generalized  lower  hybrid  drift  instability  (GLHDI) 
actually  incorporates  other  microinstabilities  as  limiting  cases. 
This  instability  can  occur  at  relatively  low  Uje  and  operates 
near  the  lower  hybrid  frequency.  This  instability  is  thought  to 
be  the  primary  reason  for  anomalous  resistivity,  and  thus 
turbulent  plasma  heating.  The  general  solution  does  not  have  an 
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analytic  representation  for  the  range  of  stability.  Therefore, 
the  limiting  cases  will  be  discussed. 

11. 1.5.1  Lower  Hybrid  Drift  Instability  (LHDI) 

The  lower  hybrid  drift  instability  (LHDI),  alternatively 
named  as  the  lower  hybrid  gradient  driven  instability  (LHGDI)  by 
Choueiri  [2,4]  occurs  when  gradients  in  the  plasma  parameters 
(inhomogeneities)  exist.  Numerically,  this  exists  as  kz  goes  to 
zero  and  \  ILg  -  l)Y  !  goes  to  zero  [24].  In  other  words,  the  waves 
in  this  instability  travel  in  a  "flute"  mode,  perpendicular  to 
the  local  magnetic  field  vector.  Thus,  in  the  MPD  thruster,  one 
would  expect  to  see  this  instability  near  the  cathode  where  the 
electrons  and  ions  move  roughly  in  the  same  direction.  The 
primary  gradients  responsible  for  this  instability  are  density, 
temperature,  and  magnetic  field  gradients.  The  instability 
itself  can  be  attributed  to  a  coupling  between  a  negative  energy 
drift  wave  and  positive  energy  lower  hybrid  wave  for  moderate 
gradients,  and  by  an  inverse  Landau  damping  caused  by  a.  resonance 
between  the  drift  wave  and  ions  which  are  near  the  drift  wave 
velocity  for  strong  gradients  [4].  This  instability's  length 
scale  is  on  the  order  of  an  electron  larmor  radius  [4], 

11. 1.5. 2  Modified  two-stream  instability  (MTS1) 

The  modified  two-stream  instability  (MTSI),  alternatively 
named  as  the  lower  hybrid  current  driven  instability  (LHCD1)  by 
Choueiri  [2,4],  is  the  instability  which  results  due  to  a 
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relative  streaming  of  charged  particles  due  to  the  E  x  B  drift  (E 
s  electric  field  vector)  or  the  presence  of  charged  beams  being 
the  major  source  of  free  energy.  This  is  also  referred  to  as  the 
kinetic  cross-field  streaming  instability  by  Wu  et.  al  [25]. 
Numerically,  this  case  is  represented  as  d(ln  n)/dx  goes  to  zero 
and  d  (In  B)/dx  goes  to  zero,  while  UJg  >>v...  In  the 
collisionless  limit,  this  wave  propagates  at  small  but  finite 
angles  from  the  local  perpendicular  to  the  local  magnetic  field 
vector.  In  the  MPD  thruster,  this  translates  into  a  wave  which 
travels  either  radially  or  axially,  but  has  a  slight  azimuthal 
component.  The  azimuthal  component  is  on  the  order  of  the  square 
root  of  the  mass  ratio  times  the  perpendicular  components.  It 
most  likely  occurs  near  the  anode  region  where  density  gradients 
and  magnetic  field  gradients  are  small,  but  electron  drift 
velocities  are  fairly  high. 

This  instability  does  not  depend  on  resonance  between  the 
wave  and  a  small  number  of  particles  in  the  velocity  distribution 
for  the  electrostatic  limit,  and  it  becomes  more  kinetic 
(resonant)  in  nature  as  electromagnetic  effects  becomes  important 
(when  Alfv6n  velocity  is  exceeded) [25] .  The  instability  results 
from  a  coupling  of  the  negative  energy  beam  of  particles  and  a 
positive  energy  lower  hybrid  wave.  As  electromagnetic  effects 
become  more  important,  a  coupling  with  an  electromagnetic  mode, 
such  as  a  whistler  mode,  may  become  significant [25] .  Although  it 
is  damped  by  electromagnetic  waves  when  operating  below  the 
Alfv6n  velocity,  it  is  not  stabilized  when  the  drifts  are  above 
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this  level. 

Lastly,  the  instability's  length  scale  is  larger  than  the 
electron  Larraor  radius.  This  means  that  it  is  possible  to 
measure  the  lower  hybrid  oscillations  of  this  mode  by  a  probe 
with  dimensions  smaller  than  an  electron  Lartnor  radius. 

11.1.6  Two-stream  instability  (TSI) 

The  two-stream  instability  (TSI),  or  Bunemann  instability, 
is  actually  a  limiting  case  of  the  MTSI,  but  as  it  lends  itself 
to  discussion  based  on  hydrodynamic  as  opposed  to  kinetic  theory, 
it  is  discussed  separately.  It  is  the  modified  two-stream 
instability  in  unmagnetized  plasmas  or  plasmas  where  the  current 
is  aligned  with  the  magnetic  field.  For  small  0,  if  kz  is 
sufficiently  large,  a  uniform  magnetized  plasma  with  a  cross 
field  drift  will  reduce  to  the  unmagnetized  case.  In  this  case, 
the  instability  will  operate  near  the  electron  plasma  frequency 
if  kr  g>>l  and  lbe>v.g.  Although  this  instability  is  not  likely  to 
be  present  in  the  MPD  thruster  due  to  the  large  azimuthal 
magnetic  fields,  its  relative  simplicity  allows  one  to  get  a 
basic  understanding  of  more  complex  microinstabilities.  This 
instability  is  used  to  demonstrate  the  technique  of  linear 
analysis  to  determine  dispersion  relations  in  the  next  section. 

11.2  Linearized  derivation  of  dispersion  relation 
Considering  a  uniform  plasma  in  which  the  ions  are 
stationary  and  the  electrons  have  drift  velocity  Vg  relative  to 
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the  ions,  a  dispersion  relation  may  be  deduced  under  the 


assumptions  of  a  cold  plasma  ( KT  =  0  for  both  species;  and  no 


magnetic  field.  Starting  with  the  equations  of  motion. 
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for  both  ions  and  electrons,  where  n^  s  equilibrium  number 
density,  vg  =  velocity  vector  of  species  s,  n.  =  number  density 
of  species  s,  and  E  s  electric  field  vector.  Letting  density, 
electric  field  and  velocity  of  each  species  be  represented  by  a 
sum  of  perturbation  and  equilibrium  values,  the  gradient  term  in 
the  ion  equation  is  dropped  since  the  ions  are  stationary.  Also 
neglecting  one  of  the  two  gradient  terms  because  the  equilibrium 
velocity,  v-, ,  is  uniform  yields  the  linearized  equations: 

=  en0E1  (II.  10) 
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where  =  first  order  ion  velocity  perturbation,  Ej  s  first 
order  electric  field  perturbation,  and  vej  s  first  order  electron 
velocity  perturbation.  Looking  for  electrostatic  waves  of  the 
form  el'kx'^°  ,  where  o>  s  wave  frequency  and  i  s  \f~T  , 

V=ik  (11.12) 

Next,  it  is  necessary  to  solve  for  both  species  density 
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perturbations  and  both  species  velocity  perturbations.  Since  the 
instability  causes  high  frequency  plasma  oscillations,  the  plasma 
approximation,  that  is  n.  -  n,,  cannot  be  used,  which  means 
Poisson's  equation  must  be  used  instead: 

(11.13) 

where  e-,  =  permitivitty  of  free  space  and  n.  ==  ion  number 
density.  This  equation  must  be  linearized.  Substituting  the 
solutions  for  the  densities  and  velocities  evaluated  earlier 
yields  the  dispersion  relation. 

_L. - i — J  (ii.i4> 

‘V.  0 >"  -u)-)cvv)-j 

l  J 

where  u>p  =  overall  plasma  frequency. 

If  the  dispersion  relation  is  multiplied  through  to  get  a 
common  denominator,  a  fourth  order  polynomial  for  c o  results. 
Theoretically  both  k  and  u>  could  be  complex.  However,  to  find  if 
a  plasma  is  unstable  one  must  first  look  at  the  situation.  If 
the  plasma  merely  exists,  disturbance  waves  must  be  considered 
across  a  whole  spectrum  due  to  the  randomness  of  nature;  a 
specific  disturbance  cannot  be  readily  predicted.  Thus  k  is 
typically  treated  as  an  independent,  real  parameter,  while  o>  is 
allowed  to  be  complex.  Solving  the  above  equation,  there  can 
exist  both  real  and  imaginary  roots  of  u  .  If  all  the  roots  are 
real,  no  instability  exists.  However,  if  complex  roots  are 
present,  they  exist  as  complex  conjugates.  The  conjugate  with  a 
negative  imaginary  term  will  indicate  an  exponentially  damped 
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wave.  The  conjugate  with  a  positive  imaginary  term  will  have  an 
exponential  growth.  These  two  effects  do  not  cancel  out  when 
added;  all  that  is  seen  will  be  exponential  growth.  Thus,  if 
complex  roots  exist,  the  plasma  is  unstable.  As  kv.  decreases 
there  exists  a  point  where  it  finally  yields  complex  roots.  What 
this  translates  to  is  that  for  the  Bunemann  instability,  for 
sufficiently  small  kv- ,  the  plasma  is  unstable.  This  means  that 
since  wavelength  is  inversely  proportional  to  k,  for  any  given 
v.j,  the  plasma  is  unstable  to  long  wavelength  oscillations.  If 
the  plasma  has  a  finite  temperature,  the  fluid  equations  cannot 
be  used  to  predict  the  instability;  kinetic  theory,  however,  can 
rake  into  account  finite  temperature  plasmas,  although  the 
analysis  is  also  more  complex.  Kinetic  theory  predicts  that 
Landau  damping  will  occur  for  v-<  <  thermal  velocity,  so  no 
instability  will  exist  if  v-j  is  too  small. 

This  instability  can  be  physically  explained  as  follows. 

Both  electrons  and  ions  have  their  own  natural  frequencies.  What 
happens  is  that  for  the  right  value  of  kvj,  these  frequencies 
coincide  due  to  the  doppler  shift  of  the  moving  electron  plasma 
fluctuations.  Also,  these  electron  oscillations  have  negative 
energy.  In  other  words,  the  electrons  have  less  kinetic  energy 
when  the  oscillations  occur  than  they  did  when  they  were  absent. 
In  contrast,  the  ions  have  positive  energy,  meaning  they  have 
more  energy  when  the  oscillations  are  present  than  when  they  are 
absent.  What  this  means  is  that  both  waves  can  continue  to  grow 
while  conservation  of  energy  is  satisfied  because  the  total 
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energy  remains  constant. 

As  an  aside,  when  a  plasma  is  excited  by  a  specific 
wavelength  disturbance,  say  by  a  microwave  beam,  oj  stays  real, 
and  k  is  considered  to  be  complex.  Thus  if  the  imaginary  part  of 
k  is  positive,  the  wave  will  be  damped,  while  if  the  imaginary 
part  is  negative,  exponential  growth  will  occur. 

II. 3  Relationships  between  Microinstabilities 

Now  that  a  foundation  for  a  physical  and  mathematical 
understanding  of  what  a  microinstability  has  been  laid, 
relationships  may  be  explored.  The  GLHDI  can  be  viewed  as  the 
most  basic  instability,  where  most  of  the  other  instabilities  are 
special  cases,  depending  on  limits  of  key  parameters.  In  fact, 
the  most  frequently  cited  frequencies  of  microinstabilities  in 
the  MPD  thruster  has  been  at  the  lower  hybrid  frequency  [1,30]. 
The  EAI  may  be  viewed  as  the  limit  of  the  GLHDI  for  T./T.  >>  1, 
while  the  IAI  is  the  case  where  T, /Te  <<  1.  The  LHCDI  or  MTSI  is 
the  instability  which  results  when  current  is  the  primary  energy 
source  for  the  instability,  where  the  TSI  is  the  limit  of  this 
instability  in  the  case  of  a  cold  plasma  with  no  magnetic  field. 
The  LHGDI  or  LHDI  is  the  instability  which  results  when  gradients 
in  plasma  properties  are  the  energy  source  for  the  instability. 

The  only  instabilities  which  are  not  exactly  limits  of  the 


GLHDI  are  the  cyclotron  instabilities.  However,  they  are 
related.  Although  they  operate  at  a  different  characteristic 
frequency  (either  the  ion  or  electron  cyclotron  frequency  rather 
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than  the  lower  hybrid  frequency),  limits  of  the  GLHDI  are  linked 
to  the  cyclotron  instabilities.  The  dispersion  relation  of  the 
ECDI  intersects  the  ion-acoustic  branch  of  the  GLHDI  when 
equation  (1 1. 5)  is  met,  and  the  DCI  actually  becomes  the  LHGDI 
when  equation  (1 1. 8)  is  met.  The  main  difference  between  all  of 
these  instabilities,  with  the  exception  of  the  ECDI,  is  the 
conditions  of  the  plasma  which  cause  them;  they  can  all  be 
analyzed  in  the  same  manner  by  using  the  linearization  method  to 
find  the  dispersion  relation  as  demonstrated  in  the  section  II. 2, 
although  particle  kinetics  rather  than  hydrodynamic  theory  must 
be  used. 

Figure  11.2  gives  a  graphical  depiction  of  the  relationships 
between  microinstabilities.  The  cube  represents  all  lower  hybrid 
instabilities  (the  GLHDI),  the  dots  represent  the  electron- 
cyclotron  drift  instability,  and  the  cylinder  represents  the 
drift-cyclotron  instability.  All  of  these  oscillations  can 
simultaneously  occur  in  the  same  parameter  space,  but  not 
necessarily  at  the  same  spatial  location.  The  colored  panels  of 
the  cube  represent  limiting  cases  of  the  GLHDI  for  the  plasma 
parameters  listed  on  the  three  axes.  The  location  of  the  colored 
sections  should  give  the  reader  an  idea  under  what  conditions 
that  particular  microinstability  can  occur.  It  may  be  noticed 
that  the  dots  and  circles  are  denser  in  some  portions  of  the  cube 
than  others.  This  is  because  theoretically,  these  instabilities 
can  occur  anywhere  in  the  parameter  space,  but  are  much  more 
likely  to  occur  where  the  circles  or  dots  are  thickest.  Lastly, 
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when  a  colored  surface  shows  a  brighter  color  where  one  of  the 
circles  intersects  that  surface,  it  indicates  a  linkage  between 
the  two  raicroinstabi 1 ities .  As  the  DCI  is  known  to  transition  to 
the  LHDI  for  large  gradients,  the  intersection  between  the 
circles  representing  the  DCI  and  the  plane  representing  the  LHDI 
is  colored  a  brighter  green  than  the  rest  of  the  LHDI  plane. 
However,  the  ECDI  has  no  direct  linkage  to  the  GLHDI  or  DCI  in 
this  parameter  space. 


Figure  II.l  -  Enclosed  current  contours  for  operation  at  6.8  kA 
and  2  g/sec  mass  flow 
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Figure  II. 2  -  Pictorial  depiction  of  microinstabilities  existing 
in  parameter  space. 


III.  Experimental  Work 


Experiments  were  performed  on  a  MW  level,  quasi -steady  self- 
field  MPD  thruster  with  a  1  millisecond  pulse  duration  at  the 
Electric  Proculsion  Laboratory  at  the  Philips  Laboratory  at 
Edwards  Air  Force  Base  in  California.  A  quasi-steady  device 
produces  a  current  pulse  which  achieves  a  stable,  constant  value 
of  current  for  a  long  enough  duration  to  approximate  steady-state 
behavior.  The  global  operating  characteristics  of  the  thruster, 
the  azimuthal  magnetic  field,  number  density  and  electron 
temperature  were  measured.  A  complete  description  of  equipment 
used  to  take  measurements  as  well  as  information  about  the 
thruster,  and  the  chamber  in  which  the  thruster  was  operated  is 
presented  in  Appendix  A. 


III.l  Data  Acquisition 

I I I. 1.1  Characterization  of  thruster  performance 

In  order  to  acquire  the  desired  properties,  it  was  first 
necessary  to  characterize  the  performance  of  the  thruster.  The 
thruster  is  fired  by  first  injecting  propellant  into  the 
thrustar,  and  subsequently  initiating  a  high-voltage  spark  in 
order  to  ionize  the  propellant.  This  ionized  propellant  provides 
a  conducting  path  for  the  high  current  discharge  which  produces 
thrust.  The  high  current  discharge  was  supplied  by  a  10  section 
LC  (Inductor-Capacitor)  pulse  forming  network  (PFN),  which  is 
capable  of  producing  a  1  millisecond  current  pulse  of  up  to  40  kA 
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and  400  V  to  the  MPD  thruster  assuming  a  perfectly  matched  load 
[34].  As  the  gas  pulse  and  spark  determine  when  and  how  the 
thruster  fires,  the  first  tests  conducted  were  the  propellant  gas 
pulse  and  spark  timing  measurements.  A  plot  of  a  typical  gas 
pulse,  with  the  spark  trigger  shown  as  a  noise  spike,  is  shown  in 
figure  III.l. 

The  gas  pulse  had  an  exponential  rise  and  decay  portion  with 
a  relatively  flat  steady  state  portion  was  reached  in 
approximately  20  milliseconds  after  triggering  the  valve  and 
lasted  approximately  65  milliseconds  .  The  gas  pulse  was  marred 
by  a  small  dip  in  the  quasi-steady  mass  flow  shortly  after  the 
quasi-steady  state  had  been  reached,  but  this  variation  was  less 
than  5%  of  the  pulse  magnitude  and  does  not  appear  to  have  much 
of  an  effect  on  the  thruster  operation. 

The  finite  rise  time  of  the  gas  pulse  does  not  cause  many 
complications  in  a  quasi-steady  firing.  However,  a  20 
millisecond  delay  in  the  initiation  of  the  spark  trigger  was 
added  to  ensure  that  the  engine  fired  while  in  the  steady  state 
portion  of  the  gas  pulse. 

Before  testing,  it  was  believed  that  raicroinstabilities 
would  be  more  noticeable  if  the  MPD  thruster  were  operated  at  a 
point  near,  but  below  onset.  In  order  to  determine  the  input 
parameters  (power  supply  voltage  and  current)  corresponding  to  a 
point  slightly  below  onset,  measurements  of  the  global  thruster 
voltage  and  current  were  conducted.  The  signals  were  acquired  on 
the  DSA-601  digital  oscilloscope,  acquiring  the  voltage  and 


35 


current  simultaneously.  Figure  I I I. 2  shows  sample  voltage  and 
current  traces  at  levels  before  and  after  onset. 

The  voltage  trace  was  characterized  by  a  large  noise  spike 
at  the  trigger  time,  followed  by  a  steady  state  region  of 
approximately  1.4  msec.  It  then  decayed  to  a  second  or  even 
third  voltage  level  depending  on  the  particular  firing.  In  the 
steady  state  region,  noise  was  very  small  for  low  power  settings. 
When  the  PFN  voltage  was  in  the  range  350-550  V,  large  voltage 
fluctuations  were  observed,  on  the  order  of  10  per  cent  of  steady 
state  voltage  or  larger.  Malliaris,  et.  al  has  indicated  that 
these  large  voltage  fluctuations  indicate  onset  [9].  These 
fluctuations  occurred  at  a  current  level  of  about  8  kA,  with 
widely  varying  voltages.  Figure  III. 3  shows  a  plot  of  global 
voltage  vs.  current,  referred  to  hereafter  as  a  V-I  curve.  In 
this  plot,  a  near  infinite  slope  was  present  at  a  current  of  8 
kA,  leading  to  the  observation  that  onset  can  be  also  defined  as 
the  region  of  discontinuity  in  a  V-I  curve.  Above  550  V,  the 
voltage  fluctuations  drop  below  10  per  cent  of  the  steady  state 
value.  Although  the  oscillations  are  greater  than  at  low  power 
levels,  it  appears  that  the  instabilities  seen  earlier  have 
saturated.  However,  at  these  power  levels,  erosion  seems  to 
become  significant. 

Current  responses  were  relatively  free  of  noise  and  had  an 
initial  overshoot  of  no  more  than  10  per  cent.  The  steady  state 
portion  was  extremely  flat,  especially  at  lower  power  levels,  and 
the  current  pulse  itself  had  exponential  growth  and  decay.  The 
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current  traces  were  almost  identical  at  all  voltages  in  the  onset 
regime.  However,  two  forms  of  traces  occurred  above  the  onset 
regime.  In  the  majority  of  cases,  the  current  increased  by  a 
steady,  but  slight  amount  over  the  steady  state  portion. 

However,  in  a  few  cases,  a  large  peak  occurred  in  the  steady 
state  region  just  before  the  current  decayed.  This  corresponds 
to  the  voltage  drop  seen  in  the  voltage  readings  on  the  same 
firings.  The  most  likely  explanation  is  that  erosion  caused 
additional  mass  from  the  copper  electrodes,  boron-nitride  from 
the  injection  plate,  etc.,  to  enter  the  thruster,  increasing  the 
number  of  charge  carriers,  and  consequently,  current.  If  the 
tests  with  these  types  of  curves  are  plotted  on  a  V-I  curve, 
there  is  a  clear  separation  of  data  from  the  majority  of  firings. 
This  set  of  data  has  a  smaller  slope  than  the  other  cluster  of 
points,  seeming  to  indicate  a  regime  of  higher  performance.  This 
makes  sense  if  the  above  explanation  is  accepted;  the  additional 
mass  of  eroded  materials  is  expelled,  increasing  thrust  at  a 
given  power  level.  Obviously,  though,  this  type  of  operation  is 
not  desired  because  parts  of  the  thruster  erode. 

Tests  were  also  made  at  1.0  and  0.5  g/sec  up  to  the  onset 
condition.  Based  on  the  V-I  curves,  an  operation  point  with  a 
current  level  of  about  6.8  kA  at  a  mass  flow  of  2  g/sec  was 
chosen  for  the  remaining  experiments.  These  conditions 
corresponded  to  a  point  slightly  below  onset. 
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III. 1.2  Testing  Performed 

Further  testing  was  designed  to  produce  spatial  maps  of 
magnetic  field  strengths  (azimuthal  direction),  electron  number 
density,  electron  temperature,  and  the  frequency  of  number 
density  fluctuations.  Measurements  were  desired  both  in  the  near 
plume  and  inside  the  thruster,  so  a  22  point  test  pattern 
involving  locations  in  the  thruster  plume  and  inside  the  thruster 
was  used.  This  test  pattern,  which  is  shown  in  Figure  III. 4,  was 
formulated  assuming  that  the  plasma  discharge  would  be  axially 
symmetric;  only  axial  and  radial  positions  were  varied. 

Measurements  of  the  magnetic  field  were  then  taken  at  each 
of  these  grid  locations  by  magnetic  probes.  Only  between  one  and 
three  measurements  were  able  to  be  taken  at  each  point  due  to 
time  constraints,  with  one  point  being  retested  after  the  entire 
sequence  to  check  hysteresis.  Magnetic  field  measurements  were 
desired  because  the  magnetic  field  strength  is  required  to 
calculate  the  cyclotron  frequencies  and  hybrid  frequencies. 

After  magnetic  probe  measurements  were  concluded, 
measurements  with  the  Langmuir  triple  probe  were  conducted.  Due 
to  the  possibility  of  probe  contamination,  only  one  firing  was 
made  at  each  grid  location,  except  at  positions  which  produced 
anomalous  signal  traces;  these. points  were  retested.  To  further 
reduce  contamination,  the  probes  were  cleaned  by  glow  discharge 
after  every  5  firings.  At  each  grid  location,  the  probe's  mean 
voltage  and  current,  fluctuations  in  the  probe's  voltage  and 
current,  and  global  voltage  and  current  were  measured 
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simultaneously.  Electron  number  density  and  electron  temperature 
measurements  were  calculated  based  on  these  measurements.  These 
two  properties  were  desired  because  they  characterize  the  plasma 
state,  while  the  plasma  fluctuations  w  -re  measured  to  determine 
if  the  plasma  oscillations  corresponded  to  a  characteristic 
frequency  of  suspected  microinstabilities.  Unfortunately,  after 
the  first  test  sequence  was  completed,  the  facility  became 
unavailable  for  further  tests. 

III. 1.3  Probes 

1 1 1. 1.3.1  Magnetic  Probes: 

Two  different  types  of  magnetic  probes  were  used  to 
measure  magnetic  fields,  although  only  one  type  worked 
successfully.  The  first  type  was  a  Hall  probe.  The  theory 
governing  the  operation  of  this  probe  is  explained  in  Appendix  B. 
This  was  the  preferred  type  of  probe  to  be  u?  1  because  of  the 
linear  relation  between  magnetic  field  strength  and  probe 
voltage,  and  it  has  been  used  with  a  great  deal  of  success  in 
steady  state  MPD  devices.  Unfortunately,  it  did  not  work  well  in 
the  quasi-steady  MPD  device  used  in  this  work.  One  disadvantage 
of  the  Hall  probe  is  that  the  output  voltage  is  on  the  order  of  a 
few  microvolts.  Extremely  large  noise  spikes  at  start-up  of  the 
quasi-steady  thruster  drowned  out  any  signal  the  Hall  probe 
produced.  This  phenomenon  will  be  discussed  further  in  section 

1 1 1. 2.1  and  in  great  detail  in  Appendix  B. 

In  order  to  get  a  discernible  signal,  an  induction  type 
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probe  with  an  analog  integrator  circuit  was  used.  The  induction 
probe  itself  produces  a  voltage  which  is  proportional  to  the 
change  in  magnetic  field  with  time.  Thus,  to  get  a  voltage  which 
is  proportional  to  the  magnetic  field  strength,  the  signal  must 
be  integrated.  This  type  of  probe  is  a  little  more  complex  on 
account  of  the  presence  of  an  integrator  circuit,  but  much  higher 
signal-to-noise  ratios  can  be  achieved.  This  type  of  device  has 
been  used  almost  exclusively  in  pulsed  devices  up  to  this  time, 
and  it  was  used  successfully  in  these  experiments.  The 
construction  and  use  of  the  induction  probe  is  discussed  in 
Appendix  B. 

1 1 1. 1.3. 2  Langmuir  Probes: 

In  order  to  measure  number  density  and  electron  temperature, 
the  triple  Langmuir  probe  was  used.  This  probe  can  measure 
voltage  and  current  simultaneously,  which  allows  simultaneous 
measurements  of  electron  temperature  and  number  density  based  on 
the  application  of  the  probe  theory.  In  these  calculations, 

Laf ramboise ' s  method  with  the  Peterson-Talbot  curve  fits  [35,36], 
and  Bohm's  thin  sheath  criteria  [35,37]  were  used  to  calculate 
electron  number  density  and  electron  temperature.  A  more 
detailed  discussion  of  probe  theory  and  construction  is  given  in 
Appendix  C. 

In  order  to  get  mean  values  of  plasma  properties  at  the  time 
in  the  pulse  when  the  frequency  data  of  plasma  fluctuations  were 
acquired,  two  separate  oscilloscopes  were  used.  The  high 
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frequency  oscilloscope  was  triggered  about  halfway  through  the 
steady  state  portion  of  the  current  pulse.  Because  the  probe 
current  is  directly  proportional  to  electron  number  density  using 
thin-sheath  probe  theory  [35].  the  frequency  spectra  of  the 
current  oscillations  can  be  assumed  to  be  identical  to  the 
frequency  spectra  of  the  number  density  oscillations.  Likewise, 
the  frequency  spectra  of  probe  voltage  oscillations  can  be 
assumed  to  be  identical  to  the  frequency  spectra  of  electron 
temperature  oscillations  because  voltage  is  directly  proportional 
to  the  electron  temperature  in  the  thin  sheath  limit  [35]. 
However,  the  absolute  magnitudes  of  the  oscillations  for  both 
plasna  properties  are  merely  proportional  to  the  measured  voltage 
and  current,  and  depend  on  the  mean  properties  of  the  plasma. 

Fast  fourier  transforms  were  used  to  convert  voltage  and  current 
oscillations  into  frequency  spectra.  Although  the  frequency 
oscillations  without  wave  vector  and  amplitude  distribution  is 
not  sufficient  to  identify  the  types  of  instabilities  [15],  they 
do  give  an  indication  of  a  range  of  instabilities  which  may  be 
occurring . 

Probe  contamination  became  a  problem  by  causing  apparently 
erroneous  mean  values  to  be  calculated  inside  the  thruster. 
Cleaning  the  Langmuir  probes  by  means  of  a  glow  discharge  helped 
to  some  extent,  but  for  measurements  inside  the  thruster, 

Langmuir  probes  will  have  to  be  replaced  after  a  few  firings  in 


the  future. 
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III. 1.3  Calibration: 

Calibration  of  the  magnetic  probes  was  accomplished  by 
placing  the  magnetic  probe  at  a  radius  just  inside  the  anode  and 
very  close  to  the  back  plate  of  the  thruster,  and  measuring 
global  current  while  firing  the  thruster.  It  was  assumed  that 
all  of  the  measured  current  would  pass  through  a  circle  enclosed 
by  the  radius  the  probe  was  placed  at.  As  a  result,  the  actual 
magnetic  field  was  determined  based  on 

B.  =  -i^£  (III  - 1 ) 

6  2  nR 

where  s  azimuthal  magnetic  field,  and  R  s  radius  the  magnetic, 
field  is  measured  at.  The  voltage  of  the  probe  was  then  acquired 
after  passing  through  the  integrator.  The  thruster  current  was 
varied  and  readings  taken  at  each  of  these  points  to  produce  a 
calibration  curve.  A  more  detailed  discussion  of  calibration  is 
given  in  Appendix  B 

No  calibration  of  the  triple  Langmuir  probe  was  required, 
although  the  frequency  response  of  the  probe  was  checked. 

However,  this  was  actually  a  check  of  the  electronics  of  the 
probe  circuit  more  than  a  check  of  the  probes  themselves,  as  the 
operational  amplifiers  used  to  amplify  the  signal  were  the 
limiting  factor,  not  the  probe  itself. 

I I I. 2  Measurements 

III. 2.1.  Magnetic  Field 

Initial  attempts  to  make  magnetic  field  measurements 
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employed  a  Hall  probe.  However,  due  to  very  low  signal-to-noise 
ratios  (SNR)  which  made  results  indiscernible,  an  induction  type 
probe  with  an  integrator  circuit  was  used  instead  to  measure  the 
magnetic  field.  After  some  initial  calibration  problems,  usable 
signals  from  this  type  of  probe  were  acquired  because  the  signal 
levels  on  this  type  of  probe  were  a  few  orders  of  magnitude  above 
the  signals  produced  by  the  hall  probe.  However,  the  initial 
start-up  spike  still  caused  some  problems.  A  sample  magnetic 
field  trace  is  shown  in  Figure  III. 5. 

This  noise  appears  in  the  start  up  and  dissipates  fairly 
rapidly  with  time.  By  the  time  the  current  pulse  begins  its 
downward  swing,  this  noise  has  largely  disappeared.  Therefore, 
when  determining  the  voltage  level  of  the  probe,  the  voltage 
level  after  the  thruster  had  finished  firing  was  defined  as  the 
reference  level,  rather  than  the  voltage  level  before  the 
thruster  fired.  This  technique  of  measurement  seemed  to  yield 
consistent  results,  within  the  range  of  noise  still  present.  A 
fairly  steady  variation  of  magnetic  field  strength  with  axial  and 
radial  positions  resulted.  In  contrast,  when  the  voltage  before 
firing  was  defined  as  the  reference  level,  the  magnetic  field 
strengths  fluctuated  wildly  with  position.  This  was  not 
consistent  with  results  observed  by  other  experimenters  [12,38- 
41]. 

Only  one  magnetic  field  measurement  was  made  at  most 
locations  due  to  time  constraints.  The  voltage  trace  of  the 
probe  was  expected  to  look  similar  to  the  global  current  trace 


43 


because  of  the  relationship  between  magnetic  field  and  current, 
defined  in  equation  (III.l).  If  the  voltage  trace  behaved  as 
expected,  the  probe  was  moved  to  the  next  location.  However,  if 
unexpected  voltage  traces  occurred,  a  few  additional  measurements 
were  made  to  determine  if  the  unexpected  trace  was  repeatable. 

At  these  locations,  the  magnetic  fields  determined  from  the 
voltage  traces  are  averaged.  A  spatial  map  of  the  average 
magnetic  field  at  each  grid  location  is  shown  in  figure  III. 6. 
Based  on  the  magnetic  field  measurements  shown  in  figure  II I. 6, 
an  enclosed  current  contcur  map  was  generated,  shown  in  figure 
II. 1.  This  map  was  used  in  chapter  II  to  determine  the  location 
of  certain  microinstabilities  because  the  contour  lines  may  be 
interpreted  as  actual  current  paths  inside  the  thruster. 

However,  the  primary  purpose  of  the  magnetic  field  measurements 
was  to  provide  data  for  the  calculation  of  some  characteristic 
frequencies  of  the  plasma.  These  frequencies  are  tabulated  in 
Table  III.l. 

It  should  be  noted  that  positions  (17)  and  (23)  in  the  test 
pattern  (shown  in  figure  I I I. 4)  are  actually  the  same  location. 

As  both  values  of  Bg  are  fairly  close,  it  appears  that  hysteresis 
is  negligible  within  the  consistency  of  thruster  firings.  The 
major  trend  noticeable  here  and  in  the  contour  maps  is  that 
increases  significantly  as  the  position  moves  towards  the  rear  of 
the  thruster,  with  the  sharpest  gradients  occurring  inside  the 
thruster  .  However,  the  region  of  strongest  magnetic  field  for  a 
given  axial  position  was  about  halfway  between  the  electrodes. 
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As  the  characteristic  frequencies  are  directly  proportional  to 
B^,  they  follow  the  same  general  trends.  The  characteristic 
frequencies  listed  in  Table  III.l  will  be  important  in  the  next 
section . 

II I. 2. 2. 2  Langmuir  Probe  Measurements 
III. 2. 2. 2.1  Number  Density  Fluctuations: 

The  triple  Langmuir  probe  was  used  to  measure  number  density 
fluctuations,  while  simultaneously  measuring  the  mean  electron 
temperature  and  electron  number  density.  However,  to  accomplish 
this,  an  additional  oscilloscope  was  required  to  record  signals  . 
on  two  different  time  scales.  Global  current  and  voltage  were 
also  monitored  at  the  same  time,  using  the  same  time  scale  used 
for  measurement  of  mean  plasma  properties.  Sample  probe  traces 
for  the  high  frequency  oscillations  are  shown  in  figure  III. 7, 
and  sample  traces  of  the  absolute  values  of  probe  current  and 
voltage  used  for  calculation  of  mean  plasma  properties  are  shown 
in  figure  1 1 1. 8.  Both  voltage  and  current  traces  from  the  probe 
were  expected  to  look  similar;  namely,  to  have  a  fast  initial 
transient,  a  steady  state  region,  and  a  slower  final  transient 
because  the  temperature  and  number  densities  will  increase  as  the 
thruster  fires,  and  return  to  zero  when  the  firing  is  over.  The 
probe  traces  behaved  according  to  expectations  in  the  plume 
region  and  were  only  slightly  different  inside  the  thruster  along 
the  cathode.  However,  the  last  six  locations  (near  the  anode  and 
inter-electrode  area)  examined  appeared  to  produce  questionable 
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data.  The  values  for  voltage  and  current  of  the  probes  no  longer 
seemed  to  have  a  steady  variation  with  position.  Instead  the 
values  were  a  great  deal  higher,  which  produced  unreasonable 
electron  temperature  and  electron  number  density  values  when  the 
probe  theory  was  applied  to  the  measured  mean  values.  In 
addition,  the  strength  of  fluctuations  seemed  to  increase 
significantly. 

Looking  at  the  fast-fourier  transforms  ( FFT )  of  the  current 
traces  at  all  locations,  there  is  a  large,  broad  band  peak 
centered  between  100  kHz  and  1.5  MHz  depending  on  the  grid 
location.  It  can  also  be  noted  that  the  magnitude  of  this  peak  , 
increases  significantly  at  the  locations  where  onset  occurred. 
However,  this  peak  is  of  relatively  low  frequency,  and  is  on  the 
order  of  the  frequency  expected  of  the  global  hydrodynamic 
instabilities  associated  with  onset  [10,15],  which  have  been 
postulated  to  be  either  macroscopic  electron  acoustic 
instabilities  [10],  Pierce  instabilities  [13,14],  or  electro¬ 
thermal  instabilities! 10] . 

However,  the  experiments  performed  in  this  paper  were  based 
upon  expectations  of  oscillations  at  the  lower  hybrid  frequency 
of  the  plasma,  determined  by  the  magnetic  field  strength  at  the 
same  location.  The  values  of  the  lower  hybrid  frequency  at  each 
location  are  compiled  in  table  III.1,  where  the  numbered 
locations  are  shown  in  figure  III. 4.  However,  a  clear  peak  in 
the  frequency  spectra  at  the  lower-hybrid  frequency  at  any 
location  in  the  thruster  was  not  observed.  There  might  have  been 
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an  indication  of  lower  hybrid  oscillations  at  one  point  ne^r  the 
cathode  at  the  exit  plane,  but  the  observed  peak  was  so  close  to 
the  noise  level  as  to  be  inconclusive. 

1 1 1. 2. 2. 2. 2  Mean  Plasma  Properties 

Mean  plasma  properties  were  obtained  by  acquiring  the 
absolute  values  of  the  current  and  voltage  of  the  probe  at  the 
time  the  plasma  fluctuations  were  measured.  A  larger  time  scale 
was  used  to  allow  average  values  of  voltage  and  current  to  be 
obtained  in  the  neighborhood  of  the  time  at  which  plasma 
fluctuations  were  recorded.  A  more  detailed  discussion  on  how 
plasma  properties  were  acquired  from  voltage  and  current  traces 
is  explained  in  Appendix  C.  Table  I I I. 2  shows  various  plasma 
properties  at  each  measured  location,  where  the  numbered 
locations  are  shown  in  figure  I I I. 4 

It  can  be  noted  that  the  last  six  locations  have  abnormally 
high  electron  temperatures,  with  position  (19)  having  an 
extremely  high  value.  The  reasons  for  these  abnormal  results 
will  be  discussed  in  section  III. 3.  As  stated  earlier,  points 
(17)  and  (23)  comprise  a  check  for  hysteresis,  giving  an 
indication  of  the  error  in  the  last  six  points.  Not  only  was 
there  a  significant  difference . between  the  values  at  points  (17) 
and  (23),  but  all  of  the  points  in  between  were  so  far  beyond  the 
range  of  expected  results  that  Laf ramboise ' s  theory  could  not  be 
used  to  calculate  Te  and  nfi.  Instead,  Bohn's  thin  sheath  theory 


was  used. 


Examining  spatial  variations  of  the  properties  for  the  good 
locations,  the  trends  generally  follow  what  would  be  expected.  A 


contour  map  of  electron  temperature  is  shown  in  Figure  1 1 1. 9. 
Examining  the  map,  there  appears  to  be  a  low  temperature  pocket 
starting  at  the  exit  plane  near  the  anode,  and  extending  radially 
outward.  In  general,  electron  temperature  increases  at  locations 
toward  the  rear  of  the  thruster,  and  increase  as  position 
approaches  the  anode.  The  sharp  gradients  in  Te  apparent  in 
Figure  III. 9  are  due  to  the  abnormally  high  results  acquired  in 
the  anode  region  inside  the  thruster. 

A  contour  map  of  electron  number  density  is  shown  in  Figure 
III. 10.  Examining  the  map,  the  number  densities  are  highest  near 
the  cathode,  with  fairly  steep  gradients  in  this  region  as  well. 

In  general,  there  is  a  strong  decrease  in  number  density  with 
radial  position,  and  a  weaker  increase  in  number  density  as  the 
position  approaches  the  rear  of  the  thruster. 

The  other  two  properties  in  the  Table  III. 2  are  shown 
primarily  for  illustrative  purposes.  The  electron  debye  length,  a 
was  calculated  according  to 


(III. 2) 


where  e0  is  the  permitivitty  of  free  space,  n  is  total  number 
density,  and  K  is  the  Boltzmann  constant.  This  is  an  important 
parameter  in  probe  theory,  and  typically  the  electron  debye 
length  is  assumed  to  be  about  the  same  as  the  total  debye  length. 
The  electron  collision  frequency  is  the  sum  of  the  electron-ion 
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and  electron-electron  collision  frequencies  using  the  formula 
[6]: 


„  n  In  A 

V  ,  =  7 - 

'  T  >■  - 

1 


(III. 3) 


where  In  A  is  the  natural  log  of  the  maximum  impact  parameter. 
In  using  this  formula,  it  is  assumed  that  n  is  approximately 
equal  to  n0  and  only  singly  charged  ions  are  present.  In 
addition,  collision  frequency  is  necessary  to  find  the  electron 
Hall  parameter.  Anomalous  conductivity  has  been  seen  to  scale 
strongly  with  the  electron  Hall  parameter  [4]. 


I I I. 3  Interpretation  of  Results 
III. 3.1  Magnetic  Field 

The  magnetic  field  results  seemed,  in  general,  quite 
reasonable  and  without  spurious  values.  Steep  gradients  occurred 
inside  the  thruster,  but  the  magnetic  field  strengths  increased 
monotonical ly  as  the  probe  moved  further  back  into  the  thruster. 
In  addition,  errors  which  would  be  caused  in  measurements  due  to 
effects  such  as  probe  ablation,  plasma  cooling,  and  displacement 
of  the  current  path  due  to  the  presence  of  the  probe  were 
examined  for  plasma  conditions  which  were  expected  in  the  near 
plume,  and  were  found  to  be  negligible.  However,  there  are  some 
factors  which  introduce  errors  into  these  measurements. 

Because  the  actual  signal  is  integrated,  the  error  resulting 
from  the  noise  should  be  reduced  because  a  positive  and  negative 
spike  will  come  close  to  cancelling  each  other  out.  However,  as 
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the  operational  amplifier  only  had  a  frequency  response  of  about 
10  MHz,  very  high  frequency,  large  amplitude  noise  spikes  may 
have  been  only  partially  eliminated.  In  other  words,  a  positive 
noise  spike  may  have  been  "seen”  by  the  integrator,  but  its 
opposing  negative  spike  which  occurred  only  a  few  nanoseconds 
later  may  have  been  missed  or  severely  attenuated.  Thus,  to  the 
integrator,  it  appears  as  if  there  is  a  net  positive  signal, 
rather  than  just  noise. 

At  some  locations,  namely  inside  the  thruster  near  the 
anode,  i.here  was  no  actual  flat  portion  of  the  trace.  Coupled 
with  rather  extreme  effects  of  electrical  noise  in  these  regions, 
there  is  a  fairly  significant  amount  of  error  in  these  regions. 

At  these  points,  three  firings  were  made  and  the  measurements 
were  averaged  in  an  attempt  to  get  a  mean  value  of  magnetic  field 
strength.  The  results  shown  in  figure  I I I. 6  and  table  III.l 
indicate  that  the  magnetic  field  at  the  anode  locations  inside 
the  thruster  still  seem  to  vary  smoothly  with  position,  but  the 
absolute  values  of  the  strength  may  be  somewhat  questionable. 

Looking  at  these  undesired  effects,  one  may  question  the 
calibration  technique,  as  the  probe  was  placed  very  close  to  the 
anode.  However,  the  voltage  traces  with  the  probe  near  the  back 
plate  were  flat  in  the  steady  state  portion  of  the  current  pulse, 
in  contrast  to  the  greatly  sloped  trace  in  the  steady  state 
portion  of  the  current  pulse  that  was  seen  during  firings  near 
the  anode  about  halfway  into  the  thruster.  When  the  probe  is  in 
these  locations,  it  appears  to  cause  a  significant  disturbance  in 
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the  plasma.  Due  to  a  lack  of  digitizers  in  the  oscilloscope,  a 
global  voltage  trace  of  the  thruster  while  making  magnetic  field 
measurements  was  not  acquired.  As  will  be  discussed  in  the  next 
section,  when  the  triple  Langmuir  probe  was  placed  in  these  same 
locations  the  global  voltage  showed  that  onset  was  occurring, 
although  no  indication  of  onset  was  apparent  when  the  probe  was 
in  other  locations. 

I I I. 3. 2  Langmuir  Probe 

Most  langmuir  probe  measurements  made  in  the  plume  and  near 
the  cathode  seemed  fairly  reasonable.  It  should  be  noted, 
however,  that  there  is  an  over  prediction  of  measured  electron 
temperatures  at  outer  radii  where  the  probe  was  not  aligned  with 
the  probe.  This  effect  will  be  discussed  later  in  this  section 
and  in  Appendix  C.  However,  all  of  the  firings  with  the  probe  in 
the  anode  region  inside  the  thruster  were  highly  questionable. 

Two  of  the  major  reasons  for  this  are  probe  contamination  and 
perturbation  of  the  plasma  by  the  probe. 

To  prevent  probe  contamination,  after  every  five  firings,  a 
glow  discharge  was  performed  in  order  to  clean  the  probes  by  ion 
bombardment.  When  the  probe  was  in  the  plume,  this  technique  was 
very  successful  in  keeping  the  contaminants  off  the  probe,  as 
evidenced  by  the  lack  of  spurious  results.  The  "good"  behavior 
of  the  signals  probably  resulted  because,  in  the  plume,  the  probe 
is  farther  away  from  the  sources  of  contaminants,  and  the  flux  of 
contaminant  particles  is  correspondingly  reduced.  However,  as 
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the  probe  was  moved  inside  the  thruster,  it  moved  closer  to  the 
sources  of  contamination,  namely  the  electrodes  and  the  back 
plate . 

Although  the  probe  voltage  and  current  traces  are 
reasonable,  the  absolute  magnitudes  of  the  probe  current  trace 
may  be  over  predicted  slightly  at  locations  near  the  back  of  the 
thruster.  The  overprediction  results  because  tungsten  from  the 
cathode  ablates,  and  may  attach  to  the  triple  probe.  This 
tungsten  coating  increases  the  effective  collection  area  of  the 
probe,  which  will  of  course,  increase  the  current  to  the  probe. 
The  coating  is  cumulative  for  a  constant  flux  of  contaminants. 

Near  the  cathode,  this  effect  may  have  been  alleviated 
somewhat  by  the  contamination  of  an  insulating  species,  boron 
nitride.  This  is  the  material  of  which  the  back  plate  is  made, 
and  the  flux  of  contaminant  particles  will  increase  as  the  back 
plate  is  approached.  Thus  the  increase  of  collecting  area 
resulting  from  the  deposition  of  tungsten  on  the  probe  electrodes 
may  have  been  reduced  by  the  deposition  of  insulation  material  on 
the  probe  as  well.  What  all  this  means,  is  that  by  the  time  the 
probe  approached  the  anode  in  the  test  sequence,  there  was  a 
great  deal  of  uncertainty  in  the  actual  collection  area  of  the 
probe.  As  the  probe  collection  area  is  a  parameter  used  in  the 
calculation  of  electron  number  density  and  electron  temperature 
from  the  voltage  and  current  readings,  a  great  deal  of 
uncertainty  was  introduced  into  the  calculated  values  of  those 
plasma  properties.  Although  glow  discharges  were  still 
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performed,  inside  the  thruster  the  mass  flux  of  contaminants  was 
high  enough  to  cause  concentrations  of  contaminants  on  the  probe 
which  were  unable  to  be  completely  removed  by  ion  bombardment. 

Another  factor  which  resulted  in  strange  results  inside  the 
thruster  was  the  perturbation  of  the  plasma  and  current  paths  by 
the  triple  Langmuir  probe.  As  mentioned  earlier,  when  the 
magnetic  probe  was  placed  in  certain  locations,  the  signal  traces 
did  not  behave  as  expected.  The  same  phenomenon  was  observed 
with  the  triple  Langmuir  probe,  especially  in  the  region  near  the 
anode  at  the  axial  location  corresponding  approximately  to  the 
cathode  tip.  This  point  shall  be  referred  to  as  the  "mid-anode", 
position.  Due  to  the  lack  of  oscilloscope  digitizers  in  the 
magnetic  probe  measurements,  monitoring  of  the  global  voltage  at 
the  same  time  as  magnetic  probe  voltage  was  not  possible. 

However,  with  the  addition  of  a  second  oscilloscope  for  the 
triple  Langmuir  probe  measurements,  global  values  were  monitored. 

When  the  triple  Langmuir  probe  was  placed  in  the  mid-anode 
location,  the  global  voltage  trace  indicated  that  onset  occurred. 
Two  other  locations  showed  indications  of  onset,  but  they  were 
not  as  "full-blown"  as  at  the  mid-anode  position.  In  addition 
the  other  two  points  of  significant  global  voltage  oscillations 
were  not  contiguous.  Looking  at  the  enclosed  current  contours 
shown  in  Figure  II. 1  earlier,  these  locations  have  many  current 
paths  pass  through  them.  When  the  exposed  wires  of  the  Langmuir 
probe  were  positioned  in  these  locations,  onset  was  induced. 
However,  when  the  insulated  portion  of  the  probe  was  in  the  same 
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location  no  significant  effect  on  the  global  voltage  was 
observed.  In  addition,  when  the  bare  wires  were  at  locations 
nearer  the  cathode  which  had  a  comparable  number  of  current 
paths,  onset  was  not  observed. 

At  the  locations  where  onset  was  induced,  the  voltage  traces 
were  as  much  as  an  order  of  magnitude  higher  than  at  other 
locations.  Probe  contamination  can  account  for  small  increases 
in  current  and/or  voltage  values,  but  not  of  this  scale. 
Obviously,  the  absolute  magnitude  of  any  number  density  or 
electron  temperature  calculations  in  these  regions  is  extremely 
suspect.  Further  investigation  of  the  perturbation  effects  of 
conducting  materials  placed  in  the  discharge  current  paths  is 
necessary  to  determine  the  reason  for  the  erroneous  probe 
measurements . 

However,  probe  contamination  will  not  affect  the  frequency 
of  the  fluctuations  to  the  probe,  merely  the  amplitude.  Thus, 
the  frequency  data  even  in  these  highly  questionable  regions  can 
probably  be  accepted,  albeit  with  some  skepticism.  In  fact, 
examining  the  frequency  data  at  these  highly  questionable  points, 
the  characteristic  frequencies  are  indicative  of  those  which  are 
expected  to  occur  at  onset. 

These  results  were  unexpected.  Previous  work  by  other 
experimenters  had  suggested  that  lower  hybrid  frequency 
oscillations  would  be  present  at  some  locations  [1/30].  However, 
calculation  of  Uie/Vti  based  on  mean  electron  temperature,  mean 
numoer  density,  and  magnetic  field  gradients,  the  values  range 
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from  .92  to  4.75.  According  to  Choueiri,  microinstabi 1 i t les  are 
supposed  to  occur  when  UJa/v..  >  1.5  [28].  As  stated  earlier, 
some  locations  will  yield  questionable  results  due  to 
contamination  and  perturbation  of  the  plasma,  and  there  may  be 
errors  of  as  much  as  a  factor  of  2  in  both  T?  and  n5  due  to  the 
non-alignment  of  the  probe  with  the  ion  flow  vector  [35] . 
Consequently,  although  the  above  ratios  appear  to  indicate  that 
microinstabilities  will  occur  at  some  locations,  all  of  these 
ratios  could  fall  below  1.5  when  experimental  error  is  taken  into 
account . 
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Table  III.l  -  Characteristic  Frequencies 


Point 


1.87 


3.35 


3.25 


5.15 
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Table  III. 2  Mean  Plasma  Properties 

Point 

WS8SSM 

n,.  ■  n'Jxio2-) 

V  0  S  A  Ju  *J  '  . 

1 

2.99 

6.76 

0.494 

7  .  36 

2 

2.72 

3.75 

0.633 

4.80 

3 

3.34 

1  .79 

1.02 

1.82 

4 

3.29 

0.672 

1 . 64 

0.736 

5 

0.258 

2.60 

0.316 

6 

3.39 

6.71 

0.528 

6.2 

7 

3.11 

4.56 

0.614 

4.83 

8 

2.79 

2.68 

0.759 

3.38 

9 

2.33 

1.45 

0.943 

2.4 

10 

2.11 

0.802 

1.21 

1.56 

11 

4.16 

6.84 

0.58 

4.81 

12 

3.03 

4.76 

0.594 

5.18 

13 

5.22 

1.83 

1.26 

1.02 

14 

4.81 

7.07 

0.613 

4.09 

15 

5.41 

14.0 

0.463 

6.67 

16 

5.29 

18.6 

0.397 

8.96 

17* 

7.25* 

1.74* 

1.52* 

0.626* 

18* 

8.37* 

1.4* 

1.82* 

0.419* 

19** 

59.9** 

3.77* 

2.96* 

0.073* 

20* 

8.83* 

2.18* 

1.5* 

0.594* 

21* 

9.76* 

4.84* 

1.06* 

1.11* 

22* 

12.2* 

4.14* 

1.27* 

0.71* 

23* 

9.79* 

3.82* 

1 . 19* 

0.882* 

*  questionable  data  point 


**  erroneous  data  point 


Figure  III.l  -  Gas  Pulse  with  Spark  trigger 
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Figure  II I. 2  -  Sample  current  and  voltage  traces  before  and  after 
onset 
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Figure  I I I. 5  -  Sample  magnetic  probe  voltage  trace 


Figure  III. 6  -  Spatial  Map  of  average  magnetic  field  with 

contours  for  6.8  kA  of  current. 
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Figure  III. 7  -  Sample  high  frequency  oscillation  signal  from  the 
triple  Langmuir  probe 
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Figure  III. 8  -  Sample  probe  voltage  and  current  traces  (absolute 
values ) 
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Figure  III. 9  -  Contour  map  of  electron  temperature  for  current  of 
6.8  kA. 


67 


IV.  Calculation  of  Plasma  Transport  Properties 
IV. 1  Bose  Method 

To  illustrate  the  effects  of  microinstabilities  on  the 
transport  properties  inside  the  thruster,  a  computer  code  which 
calculates  transport  properties  classically  according  to  kinetic 
theory  was  modified  to  account  for  microinstabilities.  The 
computer  code  used  was  developed  by  Tarit  Kumar  Bose[42,43],  and 
it  assumes  a  two-temperature,  multiple  species,  multiply  ionized 
plasma.  As  a  two  temperature  plasma  is  a  thermal  non-equilibrium 
tase  (which  is  one  reason  these  instabilities  can  occur)  it  is 
assumed  that  the  particle  temperatures  are  perturbed  from 
equilibrium.  This  assumption  is  made  because  transport 
properties  are  derived  from  a  first  order  perturbation  of  the 
equilibrium  state. 

Thus,  a  quasi-equilibrium  state  must  be  determined,  assuming 
the  state  can  be  described  in  terms  of  the  total  pressure,  p,  the 
heavy  particle  translational  temperature,  Ta,  and  the 
translational  temperature  of  the  electrons,  Tg.  Also,  it  is 
assumed  that  the  species  excitation  temperature  is  equal  to  Te . 
Mole  fractions  of  the  various  ionized  species  are  dttermined 
through  the  use  of  generalized  Saha  equations  according  to  the 
Monti-Napol itano  and  Veis  model [44 , 45] ,  which  is  based  on  the 
principle  that  the  total  species  partial  pressure  does  not  change 
due  to  chemical  reactions.  However,  it  must  be  noted  that  the 
use  of  the  generalized  Saha  equation  assumes  a  quasi -equilibrium 
state.  In  terms  of  an  MPD  plasma,  this  assumes  that  the  time  it 
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takes  for  a  particle  to  ionize  is  much  less  than  the  ion 
residence  time  in  the  thruster.  Unfortunately,  it  is  not  known 
if  this  condition  is  met  in  the  MPD  thruster  due  to  the  high  ion 
drift  velocities,  on  the  order  of  1-5  km/sec  [46].  For  the 
purpose  of  this  analysis,  it  will  be  assumed  that  the  particles 
have  time  to  ionize  before  they  leave  the  thruster.  It  can  be 
demonstrated  that  the  transport  properties  for  collisions  in 
which  one  of  the  colliding  partners  is  an  electron  and  collisions 
between  only  heavy  particles  can  be  evaluated  separately [ 42 ] . 
Viscosity,  translational  heat  conductivity,  and  electrical 
conductivity  are  those  properties  which  can  put  this  fact  to  use. 
Other  transport  properties  must  use  an  ultra-simplified  method 
[47];  however,  since  the  only  transport  property  examined  in  this 
report  is  electrical  conductivity,  this  ultra-simplified  theory 
will  not  be  discussed  further.  The  reason  for  this  emphasis  on 
electrical  conductivity  is  that  the  current  which  supplies  the 
thrust  for  an  MPD  device  is  dissipated  directly  by  electrical 
resistivity. 

To  calculate  the  electrical  conductivity,  a  fourth  order 
Sonine  polynomial  expansion  was  used.  Chapman  and  Cowling  have 
shown  that  transport  properties  can  be  expressed  exactly  as 
ratios  of  infinite  determinants  [47].  However,  these  ratios 
converge  fairly  rapidly  as  more  rows  or  columns  are  adCed  to  the 
determinant.  Thus,  one  can  get  a  good  approximation  of  the 
actual  transport  properties  using  a  ratio  of  finite  determinants. 
In  this  paper  and  in  Bose's  code,  a  fourth  order  expansion  was 
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deemed  sufficient. 

The  equation  used  to  calculate  electrical  conductivity  is 


o  = 


1  ■  n 


(IV. 1) 


[r> 


3 xs(  2iticre\l/i  jg::i 

me  j  jgj 


(IV. 2) 


where  1  q j  is  the  determinant  of  the  Sonine  polynomial,  !q<>!  is 

the  first  minor  of  the  Sonine  polynomial,  xe  is  the  electron  mole 

fraction,  and  n  is  total  number  density.  [ Dee 3 ^  is  defined  as  the 

fourth  order  electron-electron  diffusion  coefficient. 

In  order  to  calculate  the  determinants,  it  is  first 

necessary  to  calculate  the  Sonine  polynomials,  which,  in  turn, 

depend  on  the  collision  cross  sections.  Bose's  code  uses 

different  collision  potentials  for  different  type  of  collisions. 

The  only  gas  which  is  examined  in  this  study  is  argon,  as  that  is 

the  gas  that  was  used  in  the  experimental  work  of  chapter  III. 

For  argon  the  following  potentials  are  used  to  determine 

collision  cross  sections. 

Neutral-neutral:  exponential  repulsive  potential 
Electron-Neutral  :  experimental  gas-kinetic  cross  sections 
Ion-Neutral:  charge  transfer  cross  sections 
Ion-Electron:  attractive  coulomb  potential 
Ion-Ion:  repulsive  Coulomb  potential 
Electron-Electron:  repulsive  Coulomb  potential 

For  the  remainder  of  this  paper,  the  use  of  the  term  "Bose 

method"  will  refer  to  Bose's  method  of  calculation  of  transport 

properties  in  a  two  temperature  plasma  according  to  kinetic 

theory. 
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IV. 2  Caldo-Choueiri  Calculation  of  Transport  properties 

Giuliano  Caldo,  Edgar  Choueiri,  Arnold  Kelly  and  Robert  Jahn 
at  Princeton  University  have  developed  a  numerical  code  for  a  two 
fluid  case  which  takes  anomalous  transport  into  account  self- 
consistently  for  the  solution  of  a  two-dimensional  MPD  thruster 
flow  field  [28,29].  It  employs  a  modified  finite  difference 
MacCormack  code  for  the  solution  of  the  two-fluid  conservation 
equations,  while  simultaneously  solving  Maxwell's  equations. 

Several  assumptions  are  made.  First,  the  plasma  is  assumed 
fully  ionized,  so  no  neutrals  are  present.  Also,  it  is  assumed 
that  each  fluid  obeys  the  ideal  gas  law  separately.  Flow  is  als.o 
assumed  to  be  inviscid.  The  type  of  microinstability  that  is 
expected  is  the  generalized  lower  hybrid  drift  instability,  which 
operates  at  the  lower  hybrid  frequency.  Because  ion  trapping  has 
been  postulated  to  be  the  mechanism  for  non-linear  saturation  of 
this  instability,  the  calculation  of  anomalous  transport  assumes 
that  ion  trapping  is  the  means  of  saturation.  Lastly,  it  is 
assumed  that  the  flow  can  be  modeled  by  a  two  fluid  model  as 
opposed  to  kinetic  theory. 

In  addition,  calculation  of  electrical  conductivity  uses  the 
mean-free  path  formula  to  find  an  effective  value  of  conductivity 
which  accounts  for  anomalous  transport,  oeff*. 


e2n< 


(IV. 3) 


where  (va)p  =  the  anomalous  momentum  exchange  collision 


frequency.  Note  that  only  electron-ion  collisions  are  accounted 


for  in  determining  conductivity  if  this  formula  is  used: 
electron-electron  and  electron-neutral  collisions  are  neglected. 
For  the  remainder  of  the  paper,  this  method  of  calculation  of 
both  the  classical  and  effective  electrical  conductivity  will  be 
referred  to  as  the  "mean  free  path  method" .  The  anomalous 
collision  frequency  is  calculated  using  a  curve  fit  based  on 
calculations  by  Choueiri [4, 28] : 

(v/)^  3  v#.  [0.192  ♦  0.033Q.  -  0 . 2 1 2  Q  „2  -  3.27xlCT5Q93 

T  fry  41 

+  (0.00123  -  0.0158Q,  -  0.00789£Jo2)] 

IV. 3  Modification  of  Code 

The  key  differences  between  the  Bose  method  and  the  mean 

*  .  I 

free  path  method  used  by  Caldo  and  Choueiri  are  that  the  Bose 

method  uses  kinetic  theory  as  opposed  to  fluid  theory,'  a  4th 

order  Sonine  Polynomial  expansion  is  used  to  calculate  electrical 

conductivity  rather  than  the  mean-free  path  formula,  and 

electron-electron  and  electron-neutral  collisions  are  taken  into 

* 

account.  However,  Bose's  computer  code  only  calculates  point 

transport  properties.  If  Bose’s  computer  code  is  modified  for 

.# 

microinstabilities,  it  can  be  used  as  a  subroutine  in  a  numerical 

i 

model  such  as  that  developed  by  Caldo  and  Choueiri  in  order  to 
get  a  spatial  mapping  of  transport  properties.  From  this  point 
on,  the  term  "effective"  will  refer  to  properties  which  account 
for  both  classical  and  anomalous  effects  resulting  from  the 
presence  of  microinstabilities,  while  the  term  "anomalous"  will 
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refer  to  effects  stemming  only  from  microinstabilities. 

To  account  for  microinstabilities,  the  collision  cross- 
sections  for  electron-ion  collisions  were  modified.  This 
modification  involves  the  use  of  the  curve-fit  in  equation  (IV. 4) 
expressing  the  ratio  of  anomalous  collision  frequency  to 
classical  electron-ion  collision  frequency  as  a  function  of 
electron  hall  parameter.  As  Caldo  and  Choueiri  used  only  the 
electron-ion  collision  frequency  in  determining  the  electron  Hall 
parameter,  this  author  has  done  the  same  to  maintain 
compatibility  in  the  results.  In  any  event,  for  the  conditions 
measured  in  the  MPD  thruster,  the  electron-heavy  particle 
collision  frequency  is  identical  to  the  electron-ion  collision 
frequency  because  no  neutrals  are  present. 

The  lower  hybrid  wave  typically  has  a  phase  velocity  on  the 
order  of  the  ion  thermal  velocity  [4],  Because  electrical 
resistivity  is  effectively  a  measure  of  the  "drag"  on  electrons, 
usually  resulting  from  collisions,  the  interaction  of  the  wave 
with  electrons  increases  resistivity  by  increasing  the  "drag". 

In  calculating  the  effective  collision  cross  section  from  the 
collision  frequency,  the  relative  velocity  between  the  wave  and 
electrons  must  be  used.  However,  as  the  electron  thermal 
velocity  is  much  larger  than  the  wave  velocity  (ion  thermal 
velocity),  the  relative  velocity  between  the  wave  and  electron 
was  approximated  by  the  electron  thermal  velocity.  Lastly,  since 
the  electrons  are  the  particles  which  are  related  to  electrical 
conductivity,  the  number  of  electrons  has  a  direct  impact  on  the 
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collision  frequency.  Thus,  using  the  formula, 

v\w  =  Q?as  (IV.  5) 

where  Qr;r  s  anomalous  momentum  transfer  collision  cross  section, 
the  anomalous  collision  cross  section  can  be  solved  for  based  on 
the  resulting  value  of  anomalous  collision  frequency  from  the 
curve  fit.  This  collision  cross  section  is  then  superimposed  on 
the  electron-ion  collision  cross  section  to  create  an  effective 
collision  cross  section.  In  order  to  calculate  the  transport 
properties  by  the  Bose  method  for  a  two  temperature  plasma  while 
accounting  for  microinstabilities,  this  cross  section  is  used. 

In  order  to  calculate  the  transport  coefficients  according 
to  the  Bose  method,  collision  cross  section  data  for  collisions 
which  have  different  speed  and  angular  dependencies  than  that 
based  on  a  rigid-sphere  model  (the  maximum  cross  section)  are 
required  for  the  Sonine  polynomial  expansion.  These  additional 
cross  sections  are  available  for  use  in  the  classical  calculation 
of  transport,  but  this  body  of  data  is  non-existent  for  anomalous 
transport.  A  key  assumption  that  was  made  in  the  modification  of 
the  cross  sections  to  account  for  microinstabilities  is  that  the 
additional  cross  sections  used  to  determine  effective  electrical 
conductivity  have  a  similar  distribution  to  the  distribution  of 
the  classical  cross  sections.  As  the  Caldo-Choueiri  curve  fits 
produce  a  maximum  cross  section,  the  deviations  from  this  cross 
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section  were  based  upon  the  following  scaling  relation. 


In  this  equation,  Q -  - =  the  modified  "effective"  collisj-n 
cross  section  accounting  for  deviations  from  the  rigid  sphere 
model,  Q  =  the  experimentally  determined  classical  cross 

sections  which  deviate  from  the  rigid  sphere  model,  Q121.:  =  the 
maximum  classical  cross  section  (the  one  corresponding  to  a  rigid 
sphere  model},  and  Qma.:.-  =  the  cross  section  calculated  directly 
from  equations  (IV. 4)  and  (IV. 5).  It  must  be  noted  that  until 
reliable  values  of  cross  sections  used  in  the  Sonine  polynomial 
expansion  for  anomalous  cross  sections  can  be  determined 
experimentally,  the  results  shown  in  this  paper  have  a  degree  of 
uncertainty  which  cannot  be  quantified.  However,  this  is  the 
best  approximation  that  can  be  made  at  the  current  time-. 

In  addition,  it  should  be  noted  that  this  modification  of 
collision  cross-section  is  only  good  for  wave-electron 
interactions.  If  the  same  method  is  attempted  for  modification 
of  wave-ion  or  wave-neutral  interactions,  the  effective  collision 
cross  sections  can  become  extremely  large  because  ions  and 
neutrals  move  at  almost  the  same  speed  as  the  lower  hybrid  wave, 
resulting  in  a  relative  velocity  of  zero.  To  calculate  the 
effective  cross  section  from  the  Caldo-Choueiri  curve  fits,  the 
collision  frequency  must  be  divided  by  the  relative  velocity. 

This  would  introduce  an  almost  infinite  collision  cross  section. 

This  makes  sense  for  ions.  The  Caldo-Choueiri  curve  fits 


assume  that  ion  trapping  is  the  means  of  damping  the  collisions. 
As  a  result,  the  ions  do  experience  an  'nfinite  collision  cross 
section  with  the  wave  because  they  cannot  escape  from  it;  they 
are  "trapped”.  However,  neutrals  should  not  experience  this  same 
effect  because  there  is  no  net  charge  on  the  neutrals,  so  a 
potential  well  should  not  be  able  to  trap  them  as  easily. 

IV. 4  Comparison  of  Results 

IV. 4.1  Comparison  between  mean  free  path  and  Bose  methods 

To  get  an  accurate  comparison  of  the  mean  free  path  and 
Bose  methods  of  calculating  electrical  conductivity,  the  number  . 
densities  of  each  species  must  be  kept  the  same,  regardless  of 
method  of  calculation.  This  leads  to  the  use  of  only  one  method 
of  calculating  the  species  composition  of  the  gas.  Although  the 
mean- free  path  method  only  requires  the  electron  number  density, 
heavy  species  number  densities  are  required  for  the  Bose  method. 
Therefore,  the  method  of  calculating  species  composition  of  the 
gas  which  was  used  in  Bose's  computer  code  was  also  used  to 
determine  species  number  densities  for  the  mean-free  path  method 
as  well  as  the  Bose  method. 

The  electrical  conductivities  calculated  by  both  the  Bose 
and  mean-free  path  methods  were  calculated  for  a  variety  of 
conditions.  The  pressure  was  held  approximately  constant,  while 
the  electron  temperature  and  ratio  of  electron-to  heavy  particle 
temperatures  were  varied.  The  range  of  electron  temperatures 
spanned  from  5000  K  to  35000  K,  and  the  temperature  ratios 
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spanned  from  .5  to  1.8.  A  plot  of  electrical  conductivity  vs. 
electron  temperature  is  shown  in  figure  IV. 1  for  a  one- 
temperature  plasma,  Hall  parameter  of  1,  and  pressure  of  10  Pa. 

Examining  this  figure,  one  may  notice  the  discontinuities  in 
the  electrical  conductivities  at  a  few  temperatures,  regardless 
of  the  method  of  calculation  used.  As  Bose's  code  takes  into 
account  multiply  ionized  species,  there  are  some  discontinuities 
in  electrical  conductivities  at  temperatures  where  the  primary 
heavy  species  changes  ionization  levels.  This  is  the  reason  for 
local  peaks  in  the  distribution  of  electrical  conductivity  vs. 
temperature. 

Looking  at  the  case  involving  the  mean  free  path  method,  one 
can  see  that  the  difference  between  classical  and  electrical 
conductivity  appears  to  increase  with  higher  temperature. 

However,  if  the  ratio  of  effective  to  classical  conductivity  is 
examined,  as  shown  in  figure  IV. 2,  the  ratio  is  constant 
throughout  the  entire  temperature  range,  showing  the  effective 
conductivity  to  be  approximately  30  per  cent  lower  than  the 
classical  conductivity.  The  reason  that  this  ratio  does  not  vary 
is  that  there  is  no  explicit  absolute  temperature  dependence  in 
the  calculation  of  the  anomalous  collision  frequency.  Although 
electron  temperature  has  an  effect  on  the  electron-ion  collision 
frequency  used  in  determination  of  electrical  conductivity,  this 
factor  is  present  in  both  effective  and  classical  conductivity 
calculations,  and  cancels  out  when  a  ratio  is  taken. 

In  contrast,  absolute  temperature  does  play  a  significant 
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role  in  the  ratio  of  effective  to  classical  conductivity  using 
the  Bose  method  of  calculating  transport  properties.  Comparing 
these  results  to  the  results  from  the  mean-free  path  method  in 
figure  IV. 1,  one  can  see  that  the  electrical  conductivities 
follow  the  same  general  trend  regardless  of  the  method  used. 

That  is,  both  methods  show  an  increase  in  electrical  conductivity 
with  temperature  and  both  methods  show  local  peaks  resulting  from 
transition  in  ionization  state  of  the  primary  heavy  species. 
However,  the  absolute  values  of  classical  and  effective 
electrical  conductivities  calculated  by  the  Bose  method  are  both 
higher  than  for  both  the  classical  and  effective  values 
calculated  by  the  mean-free  path  method. 

This  discrepancy  is  due  to  the  methods  being  used  to 
calculate  the  conductivity.  Mitchener  and  Kruger  state  that  the 
mean-free  path  formula  gives  a  good  order  of  magnitude  estimate 
of  electrical  conductivity,  but  can  be  off  by  as  much  as  a  factor 
of  two  [32].  They  also  state  that  better  values  can  be  obtained 
by  using  Sonine  polynomial  approximations  [32].  The  higher  order 
the  Sonine  polynomial  approximation,  the  more  accurate  the 
result.  However,  the  Sonine  polynomials  converge  fairly  rapidly, 
so  the  fourth  order  approximation  used  by  Bose  in  the  calculation 
of  electrical  conductivity  is  fairly  close  to  the  exact  result. 
Examining  Figure  IV. 1  again,  it  can  be  seen  that  the 
discrepancies  between  the  classical  values  of  electrical 
conductivity  differ  by  a  factor  no  greater  than  2,  which  is 
within  the  uncertainty  of  the  mean-free  path  formula. 
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However,  the  effective  conductivity  based  on  the  Bose  method 
behaves  much  differently  than  that  calculated  by  the  mean  free 
path  method.  In  the  mean  free  path  method,  the  effective 
conductivity  diverges  from  the  classical  conductivity  slightly 
with  increasing  temperature,  but  as  explained  earlier,  this  is 
due  to  the  constant  ratio  between  effective  and  classical 
conductivities.  In  contrast,  for  the  Bose  method,  the  curves 
actually  converge.  The  effect  on  the  ratio  of  effective  to 
classical  electrical  conductivity  can  be  seen  figure  IV. 2. 

At  first  glance,  these  results  appears  strange.  The  same 
curve  fit  for  anomalous  collision  frequency  was  used  as  the  basis 
for  calculations  of  effective  conductivity  in  both  the  mean  free 
path  and  Bose  methods.  However,  the  curve  fit  had  no  absolute 
temperature  dependence.  Thus,  one  may  wonder  why  the  results 
using  the  Bose  method  show  a  clear  temperature  dependence. 

The  answer  lies  in  the  interactions  which  were  accounted  for 
in  the  Bose  method,  which  were  not  accounted  for  by  the  mean  free 
path  method.  The  mean  free  path  method  took  only  electron-ion 
interactions  into  account  in  the  determination  of  electrical 
conductivity.  One  reason  Caldo  and  Choueiri  limited  the  scope  of 
their  analysis  using  the  mean  free  path  method  is  because  of  the 
assumption  that  the  anomalous  effects  are  a  result  of  ion- 
electron  interactions  only,  provided  ion  trapping  is  the  dominant 
mode  of  saturation  of  microinstabilities  [28]. 

The  Bose  method  takes  into  account  several  other 
interactions,  though:  electron-electron,  electron-neutral , ion- 
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neutral,  ion-ion,  and  neutral-neutral  in  addition  to  electron-ion 
interactions.  It  appears  that  these  effects  "dilute"  the  effects 
of  the  anomalous  collision  frequency.  One  may  think  of  the 
effect  of  accounting  for  additional  interactions  as  one  may  think 
of  adding  dye  to  a  glass  of  water. 

Suppose  a  glass  contained  some  water,  and  some  red  dye  was 
added  to  it  so  that  whatever  the  initial  amount  of  water,  the 
color  of  the  water  would  be  the  same  immediately  after  the  dye 
was  added,  say  a  bright  red.  However,  as  more  and  more  water  is 
added,  the  dye  becomes  harder  to  see.  Using  this  analogy,  the 
dye  represents  the  effects  of  anomalous  collision  frequency,  and 
the  additional  water  represents  additional  collision  types  which 
are  accounted  for.  Obviously,  if  these  additional  collision 
effects  are  very  small  compared  to  the  original  collisions  which 
are  accounted  for,  the  effect  of  accounting  for  these  collisions 
will  not  be  very  great,  just  as  adding  only  a  little  water  to  the 
dyed  water  may  not  noticeably  change  the  color.  However,  as 
these  interactions  become  more  significant,  the  anomalous  effects 
which  are  associated  with  only  one  type  of  interaction  (electron- 
ion)  are  less  noticeable.  Using  the  above  analogy,  the  color  of 
the  water  becomes  noticeably  dimmer.  Eventually,  enough  water 
may  be  added  such  that  the  dye  cannot  be  seen  at  all. 

This  is  what  happens  at  very  high  and  at  very  low 
temperatures  in  the  MPD  thruster  plasma.  At  lower  temperatures 
(about  5000  K),  the  plasma  is  only  very  weakly  ionized;  in  other 
words,  neutrals  are  the  primary  species,  so  electron-neutral. 
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neutral-neutral,  and  ion-neutral  collisions  will  dominate  the 
plasma.  Electron-ion  collisions  are  very  insignificant  compared 
to  these  interactions  because  both  ion  and  electron  mole 
fractions  are  extremely  low.  As  a  result,  the  ratio  of  effective 
to  classical  conductivity  should  be  one.  Once  again  examining 
figure  IV. 2,  this  can  be  seen  to  be  the  case. 

However,  as  the  plasma  becomes  fully  ionized,  the  neutral 
mole  fractions  approach  zero,  and  all  collisions  involving 
neutrals  become  negligible.  In  contrast,  the  mole  fractions  of 
both  ions  and  electrons  increase  significantly.  Around  1  eV,  or 
11600  K,  the  mole  fractions  of  ions  and  electrons  are  roughly 
equal.  Thus,  at  this  point,  electron-electron,  electron-ion,  and 
ion-ion  interactions  are  significant,  although  electron-ion 
effects  should  be  dominant  due  to  the  much  larger  cross  sectional 
area  of  an  ion  as  opposed  to  an  electron.  In  addition,  as 
electrical  conductivity  involves  electrons,  the  ion-ion 
interactions  are  not  very  important  in  comparison  to  the  other 
two  interactions.  Thus,  the  effective  conductivity  should  be 
much  lower  at  this  temperature  range.  Once  again,  figure  IV. 2 
shows  that  the  ratio  of  conductivities  drops  to  a  much  lower 
value,  and  it  almost  reaches  the  level  of  the  ratio  predicted  by 
Caldo  and  Choueiri.  The  difference  between  the  two  ratios  at 
this  fully  ionized  stage  is  due  to  the  electron-electron 
interactions  accounted  for  in  the  Bose  method  which  are  not 
accounted  for  in  the  mean  free  path  method. 

However,  as  the  plasma  becomes  even  more  highly  ionized,  the 
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ratio  of  conductivities  according  to  the  Bose  method  starts  to 
increase,  reaching  about  0.92  at  about  3  eV  (34,800  K).  As 
neutrals  are  even  less  likely  at  these  temperatures,  some  type  of 
interaction  must  be  becoming  more  important  than  electron- ion 
interactions.  These  dominating  interactions  are  electron- 
electron  collisions.  As  the  plasma  becomes  more  highly  ionized, 
more  electrons  are  present  for  each  ion,  so  the  electron  mole 
fraction  increases,  while  the  heavy  particle  mole  fractions 
decrease.  Therefore,  if  the  temperature  is  raised  high  enough, 
electron-ion  interactions  will  be  so  small  compared  to  electron- 
electron  interactions  that  anomalous  effects  will  no  longer  be 
noticeable . 

The  effects  of  Hall  parameter  on  the  effective  and  classical 
conductivities  are  shown  in  Figure  IV. 3,  for  a  single  temperature 
plasma  at  10  Pa,  and  electron  temperature  of  31500  K.  Unlike 
the  classical  case,  the  Hall  parameter  has  a  noticeable  effect  on 
the  effective  electrical  conductivity  for  both  methods  of 
calculation.  As  the  Hall  parameter  increases,  the  effective 
conductivity  drops.  This  result  makes  sense  because  equation 
(IV. 4)  has  a  strong  dependence  on  Hall  parameter. 

The  ratios  of  effective  to  classical  electrical  conductivity 
fo1"  the  above  conditions  are  plotted  vs.  Hall  parameter  in  Figure 
IV. 4.  From  this  plot  it  is  clear  that  although  electrical 
conductivity  drops  significantly  with  Hall  parameter  regardless 
of  the  method  used,  this  drop  is  much  less  drastic  if  the  Bose 
method  is  used.  This  in  turn  means  that  electron-electron 
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collisions  appear  to  have  a  significant  mitigating  effect  on 
effective  conductivity  at  low  to  intermediate  values  of  Hall 
parameter,  although  this  mitigating  effect  diminishes  at  very 
high  values  of  the  Hall  parameter.  The  mean  free  path  case 
matches  the  results  shown  by  Caldo  and  Choueiri  [2,4,28],  albeit 
in  slightly  different  form. 

The  effects  of  temperature  ratio  will  be  discussed  in  the 
next  section. 


IV. 4. 2  Illustrative  Example 

Using  data  acquired  at  different  locations  inside  the 
thruster  through  the  experiments  described  in  chapter  III,  the 
Bose  code  was  applied  to  determine  electrical  conductivity.  In 
the  experimental  work,  microinstabilities  were  not  seen  at  any 
locations.  However,  assuming  microinstabilities  exist  at  these 
locations,  the  following  illustration  will  show  the  differences 
in  conductivities  which  occur  due  to  anomalous  effects  at  two 
different  locations  in  the  thruster  and  its  plume. 

The  two  locations  selected  had  the  following  properties: 
Point  #1: 

Position:  7.7  cm  from  the  back  plate,  2.5  cm  from  thruster 
axis . 

Pressure  =  22  Pa. 

Electron  Temperature  =  31500  K. 

Magnetic  Field  =  17  Gauss. 

Hall  Parameter  =  .014 

Point  #2: 

Position:  4.7  cm  from  the  back  plate,  1.5  cm  from  the 
thruster  axis. 

Pressure:  57.5  Pa 

Electron  Temperature  =  48200  K 

Magnetic  Field  =  224  Gauss. 
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Hall  Parameter  =  .128 

The  pressures  were  chosen  so  as  to  match  the  measured  electron 
number  density  assuming  a  single  temperature  plasma.  In  terms  of 
the  positions  in  relation  to  the  thruster  electrodes,  point  #1 
was  about  halfway  between  the  electrodes  in  the  radial  direction, 
and  about  3  cm  beyond  the  exit  plane  of  the  thruster.  Point  #2 
was  located  fairly  close  to  the  tip  of  the  cathode.  These  points 
were  chosen  because  they  were  close  enough  to  the  axis  where 
error  in  the  alignment  of  the  probes  would  not  be  as  great,  they 
were  positions  which  were  believed  to  have  produced  good  data, 
and  point  #2  corresponded  to  the  location  where  there  is  a  very  - 
slight  possibility  that  some  lower  hybrid  oscillations  may  have 
occurred. 

Figure  IV. 5  shows  plots  of  the  classical  conductivities  vs 
the  temperature  ratio,  Tg/T„.  As  can  be  seen  for  the  classically 
calculated  values,  deviations  of  temperature  ratio  from  unity 
result  in  higher  electrical  conductivities.  In  addition,  the 
conductivities  calculated  by  the  different  methods  differ  by  a 
factor  of  approximately  two.  Figure  IV. 6  shows  the  effective 
conductivities  vs.  the  temperature  ratio.  Figures  IV. 5  and  IV. 6 
are  virtually  identical  with  the  exception  that  all  values  of 
conductivity  for  both  thruster. locations  are  lower  for  the 
effective  cases. 

To  see  the  effect  that  temperature  ratio  has  on  effective 
conductivity,  one  may  look  at  Figure  IV. 7,  a  plot  of  the  ratio  of 
effective  to  classical  conductivity  vs.  temperature  ratio  for 
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both  thruster  locations.  As  can  be  seen,  neither  method  shows 
much  variation  at  these  locations  for  either  case.  The  only 
noticeable  effect  is  that  the  mean  free  path  method  predicts  that 
conductivity  will  be  approximately  12  per  cent  lower  than  for  the 
Bose  method  at  these  points.  Higher  or  lower  temperature  ratios 
will,  of  course  produce  greater  effects  on  these  ratios  of 
conductivity;  however,  the  temperature  ratios  shown  are 
indicative  of  those  occurring  in  an  MPD  thruster.  Although 
temperature  ratio  does  have  a  small  effect  on  the  effective 
conductivity,  examination  of  figure  IV. 5  seems  to  indicate  that 
the  classical  effects  of  temperature  ratio  are  much  more 
significant  than  any  effects  temperature  ratio  has  on  anomalous 
transport . 

One  may  note  that  both  of  these  locations  are  characterized 
by  fairly  low  Hall  parameters  and  rather  high  electron 
temperatures  and  number  densities.  As  anomalous  transport  seems 
most  strongly  related  to  Hall  parameter,  it  appears  that  in  high 
power  devices  such  as  the  one  used  in  this  work,  anomalous 
transport  may  not  have  much  of  an  effect.  However,  this  analysis 
seems  to  show  that  the  observable  effects  of  anomalous  transport 
at  high  power  levels  are  less  if  the  Bose  method  is  used  than  if 
the  mean  free  path  method  is  u$ed.  In  a  spatial  sense,  however, 
effective  conductivity  calculated  by  the  Bose  method  seems  to 
exhibit  more  of  a  spatial  dependence  th an  that  calculated  by  the 
mean  free  path  method.  This  dependence  is  evidenced  by  the 
larger  separation  between  curves  representing  different  locations 
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in  Figure  IV. 7  for  the  Bose  case  than  for  the  mean  free  path 
case. 
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Figure  IV. 1  -  Electrical  Conductivity  vs.  Electron 
Temperature  for  a  one  temperature  plasma  at  a  pressure  of  10 
Pa.  Results  for  both  effective  and  classical  results  based 
on  both  the  Bose  and  mean  free  path  methods  are  shown. 
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Figure  IV. 2  -  Ratio  of  ef fective/classicai  electrical 
conductivity  vs.  electron  temperature  for  both  the  Bose  and  mean 
free  path  methods. 


Figure  IV. 3  -  Electrical  conductivity  vs.  Hall  parameter  for 
a  one  temperature  plasma  at  a  pressure  of  10  Pa.  Results  for 
both  the  Bose  and  mean- free  path  methods  are  shown. 
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Figure  IV. 4  -  Ratio  of  effective  to  classical  conductivity  vs. 
Hall  parameter  for  pressure  =  10  Pa,  single  temperature  plasma  at 
31500  K.  Both  the  Bose  method  and  mean  free  path  method  are 
shown . 


Figure  V.5  -  Classical  electrical  conductivities  vs.  ratio  of 
electron  to  heavy  particle  temperature  by  both  the  Bose 
method  and  mean  free  path  method.  Point  1  and  Point  2  are 
two  different  locations  in  the  thruster. 
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Figure  IV. 6  -  Effective  electrical  conductivity  vs.  ratio  of 
electron  to  heavy  particle  temperature  by  both  Bose  and  mean  free 
path  methods.  Point  1  and  Point  2  are  two  different  locations  in 
the  thruster. 
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V.  Discussion  and  Recommendations 
V . 1  Tentative  Conclusions 

Some  tentative  conclusions  may  be  drawn  from  both  the 
experimental  and  theoretical  results  of  this  work. 

1)  Microinstabilities  do  not  seem  to  occur  anywhere  in  the 
MW  level  thruster  at  operating  conditions  slightly  below  onset; 
however,  much  lower  frequency  instabilities  do  appear  to  exist, 
and  are  probably  related  to  onset.  Investigations  at  different 
power  levels  and  repeatability  studies  are  required  to  determine 
if  the  lack  of  microinstabilities  is  a  function  of  high  power 
devices  in  general,  a  function  of  operating  at  conditions  closer 
to  onset  than  attempted  in  previous  work  [1,30],  or  just  an 
erroneous  test  run. 

2)  Unlike  the  situation  where  a  "stinger"  is  introduced 
into  an  MPD  flow  field,  which  delays  onset[40],  it  appears  that 
the  insertion  of  electrical  conductors  into  the  flow  field  near 
the  anode  acts  to  induce  onset. 

3)  Effective  conductivity  appears  to  have  an  indirect 
temperature  dependence  which  is  greater  than  its  direct 
dependence  on  temperature  ratio  and  Hall  parameter.  This 
temperature  dependence  results  from  a  variance  in  species  number 
densities  with  temperature  and  does  not  seem  to  arise  to  unless 
electron-electron  and  electron-neutral  interactions  are  taken 
into  account. 

4)  Finally,  for  conditions  which  appear  to  be  present  inside 
a  MW  level  MPD  thruster  and  in  its  near  plume,  the  effect  of 


anomalous  transport  induced  by  lower  hybrid  micro i nstabi 1 1 t xes 
seems  to  be  less  significant  than  had  been  anticipated.  Assuming 
that  ionization  occurs  according  to  classical  methods,  if 
microinstabilities  are  present  there  is  only  a  small  reduction  in 
conductivity  for  conditions  present  in  a  MW  level  device,  on  the 
order  of  a  5  per  cent.  It  must  be  kept  in  mind,  however,  that  if 
ionization  does  not  occur  according  to  classical  methods,  the 
effect  of  anomalous  transport  may  be  as  great  as  anticipated. 

V.2  Recommendations 

In  the  realm  of  experimental  work,  many  improvements  can  be 
made.  First  and  foremost,  the  effect  of  power  level  on 
microinstabilities  needs  to  be  investigated,  as  well  as 
verification  of  the  results  presented  in  this  paper.  Secondly,  a 
means  of  aligning  a  Langmuir  probe  with  the  flow  will  greatly 
increase  the  accuracy  of  measurements  of  mean  electron 
temperature  and  electron  number  density.  Thirdly,  a  study  of 
the  perturbation  effects  of  an  electrical  conductor  placed  in  the 
current  paths  of  the  thruster  would  help  make  evaluation  of  any 
interior  thruster  data  more  accurate.  Fourth,  the  use  of  a  large 
number  of  replacement  probes  when  taking  measurements  inside  the 
thruster  would  greatly  reduce  the  problems  involved  with  probe 
contamination.  Finally,  if  microinstabilities  are  detected  at 
different  power  levels  in  this  device,  phase  and  wavelength  data 
would  complement  the  frequency  data  to  perhaps  identify 
specifically  which  microinstability  was  operating. 
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On  the  theoretical  side  of  this  paper,  although  the  Bose 
method  shows  clearly  tne  effect  of  accounting  for  different  types 
of  collisions  throughout  the  entire  range  of  temperatures  likely 
to  be  seer  in  an  MPD  thruster,  a  small  modification  of  the  mean 
free  path  method  used  by  Caldo  and  Choueiri  method  could  increase 
the  applicability  of  their  results  while  increasing  the  work 
infinitesimally.  This  modification  would  be  to  account  for 
electron-electron  collisions  in  the  mean-free  path  formula.  Rule 
of  thumb  formulas  for  electron-electron  collision  frequency  are 
available  from  many  sources,  such  as  Mitchener  and  Kruger  [32], 
and  Chen[6].  The  electron-ion  and  electron-electron  collision 
frequencies  may  be  added  by  superposition  to  yield  an  electron 
collision  frequency.  However,  just  adding  this  term  will  only 
work  well  only  in  the  fully,  but  singly,  ionized  regime.  In 
addition,  if  absolute  values  of  conductivity  are  desired,  whether 
classical  or  effective,  the  values  mav  still  be  off  by  as  much  as 
a  factor  of  2. 

If  only  the  relationship  between  anomalous  and  classical 
transport  is  desired  in  all  regimes,  not  absolute  values,  a 
compromise  between  the  Bose  method  and  Caldo  method  can  be  used 
to  reduce  the  work  involved  in  the  Bose  method.  If  the  electron- 
electron  collisions  are  accounted  for  as  explained  in  the  above 
paragraph,  the  mean  free  path  formula  for  conductivity  can  be 
used  and  produce  a  good  deal  of  accuracy  in  determining  the 
effect  of  anomalous  effects  on  transport.  However,  the  mole 
fractions  of  neutrals  and  multiply  ionized  ions  must  be 


93 


calculated  by  some  means.  As  MPD  thrusters  almost  always  operate 
in  the  fully  ionized  regime,  neutrals  can  be  neglected  for 
practical  comparisons.  The  multiply  ionized  mole  fractions  are 
necessary  because  they  are  directly  linked  to  increases  in 
electron  mole  fractions  at  higher  temperatures. 

Unfortunately,  if  absolute  values  of  electrical 
conductivity  are  desired  as  well  as  determining  the  effects  of 
anomalous  transport,  the  method  used  in  this  paper  is  probably 
reguired.  It  takes  a  good  deal  more  time  to  calculate  results, 
as  well  as  being  difficult  to  understand,  but  it  accounts  for 
almost  all  interactions  while  providing  fairly  accurate  results 
for  both  classical  and  effective  electrical  conductivity. 

If  some  slight  modifications  are  made  to  this  method,  the 
effects  of  anomalous  transport  on  other  transport  properties 
which  involve  electron-ion  interactions,  such  as  thermal 
conductivity  and  diffusion,  could  probably  be  accounted  for 
fairly  easily.  This  is  because  the  Bose  method  calculates  all 
transport  properties  based  on  collision  cross  sections,  and  the 
cross-sections  for  electron-ion  collisions  have  already  been 
modified  for  anomalous  effects  to  find  electrical  conductivity. 
The  calculation  of  other  transport  properties  was  not  attempted 
in  this  paper  because  one  of  the  major  aims  was  to  compare  the 
effects  of  the  method  of  calculation  on  determining  the  effects 
of  anomalous  transport,  and  this  author  has  not  seen  any  work 
done  by  others  regarding  the  effect  of  anomalous  transport  on 
transport  properties  other  than  electrical  conductivity  in  an  MPD 
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thruster.  However,  for  some  quantities  which  rely  on  collisions 
transferring  energy  rather  than  momentum,  a  different  curve  fit 
must  be  used.  Such  a  curve  fit  is  also  shown  in  Caldo  and 
Choueiri's  paper  [28,29],  so  modification  of  this  method  would 
not  be  too  difficult.  Using  this  curve  fit,  they  have  examined 
ion  heating  due  to  anomalous  transport.  To  make  these  type  of 
calculations  may  be  something  for  the  interested  researcher  to 
attempt . 

To  have  more  confidence  in  the  results  from  the  method  used 
in  this  paper,  it  would  be  necessary  to  get  cross  section  data 
while  anomalous  effects  (i.e  microinstabilities)  are  present, 
experimentally.  Although  microinstabi 1 i t ies  were  not  seen  at  the 
higher  power  levels  examined  in  this  experimental  work,  others 
[17]  at  Princeton  university  have  seen  them  at  lower  power 
levels,  and  Caldo  and  Choueiri  have  predicted  where  they  should 
occur  theoretically.  Collision  cross  sections  to  account  for 
anomalous  effects  should  be  measured  in  the  regions  shown  by  the 
Princeton  work,  both  experimental  and  theoretical,  in  order  to 
more  accurately  predict  anomalous  transport  properties. 

Lastly,  the  means  of  ionization  in  the  MPD  plasma  needs  to 
be  determined.  If  it  appears  to  occur  according  to  classical 
methods,  no  further  modification  to  this  method  other  than  what 
is  mentioned  above  should  be  necessary.  However,  if  anomalous 
ionization  occurs,  the  portion  of  the  Bose  method  which 
calculates  the  quasi -equilibrium  species  compositions  must  be 
modified  to  account  for  anomalous  ionization. 
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APPENDIX  A 
Faci 1 i ty 

The  pulsed  mul t i -megawatt  MPD  facility,  known  as  chamber  *2 
at  the  United  States  Air  Force  Philips  Lab-Electric  Propulsion 
Laboratory,  consists  of  a  power  supply  and  conditioning  system, 
vacuum  pumping  equipment,  propellant  flow  system,  and  a  pulsed 
magnetoplasmadynamic  (MPD)  thruster. 

A. I  Vacuum  Equipment: 

The  chamber  is  a  2.438  m  (8  ft)  diameter  by  3.658  m  (12  ft) 
long  stainless  steel  cylinder.  It  has  one  door  and  6  view  ports. 
Four  of  these  viewports  are  transparent,  allowing  for  visual 
observation  of  the  thruster  both  directly  and  with  the  aid  of  a 
mirror  placed  inside  the  chamber.  The  other  two  are  occupied  by 
interfaces  which  equipment  is  plugged  into.  A  diagram  of  the 
setup  is  shown  in  figure  A.l  [34]. 

The  equipment  to  produce  the  near- vacuum  in  the  chamber  are 
a  Stokes  412H-10  Mechanical  Pump,  a  Roots  615  RGS  Blower,  and  two 
Varian  0185  .3048  m  (12  inch)  diffusion  pumps.  These  pumps  can 
reduce  the  chamber  pressure  to  0.006666  Pa  (5  x  10*:  Torr). 

After  .6  g  of  argon  is  exhausted  by  the  thruster,  the  chamber  can 
return  to  this  pressure  in  5  minutes  [34], 

A. 2  Propellant  Feed  system: 

The  propellant  system  consists  of  a  9.232  m^  (326  ft^)  argon 
T-bottle  and  regulator,  with  .793  m^  (28  ft^)  bottles  of 
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hydrogen,  helium,  neon,  xenon,  and  krypton  also  plumbed  into  the 
system  [1].  The  propellant  lines  feed  into  a  .02367  or  (.8359 
ftJ)  plenum  located  outside  the  chamber.  A  thermometer  at  the 
plenum  takes  the  temperature  of  the  argon.  The  pressure  drop  per 
gas  pulse  is  also  measured  here,  giving  the  mass  flow  at  the 
plenum.  A  solenoid  valve  allows  the  gas  to  enter  the  thruster 
when  firing. 


A. 3  Power  Supply  and  Conditioning: 

A  Del  Electronics  Corp.  HPS-1-8000-3  power  supply  provides 
up  to  8  kW  of  Power  to  a  Pulse  Forming  Network  (PFN).  The  PFN  is 
a  ten-section  LC  (Inductor-Capacitor)  Network  with  a  nominal 
output  impedance  of  0.01  A.  Each  section  consists  of  three  2000 
pF  Maxwell  33800  capacitors  connected  in  parallel.  With  the 
addition  of  a  5  turn,  0.53  pH  inductor  in  each  section,  the 
network  will  produce  a  1  millisecond  current  pulse  at  up  to  40  kA 
and  400  V  to  the  thruster  assuming  a  perfectly  matched  load  [34] . 
The  power  supply  is  isolated  from  the  PFN  once  the  PFN  has 
acquired  the  desired  voltage  level.  A  schematic  diagram  of  the 
pulse  forming  network  is  shown  in  figure  A. 2. 

A  high  voltage  power  supply  (low  current)  provides  the 
voltage  for  the  spark  trigger.  Typically,  it  is  operated  at 
between  700-1000  volts.  It  is  connected  to  a  tungsten  wire  and 
the  cathode  of  the  thruster  by  means  of  a  mechanical  high  voltage 
relay  when  the  thruster  is  fired. 
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A. 4  Timing  and  control: 

Timing  and  control  are  accomplished  through  the  use  of  a 
pushbutton  control  panel,  Tektronix  PEP-301  IBM  compatible 
microcomputer,  and  two  Rapid  Systems  4000  timing  boxes.  The  PFN 
is  connected  to  the  high  voltage  power  supply  by  a  relay 
controlled  through  pushbuttons  on  the  control  panel.  A 
potentiometer  is  used  to  control  the  power  supply  voltage,  and  in 
turn  the  PFN  voltage.  Once  the  PFN  is  at  the  desired  level,  the 
power  supply  is  disconnected. 

The  microcomputer  has  two  pulse  types  stored  in  a  file,  one 
for  the  gas  pulse,  and  another  for  the  spark  trigger.  The  gas 
pulse  is  eighty  milliseconds  long,  and  the  spark  trigger  pulse 
lasts  two  seconds.  The  spark  trigger  also  incorporates  a  20 
millisecond  delay  after  the  gas  pulse  is  initiated.  These  files 
are  downloaded  into  the  timing  boxes  which  produce  the  desired 
pulses.  Both  pulses  are  initiated  when  the  fire  button  is 
depressed. 


A. 5  Thruster 

The  thruster  itself  is  a  self-field  magnetoplasmadynamic 
(MPD)  thruster.  The  annular  copper  anode  has  a  10.2  cm  (4  inch) 
outer  diameter  and  a  7.62  cm  (3  in)  inner  diameter,  with  a  length 
from  injection  plate  to  end  of  3.81  cm  (1.5  in).  The  boron- 
nitride  injection  plate  fits  snugly  against  the  anode  surface, 
and  has  a  half  inch  radial  displacement  from  the  1.27  cm  (0.5  in) 
diameter,  thoriated  tungsten  cathode.  The  tip  of  the  cathode  is 
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.635  cm  (.25  in)  from  the  end  of  the  anode  in  an  axial  direction. 
Some  photographs  of  the  thruster  are  shown  in  figure  A. 3. 

The  entire  thruster  is  imbedded  in  a  plexiglass  mounting 
which  may  swing  in  the  axial  direction  or  remain  fixed  depending 
on  the  operator's  preferences.  An  accelerometer  may  be  imbedded 
in  the  plexiglass  to  make  measurements  as  well.  The  propellant 
lines  are  inserted  into  the  rear  of  the  thruster.  The  PFN  is 
connected  to  the  thruster  by  means  of  a  copper  cross  attached  to 
the  cathode  with  a  screw,  and  loops  of  copper  wire  pressed 
against  the  anode  by  a  metal  p^ate,  attached  to  the  anode  with 
screws.  In  addition,  a  tungsten  wire  protrudes  through  the 
injection  plate  to  provide  the  high  voltage  spark  trigger  for 
thruster  ignition.  It  forms  the  spark  with  the  cathode. 

A. 6  Measurement  Devices 

Various  forms  of  instrumentation  were  used.  Chamber 
pressure  is  measured  with  a  Varian  843  vacuum  ionization  gauge. 
Pearson  Electronics  301x  Pulse  transformers  measure  the  current 
flowing  through  the  thruster.  Voltage  across  the  thruster  is 
measured  by  1000:1  voltage  probes  attached  to  the  PFN.  A 
Tektronix  DSA-601  digital  signal  analyzer  and  a  Tektronix  11403 
digitizing  oscilloscope  can  be  used  to  capture  data  from  the 
instrumentation.  Signal  conditioning  for  the  various  measurement 
devices  is  aided  by  the  use  of  Tektronix  AM501  operational 
amplifiers . 
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A  Klinger  stepper  motor  assembly  allows  positioning  of 
various  measurement  devices  inside  of  the  chamber.  The  arm  is 
level  with  the  thruster  axis  and  can  be  translated  and  rotated. 
Each  translational  step  is  2.54  x  10"'  m,  and  each  step  in 
rotational  mode  is  .001  deg. 

A. 7  Method  of  Operation 

The  firing  procedure  begins  by  charging  the  PFN  through  a 
connection  with  the  power  supply  until  a  desired  PFN  voltage  is 
reached,  at  which  time  the  PFN  is  then  disconnected  from  the 
power  supply.  The  voltage  present  in  the  PFN  is  directly  applied 
to  the  anode  and  cathode  of  the  thruster.  The  operator  also 
ensures  that  the  pumps  have  brought  the  chamber  pressure  down  to 
approximately  5  x  10*^  torr. 

When  the  fire  button  is  depressed,  the  timing  boxes  initiate 
their  pulses.  The  gas  pulse  signal  opens  the  solenoid  valve,  and 
twenty  milliseconds  later,  the  spark  pulse  is  initiated.  The 
resultant  gas  pulse  has  exponential  growth  and  decay.  With  the 
additional  lag  time  resulting  from  the  mechanical  relay  between 
the  high  voltage  power  supply  and  thruster,  the  spark  occurs 
when  the  gas  pulse  has  reached  steady  state.  This  spark  ionizes 
the  plasma  initially,  allowing  the  higher  currents  from  the  PFN 
to  flow  between  the  anode  and  cathode. 

Current  measurements  are  taken  using  the  Pearson  coils. 

These  are  low  freguency  response  devices;  therefore,  most  noise 
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is  filtered  out  of  the  current  trace.  However,  the  voltage 
probes  have  a  high  frequency  response,  and  onset  can  be  detected 
by  looking  at  the  voltage  fluctuations  in  the  steady  state  region 
of  the  voltage  trace. 


Network 


B-  1 


APPENDIX  B 
Magnetic  Probes 

In  order  to  calculate  the  lower  hybrid  frequency  of  the 
plasma,  it  is  necessary  to  know  the  magnetic  field  at  a 
particular  point.  To  accomplish  tnis,  two  different  types  of 
probes  were  used:  a  Hall  probe  and  an  induction  probe. 

B.l  Hall  Probe 

B.1.1  Theory 

The  Hall  generator  works  on  the  principle  known  as  a  Hall 
effect.  Consider  a  flat  plate  conductor  placed  in  a  magnetic 
field  with  the  normal  to  the  surface  of  the  conductor  aligned 
with  the  magnetic  field  lines.  If  a  current  is  simultaneously 
applied  through  the  plate  at  opposite  ends,  a  potential 
difference  is  created  in  the  plate  on  the  sides  parallel  to  the 
direction  of  current  flow. 

The  reason  for  this  is  fairly  simple.  The  force  generated 
by  the  current  flowing  through  the  magnetic  field  is 

F  »  IxB  (B.l) 

where  F  a  magnetic  force,  I  a  current  vector,  and  B  *  magnetic 
field  vector.  This  force  is  applied  to  the  charged  particles  in 
the  conductor.  However,  since  the  ions  do  not  move  easily,  the 
electrons  are  assumed  to  be  the  only  particles  that  move.  Since 
this  force  is  applied  perpendicular  to  the  current  path,  the 
electrons  are  forced  to  move  to  one  of  the  edges  of  the  plate 
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which  is  parallel  to  the  current;  they  cannot  travel  any  further 
from  the  path  than  the  edge  of  the  plate.  The  predominance  of 
electrons  on  one  side  and  the  absence  of  those  electrons  on  the 
other  creates  a  potential  difference  on  those  sides  of  the  plate. 
However,  the  potential  difference  will  make  those  electrons  want 
to  travel  back  to  the  other  side  of  the  plate.  The  only  force 
preventing  that  is  the  magnetic  force.  At  equilibrium 
conditions,  the  magnetic  force  will  equal  the  electromotive  force 
produced  by  the  potential  difference.  Thus,  the  potential 
difference  is  directly  proportional  to  the  magnetic  field  in  the 
perpendicular  direction.  This  potential  difference  can  be 
measured  in  many  ways,  but  if  it  is  measured  by  an  oscilloscope, 
a  trace  of  magnetic  field  vs.  time  can  be  acquired  as  well. 

The  Hall  probe  has  been  used  successfully  in  steady  state 
MPD  plasmas  at  the  kW  level.  Its  primary  advantage  is  that  the  B 
field  strength  may  be  calculated  directly  from  the  observed 
voltage.  However,  calibration  is  required  in  a  known  magnetic 
field  in  order  to  use  it  later  in  fields  of  unknown  strength. 

B.1.2  Experimental  Use 

In  the  experiments  discussed  in  chapter  III,  a  Bell  BH-2C5 
Hall  generator  was  used.  A  factory  calibration  came  with  the 
probe  (about  .91  pV/Gauss  at  100  mA  current),  but  a  calibration 
performed  in  the  laboratory  using  a  Helmholtz  Coil  operated 
between  10  and  50  gauss  showed  a  relationship  of  .87  pV/Gauss  at 
100  mA.  It  should  be  noted  that  one  disadvantage  of  the  Hall 
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probe  is  that  the  probe  sensitivity  decreases  with  temperature. 
One  reason  for  the  discrepancy  in  calibrations  is  that  isothermal 
operation  of  the  probe  was  not  possible  in  the  laboratory;  the 
current  running  through  the  probe  heated  it  slightly  above  the 
temperature  at  which  it  was  calibrated  at  the  factory.  This 
heating  would  decrease  the  sensitivity  slightly,  which  is 
accounted  for  in  the  reduction  of  slope  of  the  calibration 
relation. 

The  first  attempts  to  measure  the  magnetic  field  strength  in 
the  azimuthal  direction  involved  the  use  of  a  Hall  probe. 

Another  disadvantage  of  the  Hall  probe  is  its  extremely  low 
signal  voltage,  on  the  order  of  microvolts.  Unfortunately,  the 
electrical  noise  resulting  from  the  start-up  of  the  thruster  was 
an  order  of  magnitude  higher  than  the  signal  which  was  observed 
during  calibration  of  the  Hall  probe. 

To  ensure  that  the  unexpected  signal  was  due  to  extraneous 
noise  rather  than  an  actual  magnetic  field  fluctuation,  the  probe 
was  removed  from  its  location  in  the  near-plume  of  the  thruster 
and  placed  approximately  1.5m  directly  to  the  side  of  the 
thruster.  At  this  location,  the  magnetic  field  should  have  been 
almost  non-existent.  The  thruster  was  then  fired  and  the  signal 
was  compared  to  the  signal  where  the  noise  was  first  observed. 

Almost  the  same  signal  resulted.  The  noise  eventually 
disappeared,  but  the  settling  time  of  the  noise  was  on  the  order 
of  a  millisecond,  which  happened  to  be  the  length  of  the  steady 
state  portion  of  the  current  pulse.  Thus,  only  the  ending 
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transient  of  the  magnetic  field  signal  was  discernible,  which  was 
of  no  use  in  determining  the  magnetic  field  strength  during  the 
steady-state  portion  of  the  pulse. 

In  an  attempt  to  eliminate  this  electrical  noise,  the  probe 
was  placed  outside  of  the  chamber  and  the  thruster  was  fired. 

The  large  noise  spike  occurred,  but  it  damped  out  well  within  the 
duration  of  the  current  pulse.  The  noise  was  at  ar.  acceptable 
level  during  the  latter  half  of  the  steady-state  region  of  the 
pulse.  Because  the  chamber  was,  in  essence,  a  large 
electrostatic  shield,  it  was  hoped  that  additional  shielding  of 
the  electrical  leads  of  the  probe  could  reduce  the  noise  to  an 
acceptable  level.  This  additional  shielding  took  the  form  of  a 
layer  of  aluminum  foil  connected  to  ground,  wrapped  around  all  of 
the  coaxial  cables.  In  addition,  wires  were  twisted  in  an 
attempt  to  eliminate  ground  loops,  which  could  also  introduce 
spurious  magnetic  field  signals.  Although  these  actions  reduced 
noise  slightly,  the  effects  were  negligible;  these  actions  still 
did  not  reduce  noise  to  an  acceptable  level. 

All  of  the  leads  were  attached  to  a  the  interior  wire  of  a 
coaxial  cable.  It  might  help  to  acquire  trlaxial  cable  and 
attach  the  leads  to  both  the  interior  cable  and  the  first  coaxial 
layer  of  wire,  with  the  third  layer  used  as  shielding.  In 
addition,  about  10  inches  of  the  probe  was  unshielded  except  for 
a  twisting  of  the  voltage  and  current  leads  about  themselves,  but 
this  problem  could  not  be  easily  eliminated.  A  better  probe 
design  might  try  to  account  for  the  need  for  shielding  of  the 
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leads  which  are  directly  attached  to  the  Hall  generator. 

It  appears  that  this  noise  is  not  produced  by  the  quasi¬ 
steady  operation  of  the  thruster,  but  instead  is  caused  by  the 
high  voltage  spark  which  breaks  down  the  gas  to  allow  current  to 
flow  in  the  quasi-steady  pulse.  If  a  less  noisy  means  of  start¬ 
up  were  devised,  it  is  conceivable  that  the  Hall  probe  could  be 
used  satisfactorily  in  a  pulsed  device  even  without  the 
additional  shielding.  The  premise  that  the  spark  ignition  of  the 
engine  is  what  produces  the  electrical  noise  seems  to  explain  why 
Hall  probes  have  been  used  to  measure  magnetic  fields 
satisfactorily  in  steady  state  MPD  devices;  the  time  at  which 
magnetic  fields  are  measured  are  typically  long  after  the  start¬ 
up.  These  transient  effects  of  the  spark  trigger  may  also 
explain  why  this  author  has  not  seen  any  mention  of  Hall  probes 
being  used  in  pulsed  MPD  devices;  induct ion- type  magnetic  probes 
are  the  only  type  mentioned.  Since  the  Hall  probe  did  not 
produce  any  usable  results,  an  induction-type  probe  with  an 
analog  integrator  circuit  was  constructed. 

B.2  Induction  Probe 

B . 2 . 1  Theory 

An  induction  probe  is  very  simple  in  its  construction;  it  is 
merely  a  small  cylindrical  coil  of  wire.  In  fact,  this  coil  is 
so  small,  it  is  assumed  that  the  magnetic  field  strength  does  not 
vary  over  its  cross  section.  As  a  result,  the  voltage  induced  by 
the  introduction  of  a  magnetic  field  through  the  probe  can  be 
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represented  by  [48] 

V.  =  iC)  (  B.  2 ) 

d~ 

where  V.  2  induced  voltage,  Aa.-<  =  effective  cross  sectional  area 
of  the  coil,  and  3X  2  magnetic  field  perpendicular  to  the  coil. 
In  this  relation,  the  voltage  is  proportional  to  the  time 
derivative  of  B,  rather  than  directly  proportional  to  B.  As  a 
result,  the  signal  must  be  passed  through  an  integrator  to  get  B 
as  a  function  of  voltage.  It  is  not  necessary  to  calculate  the 
effective  cross  sectional  area  of  the  probe;  if  the  probe  is 
calibrated  in  a  known  magnetic  field  while  the  signal  is  passed 
through  the  integrator,  a  calibration  curve  for  B(V,)  can  be 
produced.  In  fact,  this  method  is  preferred.  Not  only  is  it 
hard  to  get  accurate  measurements  of  the  probe  diameter  for  small 
cross-sectional  areas,  but  the  presence  of  the  integrator  may 
alter  the  signal  slightly.  As  long  as  the  signal  is  calibrated 
with  the  integrator  in  the  circuit,  anomalies  due  to  the 
integrator's  presence  will  be  accounted  for. 

In  general,  induction  probes  can  produce  voltages  much 
higher  than  the  electrical  noise  level  of  the  plasma  [48],  but 
the  actual  level  of  voltage  is  a  function  of  the  effective  area 
of  the  coil  and  dB/dt.  One  way  to  increase  the  signal  level  is 
to  increase  the  effective  area  of  the  coil  cross  section.  This 
can  be  done  in  one  of  two  ways.  First,  the  coil  size  can  be 
increased.  Unfortunately,  if  this  route  is  attempted,  the 
spatial  resolution  will  decrease.  The  second  way  is  to  increase 
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the  number  of  turns  of  wire  about  the  core.  However,  one  effect 
of  adding  more  turns  of  wire  is  to  decrease  the  frequency 
response  of  the  probe. 

B.2.2  Experimental  Use 

The  induction  probe  used  in  the  measurements  in  chapter  III 
was  1/16"  long  by  1/16"  dia,  with  75  turns  of  wire.  This  design 
produced  a  fairly  high  signal  to  noise  ratio  with  a  frequency 
response  on  the  order  of  100s  of  MHz.  The  frequency  response  of 
the  probe  was  not  very  important,  though,  because  only  the  mean 
field  strength  in  the  steady  state  region  of  the  pulse  was 
required.  In  addition,  the  integrator  used  an  operational 
amplifier  which  had  an  effective  cutoff  frequency  of  10  MHz. 

Thus,  the  circuit  was  limited  in  response  by  the  operational 
amplifier,  not  the  probe. 

Unfortunately,  this  high  frequency  signal  attenuation  causec 
some  error  in  integration.  Since  very  high  frequency 
fluctuations  and  spikes  are  attenuated,  the  integrated  signal 
will  neglect  these  fluctuations,  and  the  level  will  be  either 
lower  or  higher  than  what  is  actually  seen.  Although  this  has 
the  benefit  of  eliminating  noise,  which  theoretically  would 
average  out  to  zero  over  the  integration,  there  may  be  some  very 
fast  rise  times  in  the  magnetic  field  strength  in  the  initial 
portion  of  the  current  pulse  which  may  be  attenuated  or 
eliminated  entirely  producing  a  small  error  in  measurement. 

Obviously,  the  stronger  the  magnetic  field,  the  stronger  the 
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signal.  On  the  other  hand,  the  time  it  takes  for  the  B-field  to 
reach  that  strength  also  comes  into  play.  Based  on  this 
dependence,  one  can  see  that  it  is  not  really  possible  to  apply 
an  induction  type-probe  in  a  steady  state  device;  theoretically, 
there  is  no  change  in  magnetic  field  with  time.  However,  it  is 
ideal  for  a  pulsed  discharge.  As  the  rise  time  of  the  current  in 
the  quasi-steady  thruster  used  in  chapter  III  is  on  the  order  of 
microseconds,  dB/dt  is  fairly  large  at  most  locations  in  the 
thruster  and  plume.  However,  this  dependence  on  dB/dt  caused 
some  problems  in  the  probe  calibration. 

Initially,  calibration  was  attempted  in  the  Helmholtz  coil 
as  was  done  for  Hall  probe  calibration.  However,  the  power 
supply  which  powered  the  Helmholtz  coil  proved  to  be  a  limiting 
factor.  First,  the  power  supply  could  only  operate  up  to  5  A, 
which  limited  the  produced  magnetic  field  to  50  Gauss  based  on 
the  calibration  for  the  Helmholtz  coil.  The  magnitude  of  the 
known  field  posed  no  problem;  plume  measurements  are  within  this 
range,  and  the  signal  was  expected  to  be  larger  closer  to  the 
thruster.  Unfortunately,  the  limiting  factor  was  the  time  it 
took  for  the  power  supply  to  reach  the  desired  current  level:  on 
the  order  of  several  milliseconds,  depending  on  the  current.  As 
this  is  about  3  orders  of  magnitude  slower  than  the  rise  time 
inside  the  thruster,  which  is  what  the  probe  was  designed  for,  a 
signal  was  not  discernible  for  calibration  in  the  Helmholtz  coll. 

To  overcome  this  problem,  the  probe  was  calibrated  in  the 
thruster  itself.  To  do  this,  it  was  assumed  that  if  the  probe 
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was  placed  at  the  back  plate  of  the  thruster  near  the  anode,  all 
of  the  current  going  to  the  thruster  would  be  enclosed  by  a 
circle  corresponding  to  the  radius  of  the  probe  location. 

Although  this  assumption  is  very  reasonable,  there  is  an  error  on 
the  order  of  about  1%  on  the  enclosed  current.  To  find  the 
magnetic  field  strength  at  the  probe  location,  the  following 
formula  was  used: 

( B.  3 ) 

6  2n  R 

where  Se  h  magnetic  field  in  azimuthal  direction  and  R  =  the 
radial  distance  from  the  thruster  axis.  The  magnetic  field  was 
varied  by  changing  the  input  current  to  the  thruster.  The 
resultant  calibration  curve  was  1.869  mV/Gauss,  significantly 
larger  than  that  for  the  Hall  probe. 

Although  this  method  had  some  error  involved  in  the  "known" 
magnetic  field,  it  also  had  some  benefits  that  calibration  in  the 
Helmholtz  coil  did  not.  First  and  foremost,  calibration  of  the 
probe  through  almost  the  entire  range  of  magnetic  field  strengths 
was  possible.  The  Helmholtz  coil  produced  very  accurate 
calibration  data  up  to  50  Gauss,  but  much  larger  magnetic  fields 
were  expected.  If  the  calibration  curve  was  completely  linear, 
the  Helmholtz  data  would  have  been  adequate,  but  if  it  wasn't, 
there  would  have  been  significant  error  in  the  range  of  high 
magnetic  field  strengths.  By  calibrating  throughout  the  entire 
range  of  expected  field  strengths,  this  source  of  error  was 
eliminated.  The  additional  benefit  is  that  the  signal  strength 
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in  comparison  to  the  noise  level  could  be  observed  prior  to 
actual  magnetic  field  measurements. 
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APPENDIX  C 
Langmuir  Probes 

The  triple  Langmuir  probe  is  a  measurement  device  which  can 
be  used  to  determine  number  density  and  electron  temperature  in 
the  plasma.  Before  going  into  the  configuration  which  was  used 
in  the  experiments,  it  is  important  to  discuss  the  theory  of 
Langmuir  probes  in  general. 

C.l  Theory 

The  simplest  form  of  a  Langmuir  probe  is  simply  an  exposed 
wire.  In  order  to  measure  electron  temperature  and  number 
density,  the  probe  must  be  "floating".  In  other  words,  no  net 
current  is  transferred  between  the  probe  and  the  plasma.  Since 
current  is  the  net  movement  of  charge  into  the  wire,  in  order  to 
make  the  probe  "float",  the  number  of  electrons  and  ions  absorbed 
into  the  probe  must  equal  the  amount  leaving  the  probe.  At  first 
glance,  it  would  appear  that  this  can  be  accomplished  be  setting 
the  probe  at  a  potential  which  is  the  same  as  the  plasma. 

Although  doing  so  does  not  produce  a  floating  probe,  it  is  useful 
to  understand  what  would  happen  in  such  a  case.  With  no 
potential  difference,  and  consequently,  no  electric  field 
present,  electrostatic  forces  would  not  cause  a  predominance  of 
one  type  of  charged  particle  at  the  probe  surface.  Logically, 
the  only  particles  which  would  impact  on  the  probe  surface  would 
result  from  the  random  thermal  motion  of  the  plasma.  Since  the 
number  of  particles  which  would  impact  on  the  probe  is  a  function 
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of  the  number  of  particles  present  locally  in  the  probe's 
vicinity,  and  also  a  function  of  their  temperature,  the  measured 
current  and  potential  can  indicate  the  local  number  density  and 
electron  temperature. 

Unfortunately,  there  is  a  complication.  Since  electrons 
are  so  much  lighter  than  ions,  for  a  given  temperature  in  an 
isothermal  plasma,  more  electrons  than  ions  will  impact  the  probe 
surface.  This  is  why  merely  setting  the  probe  at  the  plasma 
potential  will  not  result  in  a  floating  probe.  In  order  to  make 
it  float,  the  probe  must  be  biased  at  a  negative  potential  with 
respect  to  the  plasma;  this  negative  potential  will  repel  enough, 
electrons  through  electrostatic  forces  to  allow  an  equal  number 
of  electrons  and  ions  to  impinge  on  the  probe  surface. 

However,  the  presence  of  a  sheath  leads  to  further 
complications  of  the  situation.  One  of  the  properties  of  a 
plasma  is  to  maintain  macroscopic  charge  neutrality,  and  the 
formation  of  sheaths  about  potential  sources  is  the  mechanism  the 
plasma  uses  to  accomplish  this  end.  The  sheath  is  not  a  problem; 
particles  move  freely  outside  the  sheath  region,  but  the  electric 
fields  present  inside  the  sheath  prevent  too  many  electrons  from 
hitting  the  surface.  This  results  in  the  plasma  remaining 
largely  undisturbed  by  the  presence  of  the  probe,  while  still 
allowing  the  probe  to  measure  number  density  and  electron 
temperature.  The  problems  arise  because  the  measured  current, 
specifically  the  number  of  ions  hitting  the  probe,  depends  on  the 
size  of  the  sheath  relative  to  the  size  of  the  probe.  Different 
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theories  must  be  used  to  relate  current  to  electron  temperature 
and  number  density  depending  on  the  relative  sizes  of  the  mean- 
free  path,  sheath  thickness,  and  probe  size. 

In  the  MW  level  MPD  thruster,  the  number  densities  and 
electron  temperatures  have  been  shown  to  produce  mean  free  paths 
on  the  order  of  a  mm,  and  debye  lengths  around  a  pm.  Since  it  is 
very  difficult  to  get  a  wire  which  is  on  the  order  of  a  debye 
length  which  will  retain  any  durability,  the  probe  theory  which 
is  used  is  based  upon  the  assumption  of  a  thin,  collisionless 
sheath.  In  thin  sheath  theory,  the  assumption  can  be  made  that 
the  ion  current  is  independent  of  sheath  potential. 

For  the  triple  probe,  all  of  the  probes  are  biased  with  a 
negative  potential  relative  to  the  plasma.  Figure  C.l  shows  the 
configuration  of  the  triple  probe  the  key  potential  differences. 
Electrode  2  is  floating,  while  there  is  a  current  path  flowing 
into  electrode  1  and  out  of  electrode  3.  No  net  charge  is  lost 
from  the  plasma  because  the  charge  that  is  lost  by  electrode  1  is 
regained  at  electrode  3.  Using  Kirchoff's  Law,  equations  can  be 
written  for  the  three  electrodes. 
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where 


Js0  =  electron  current  at  sheath  edge  =  nv 

J:  =  ion  flux  at  electrode  surface  =  exp(-0.5)en9 

For  these  equations,  Vj  is  the  potential  of  electrode  j ,  V?  * 
plasma  potential,  %j  5  dimensionless  potential  of  electrode  j, 
and  J;  is  based  on  Bohra’s  sheath  criteria. 

Using  this  system  of  equations,  Te  and  n^  can  be  solved  for 
if  the  following  assumptions  are  made: 

1)  Jj  is  independent  of  in  the  thin  sheath  limit. 

2)  A,  .  Aj 

3)  is  set  by  the  experimenter. 

Solving  this  system  of  equations  yields  the  following  relations 
between  n^T^V^,  I,  and  Aj : 
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■1  =  exp 


XT*  ) 


{  C  .  4  ) 


^S*  exp  ;  0 . 5 ) 


n,  = 


KTe 

\j  W. 


( C .  5 ) 


‘  ll-exp(x^'X^) ) 


Equation  (C.4)  must  be  solved  implicitly  to  find  T 


C.2  Experimental  Use 

The  probe  used  in  the  experiments  in  chapter  III  consisted 
of  three  0.125  mm  dia  tungsten  wires  aligned  parallel  to  each 
other  having  an  exposed  length  of  4  nun.  liie  ends  of  the  wires 
were  aligned  in  axial  position  as  well.  Insulation  of  the  wires 
was  provided  by  alumina  tubes  of  0.25  mm  inner  dia  and  1  mm  outer 
diameter.  EPOTEK  high  temperature  epoxy  was  used  to  attach  the 
wires  to  the  alumina  tubes,  fixing  the  area  of  exposed  wire  at  Aj 
=  0.5n  nun1.  The  alumina  tubes  were  attached  to  a  4  hole  alumina 
holder  with  1.5  mm  dia  holes  and  6  mm  outer  diameter  by 
Ceramabond  671  ceramic  paste.  In  addition,  the  tungsten  wires 
were  attached  to  coaxial  cables  at  the  base  of  the  alumina  holder 
to  provide  electrostatic  shielding,  and  to  allow  the  use  of  BNC 
connectors  for  the  vacuum  feedthrough.  The  probe  itself  was 
housed  in  a  plexiglass  holder  attached  to  a  translational  X-Y 
table,  where  the  tip  of  the  probe  was  approximately  19-20  cm  from 
the  probe  holder.  The  inter  electrode  spacing  was  3  mm  between  1 
and  2,  1 . 5  mm  between  2  and  3,  and  2.5  ram  between  1  and  3  in  a 
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triangular  configuration. 

The  probe  was  designed  to  meet  the  following  criteria. 
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where  r5  s  radius  of  the  probe,  ktl  a  ion-ion  mean  free  path,  xa 
end  effect  parameter,  Zi  3  charge  of  the  ion,  Lp  a  length  of  the 
probe,  and  S  a  inter-electrode  spacing.  Each  if  these  equations 
allow  one  or  more  assumptions  to  be  made  if  they  are  satisfied. 
(C.6)  is  the  condition  which  is  necessary  to  use  thin  sheath 
theory,  (C.7)  assures  that  the  sheath  is  collisionless,  (C.8) 
allows  end  effects  to  be  ignored,  (C.9)  assures  that  ths  sheaths 
surrounding  the  electrodes  do  not  significantly  interfere  with 
each  other,  (C.10)  makes  sure  that  the  presence  of  a  magnetic 
field  does  not  significantly  hinder  the  movement  of  electrons  to 
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the  probe,  and  (C.ll)  insures  that  a  finite  electron  temperature 
can  be  calculated  from  V.;,  and  that  the  electron  flux  is 
negligible  compared  to  the  ion  flux.  The  probe  described  above 
was  designed  to  meet  these  conditions  based  on  expected 
conditions  in  the  thruster  plume. 

It  must  be  noted  that  the  underlying  assumption  which  allows 
the  use  of  thin  sheath  theory  is  that  the  ion  flux  is  independent 
of  X-  •  However,  the  functional  dependence  of  J.  depends  on 
several  factors.  The  conditions  in  (C.6-11)  allow  some  of  these 
functional  dependencies  to  be  assumed  negligible,  but  x,  can 
still  come  into  play.  The  ratio  of  ion  velocity  which  is 
perpendicular  to  the  probe  axis  over  the  most  probable  ion 
thermal  velocity  also  plays  a  role.  In  addition,  electron 
retarding  current  can  be  significantly  affected  by  magnetic 
fields,  electron  drift,  and  electron-ion  collisions  [35]. 

However,  in  the  experimental  work,  Vjj  was  set  rather  high  so 
that  the  electron  flux  at  electrode  3  is  negligible  compared  to 
the  ion  flux.  This  allows  the  formation  of  an  explicit  . relation 
for  V,:. 
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KT ' 
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where  I.}  s  ion  current  to  electrode  j ,  Aj  s  collection  area  of 
electrode  j,  and  AVconr  m  contact  potential  difference.  The 
first  term  reduces  to  In  2  in  the  thin  sheath  limit,  the  second 
term  takes  into  account  the  effect  of  unequal  probe  areas,  the 
third  takes  into  account  deviation  of  the  electron  flux  from 
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AjJ^gSX p(~Xj)  .and  the  last  term  takes  into  account  the  contact 
potential  difference  when  the  probe  is  contaminated.  To  use  the 
simplified  theory,  the  assumption  is  made  that  the  last  2  terms 
are  negligible  for  a  clean  probe,  and  the  probe  was  designed  so 
that  electrodes  have  equal  areas.  This  means  the  second  term  is 
negligible  as  well.  Thus,  the  simple  relation  eV^/KTg  =  In  2, 
from  which  electron  temperature  can  be  directly  calculated,  is 
all  that  remains. 


C.3  Laf ramboise ' s  Exact  Theory 
For  locations  where  X3  >3  and  r,/  kc  is  between  5  and  100, 
Laf ramboise ' s  theory  accompanied  by  the  Peterson-Talbot  curve 
fits  [36]  may  be  used.  This  theory  does  take  into  account  sheath 
potential  and  the  results  are  considered  exact.  The  curve  fits 
are: 

Ji(Xj)  =  iVX;]“  (C.  13) 

where 


This  theory  also  makes  the  assumption  that  only  one  ionized 
species  exists  and  end  effects  and  sheath  interactions  are 
negligible.  To  make  this  probe  theory  work,  the  following 
equations  must  be  solved  simultaneously. 
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Implicit  Relation  for  electron  temperature  : 


1-P 3^-i  -  *  U-P„( 
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v%>  1 

d*j 


and  xt  a  dimensionless  floating  potential 


From  the  Peterson-Talbot  curve-fit: 


2a 

5P+Xf 


(C.16) 


The  condition  based  on  equal  electron  and  ion  currents  for  a 
floating  probe: 


[£P+Xf] 


■exp  (  ~2%()  =  0 


(C. 17 ) 


The  number  density: 


__ _  ‘'•prob* 

=  1  _ A, _  (C.  18) 

eN  kr*  (t5p+(x.«-Xdsi  +Xf]a-[Sp>xf]exp[xdi-Xdj3) 

When  these  equations  are  solved  numerically,  the  last  check  is  to 
make  sure  that  r./  Xr  is  consistent.  This  must  be  updated  to  a 
new  value  based  on  the  ne  and  Tg  calculated  above,  and  the 
calculations  must  be  repeated  until  convergence  occurs. 
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The  la?  point  which  must  be  made  is  that  the  exact  theory 
assumes  that  the  probe  is  aligned  with  the  flow.  A  non-alignment 
of  the  probe  with  the  flow  car.  lead  to  an  error  on  the  order  of  2 
in  both  Tg  and  n,.  However,  this  situation  will  also  result  no 
matter  which  theory  is  used,  because  sheath  theory  has  not 
progressed  far  enough  to  account  for  the  effects  of  an  asymmetric 
sheath.  In  general,  Laf ramboise ' s  theory  may  be  considered  exact 
between  r3/  Xc  of  5  and  100,  but  it  is  still  valid  for  values 
above  100,  although  there  is  a  small  error  involved.  Another 
method  may  be  used  above  this  area  as  well,  but  it  is  felt  that 
the  error  in  using  Laf ramboise * s  theory  is  well  within  the 
experimental  error  for  the  Triple  probe  in  the  MPD  temperature 
and  density  regime.  It  should  be  noted  that  for  x3  <  3, 
numerical  solutions  for  the  above  equations  do  not  converge  self- 
consistently  with  r,/  Xr  .  Thus  at  these  locations,  which 
corresponded  to  locations  inside  the  thruster  for  the  experiments 
performed  in  chapter  III,  simple  thin  sheath  theory  is  used, 
although  the  error  is  more  significant.  This  sacrifice  of 
accuracy  is  necessary  in  order  to  get  any  values  in  this  region. 
In  addition,  thin  sheath  theory  is  valid  to  much  lower  potentials 
because  it  assumes  ion  flux  is  independent  of  potential. 

Laf ramboise ' s  theory  could  be  used  in  this  region  of  the  thruster 
if  were  increased  to  a  much  higher  value,  but  this  was  not 
done  in  the  experiments  discussed  in  this  paper. 
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Figure  C.l  - 
potentials 


Schematic  diagram  of  Langmuir  probe  electrodes  and 
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APPENDIX  D 

Documentation  of  Computer  Code 
D.l  Introduction 

This  program  package  for  calculation  of  thermophysical  and 
transport  properties  of  a  two-temperature  plasma  was  written  by 
Dr.-Ing.  Tarit  Kumar  Bose,  Professor  of  Aeronautical  Engineering, 
Indian  Institute  of  Technology,  Madras,  India  with  the  help  of 
his  students  in  1987.  The  present  modification  of  this  code  to 
account  for  microinstabilities  was  accomplished  by  Erik  C. 

Bowman,  a  graduate  student  of  Dr.  S.N.B.  Murthy  in  the  School  of 
Aeronautical  and  Astronaut ical  Engineering  at  Purdue  University 
in  West  Lafayette,  Indiana. 

One  of  the  earlier  versions  of  the  program  was  written  in 
connection  with  the  doctoral  work  of  R.V.  Seeniraj,  and  the  files 
are  designated  with  the  words  SEEN  or  SEENIRAJ.  The  Seeniraj 
files  are  for  noble  gases  only,  but  simpler  forms  of  the  mixture 
properties  are  used;  for  atom-ion  collision,  however,  the 
charge-transfer  collision  data  are  used,  and  for  atom-atom 
collisions  the  exponential  repulsive  potential  data  (Monchik 
tables),  which  are  available  for  noble  gases,  are  used.  Later 
programs  may  be  used  for  the  calculations  for  more  general 
dissociated  and  ionized  gas  mixtures,  and  they  use  Lennard-Jones 
6-12  potential  for  both  atom-  atom  and  atom-ion  collisions.  They 
also  use  more  complex  transport  property  mixing  rules.  Because  of 
using  complex  property  mixing  rules,  which  require  the  solution 
of  the  ratio  of  two  large  valued  determinants  to  get  a  small 
value  of  the  ratio,  they  are  known  to  give  erratic  results  under 
certain  circumstances.  For  example,  for  a  high  degree  of 
ionization,  the  heat  conductivity  coefficient  for  heavy  particles 
gives  a  very  small  negative  value.  Similarly  the  electrical 
conductivity  may  become  negative  for  a  small  degrees  of 
ionization.  These  later  programs  also  allow  calculation  of 
thermophysical  and  transport  properties  for  almost  thirty 
different  primary  components,  by  appropriately  combining  the  data 
for  these  different  components  in  the  data  file  PLASMA2.DAT. 

Therefore,  the  user  has  a  choice.  If  calculation  of  the 
properties  of  noble  gases  as  accurately  as  possible  at  all 
possible  temperatures  and  pressures  is  desired,  use  the  SEENIRAJ 
programs.  The  later  programs  can  also  be  used  for  the 
calculation  of  the  properties  of  noble  gases.  However,  for  any 
other  gas  plasma,  only  the  later  programs  are  to  be  used. 

All  of  the  programs  in  this  package  have  been  compiled, 
linked  and  run  with  the  help  of  Microsoft  software  (Microsoft 
FORTRAN  77  version  V3.13  dated  Aug.  5,  1983  +  Object linker 
version  V2. 01)  and  IBM  software  (IBM  Professional  Fortran  version 
1.21  and  corresponding  linker).  The  Microsoft  FORTRAN  compiler 
has  the  limitation  of  block  space  for  both  instructions  and  data, 
and  may  require  splitting  of  the  file  PLASMA1.FOR. 
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Alternatively,  you  may  compile  everything  with  the  help  of  IBM 
software  with./B  option.  The  1992  modification  used  Microsoft 
FORTRAN  version  5.1;  PLASMA1.F0R  did  NOT  have  to  be  split  if  this 
compiler  is  used.  A  memory  model  designated  LARGE  is  recommended 
if  this  compiler  is  used,  however.  Compilation  was  also 
attempted  on  Lahey  Personal  FORTRAN  version  2.0.  Although  it 
could  be  compiled,  it  could  not  be  linked  due  to  a  64K  limit  on 
the  sum  of  all  object  segments  to  be  linked.  After  contacting 
Lahey,  it  appears  that  even  later  versions  of  the  compiler  may 
have  this  limitation,  and  it  is  recommended  that  Microsoft  or  IBM 
compilers  should  be  used,  unless  the  user  has  information  on  how 
to  get  around  the  limit  on  block  size.  For  the  purpose  of 
compatibility  with  both  of  the  above  software,  initially  an 
asterisk  refers  to  the  console  terminal  with  keyboard,  but  after 
calling  the  program  OUTPLT  by  the  main  program,  the  console 
terminal  with  the  keyboard  is  referred  to  as  unit  4  and  the 
output  file  (con  =  console;  prn  =  printer;  and  any  other  disk 
output  file)  is  referred  to  as  unit  6.  While  running  the 
programs,  the  input  data  files  PLASMA1.DAT  and  PLASMA2.DAT  must 
remain  in  the  default  disk  or  hard  drive  directory. 

There  are  a  number  of  test  programs  (main  programs)  being 
supplied  in  this  program  package,  which  can  be  compiled,  linked 
and  executed  to  familiarize  the  user  with  the  programs.  In  actual 
practice  the  main  program  in  the  particular  combination  has  to  be 
replaced  by  a  user  supplied  main  program. 

For  any  questions  about  the  methodology,  the  following 
references  may  be  consulted: 

1.  D.  Kannappan  and  T.K.  Bose,  Transport  properties  of 
two-temperature  argon  plasma.  Physics  of  Fluids,  Vol .  20,  No. 

10,  Part  1,  Oct.  1977,  pp.  1668-73. 

2.  D.  Kannappan  and  T.K.  Bose,  Transport  properties  of  two- 
temperature  helium  plasma,  Physics  of  Fluids.  Vol.  23,  No.  7, 

July  1980,  pp.  1473-74. 

3.  T.K.  Bose  and  R.V.  Seeniraj,  On  reactive  conductivity 
coefficient  of  multiple  ionized  two- temperature  argon 

plasma, 

Waerme-  und  Stof fuebertragung.  1985,  pp.  3-8. 

4.  T.K.  Bose,  Thermodynamic  analysis  of  a  seeded  magnetogasdyn- 
amic  combustion  plasma,  AIAA/ASME  4th  Joint  Thermophysics  and 
Heat  Transfer  Conf.,  Boston/Mass.,  1986,  AIAA-86-1333. 

5.  T.K.  Bose,  Thermophysical  and  transport'  properties  of  multi¬ 
ple  component  gas  plasma  at  multiple  temperatures.  Progress 
in  Aerospace  Sciences,  1987  (in  press). 

***************************************************************** 


D. 2  List  of  Files  and  Subroutines 


D. 2 . 1  Data  files ; 

(to  be  stored  in  default  disk  during  execution): 

PLASMA1.DAT  :  file  containing  collision  cross-section  data. 

PLASMA2.DAT  :  file  containing  individual  gas  data: 

1.  Ionization  energy,  statistical  weight  and 

excitation  energy  (referred  to  as  Part  I) 

2.  gas-kinetic  cross-section  for  electron-atom 

collision  (referred  to  as  Part  II) 

3.  enthalpy-entropy  data  (in  kJ/kmole  for 

temperature  upto  15000K) .  (referred  to  as 

Part  III) 

D.2.2  Subroutine  files : 

PLASMA1 . FOR  :  principal  subroutines  to  calculate  two- temperature 
plasma  properties . 

PLASMA2.F0R  :  auxiliary  subroutines  to  organize  calculation  of 
properties  containing  atoms  and/or  molecules. 

AIPLASMA. FOR:  auxiliary  subroutine  to  organize  calculation  of 
properties  of  air  up  to  15000K. 

SEENIRAJ . FOR:  auxiliary  subroutines  to  organize  calculation  of 
properties  of  noble  gases  by  previous  method. 

SEENMICR. FOR:  auxiliary  subroutines  to  organize  calculation  of 

properties  of  noble  gases  by  previous  method  while 
accounting  for  the  effects  of  anomalous  electrical 
conductivity  resulting  from  the  presence  of 
lower  hybrid  microinstabilities. 

D.2.3  Demonstration  main  routine  files: 

SEENTEST . FOR:  to  calculate  properties  of  noble  gases  by  earlier 
method  (link:  SEENTEST+SEENIRAJ+PLASMA1 ) . 

NOBLE. FOR  :  to  calculate  properties  of  noble  gases  by  later 
method  (link:  NOBLE +PLASMA1+PLASMA2) . 

AIRTEST . FOR  :  to  calculate  properties  of  air  plasma(up  to  15000K) 
(link:  AIRTEST+AIPLASMA+PLASMA1 ) . 

ATOMTEST . FOR :  to  calculate  properties  of  atomic  and  ionic  air 
plasma  (beyond  15000K  to  50000K) ( link:  ATOMTEST+ 
PLASMA1+PLASMA2 ) . 

MGDTEST . FOR  :  to  calculate  properties  of  seeded  combustion  plasma 
(upto  6000K) ( link:  MGDTEST+PLASMA1+PLASMA2 ) . 

MICRTEST  :  to  calculate  electrical  conductivity  both 

classically  and  including  anomalous  transport  due 
to  lower  hybrid  microinstabilities.  The  other 
transport  properties  are  only  calculated 
classically  by  the  SEENIRAJ  method. 

While  the  above  demonstration  programs  are  all  interactive  in 

nature,  and  the  results  are  either  displayed  on  the  CRT  terminal 

or  printed  on  the  printer,  some  additional  information  follows 

which  may  be  quite  useful: 
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AIRTEST . FOR:This  is  to  determine  air  plasma  composition  for  a 
single  temperature  plasma.  Calculation  is  done  in 
three  ranges  up  to  1500QK. 

Range  I:  N2,02,N,0,N0 

Range  II:  N2 , 0+ . e- , N, 0, 02 , NO, N+ 

Range  III:  N+ , 0+ , e- , N, 0 

ATOMTEST . FOR :  This  is  for  a  multi-temperature  air  plasma  in  the 
temperature  range  beyond  that  in  AIRTEST.  The 
species  being  considered  are: 
e-,N+,0+,N,0,N++,0++ 

With  increasing  temperature,  ionization  level  is 
automatically  changed  and  the  correspondingly 
higher  ionized  nitrogen  and  oxygen  ions  are 
obtained . 

MGDTEST.FOR  :  This  is  for  calculation  of  the  plasma  for 

combustion  of  different  fuel  gas  (Water  gas,  Lurgi, 
Koppers-Totzek ,  Producer  gas  and  methane)  with 
oxygen  enriched  air  and  alkali  metal  seed  (either 
pure  metal  or  specific  salt).  The  following  gas 
components  are  considered: 

C02 , 02 , N2 , H2 , K+ , e- , CO, 0, N, NO, H, H20  and  CH4 

NOBLE. FOR  :  For  the  calculation  of  multi- temperature  noble  gas 
plasmas  ( He , Ne , Ar , Xe, Kr )  at  different  temperatures 
and  pressures,  but  reactive  heat  conduction  is  not 
considered . 

SEENTEST . FOR  For  calculation  of  multi -temperature  noble  gas 
plasma  ( He, Ne, Ar, Xe, Kr )  at  different  temperatures 
and  pressures,  including  reactive  heat  conduction. 
***************************************************************** 

D.3  Content  of  each  file 

PLASMA1.DAT  :  non-dimensional  collision  cross-section  tables  for 
various  potentials  ( Lennard- Jones  6-12  potential, 
Monchik  tables  for  exponential  repulsive 
potential,  and  repulsive  and  attractive  Coulomb 

potential ) . 

PLASMA2.DAT  :  data  (Table  1)  for  various  gases  (Table  2). 

PLASMA1 . FOR  :  contains  the  subroutines  as  follows: 

RDITAB  :  to  read  data  from  PLASMA1.DAT 
RDECCS  :  to  read  .data  from  PLASMA2.DAT 
ELAT  :  to  calculate  electron-atom  collision 
cross-sections . 

LENJON  :  to  calculate  collision  cross-sections  for 
Lennard-Jones  6-12  potential 
MONCH  :  to  calculate  collision  cross-sections  for 
exponential  repulsive  potential  (Monchik 
tables) 

CRGCRS  :  to  calculate  collision  cross-sections  for 
attractive  or  repulsive  Coulomb  potential 
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IONAT  : 

EQUCOM  : 
XEGRDT  : 


DETERM  : 
AILGTB  : 


SOLEQ  : 
OUTPUT  : 
BROKAW  : 


TTGPTP  : 
PLASMA2 . FOR  :  ATOMPL  : 


MOLGAS  : 


AI PLASMA. FOR:  AIRPR1  : 

AIRPR2  : 


APRHLP  : 
SEEN I RAJ . FOR : RDSEEN  : 


ELATOM  : 
INSPPR  : 
ECMTTP: 
FCMTTP  : 


to  calculate  charge-transf er  cross-section 
for  atom-ion  collisions 
to  evaluate  equilibrium  mole  fraction 
to  evaluate  the  partial  derivative  with 
respect  to  Te  and  Th  (electron/heavy 
temperatures ) (used  for  single  species 
plasma  only) 

to  evaluate  determinants 

to  interpolate  in  tables  by  Aitken-Neville 
or  Lagrange  method  (WARNING:  an  error  in 
interpolation  was  observed  when  using  this 
subroutine  in  1992.  Unfortunately,  not 
enough  was  known  about  this  interpolation 
method  to  correct  it.  It  appears  that  the 
right  magnitudes  are  given,  but  the  sign 
is  SOMETIMES  wrong) 

to  solve  simultaneous  linear  equations 
to  define  output  files 

to  calculate  mixture  transport  properties 
and  internal  heat  conductivity  coefficient 
for  molecules  by  the  Brokaw  formulas 
general  organization  routine  to  calculate 
transport  properties  of  the  plasma 
general  organization  routine  to  calculate 
thermophysical  and  transport  properties 
for  gas  plasmas  with  multi-species  atomic 
components  (change  of  ionization  level  for 
individual  components;  reactive  heat 
conductivity  coefficient  included,  but 
radiation  not  included) 

general  organization  routine  to  calculate 
thermophysical  and  transport  properties 
for  gas  plasmas  in  which  there  is  at  least 
one  molecular  component  (change  of  ioniza 
tion  level  not  allowed;  reactive  heat 
conductivity  coefficient  included,  but 
radiation  not  included),  (valid  upto  15, 
000K) 

initialization  subroutine  for  air  plasma 

air  properties  calculation  subroutine 

(valid  upto  15000K) 

auxiliary  subroutine  to  APRHLP 

to  read  from  PLASMA1.DAT  and  values  for  a 

single  noble  gas  from  the  file  PLASMA2.DAT 

for  use  in  Seeniraj  routines 

to  calculate  electron-atom  collision  cross 

section 

to  calculate  for  individual  gas  components 
the  transport  properties 
equilibrium  composition  for  noble  gas 
plasmas 

frozen  composition  mixture 
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enthalpy-entropy  for  noble  gas  plasma 
TPSSIP  :  organization  subroutine  for  calculation  of 
transport  properties  of  a  noble  gas  plasma 
MGDTEST.FOR  SEEDCP  :  organization  subroutine  for  selected 

seeded  combustion  plasma 

***************************************************************** 


For  noble  gas  plasma  calculations  ( SEENTEST , MICRTEST ,  or  NOBLE), 
a  guess  of  the  electron  mole  fraction  for  electrons,  xe,  is 
required.  This  may  be  done  by  starting  from  a  sufficiently  high 
temperature  to  one,  where  error  is  found.  As  an  example,  at  1  bar 
and  electron  temperature  =  heavy  particle  temperature,  the 
following  guess  values  of  xe  for  different  noble  gases  are  given: 

He:  7 . 587e-7  at  8000K;  Ne :  8.631e-8  at  6000K;  Ar:  8.713e-7  at 
5000K;  Xe:  4.560e-8  at  4000K;  and  Kr:  2.569e-4  at  4000K. 

SEENTEST  calculates  the  transport  properties  over  a  range  of 
electron  temperatures,  pressures,  and  temperature  ratios.  The 
inputs  require  the  following  units,  and/or  ranges: 

Number  of  Temperatures:  This  is  the  number  of  temperatures  at 
(No.  of  Temp)  which  the  calculation  of  transport 

properties  are  desired.  The  number  can 
range  between  1  to  32767,  but  they 
increase  linearly  according  to  the 
temperature  difference  entered. 

Initial  Temperature:  This  is  the  first  temperature  at  which 

(Tempi)  calculation  of  transport  properties  are 

desired.  The  units  are  in  Kelvin. 
Temperature  Difference:  This  is  the  temperature  increment 
(Temp.  Diff.)  between  the  temperatures  transport 

properties  are  calculated  at.  The  units 
are  in  Kelvin. 

Number  of  Temperature  Ratios:  This  is  the  total  number  of 
(No.  of  Temp.  Rat.)  temperature  ratios  at  which 

transport  properties  are  desired. 
However,  there  is  no  input  for  a 
starting  temperature  ratio.  It 
assumes  the  initial  value  is 
one  and  will  increase  from  there  by 
1  each  time. 

Number  of  Pressures  :  This  is  the  total  number  of  pressures  at 

(No.  of  Press.)  which  calculations  are  desired.  This 

number  can  range  from  1  to  32767. 

This  is  the  first  pressure  at  which 
transport  properties  are  desired.  The 
units  are  in  bar(l  bar  =  10**5  Pa) 


Initial  Pressure 
( Pressl ) 
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Pressure  Factor  :  The  program  increases  pressure 

(P-Factor)  logarithmically  rather  than  linearly,  so 

P-factor  represents  the  factor  the 
previous  pressure  is  multiplied  by  to 
get  the  next  value  of  pressure  at  which 
calculations  are  to  be  performed. 

MICRTEST  has  all  of  the  values  assigned  internally.  If  the  user 
would  like  to  change  them,  the  following  variable  names  are 
associated  with  the  above  properties.  The  same  variable  names 
are  used  in  SEENTEST. 

DTE:  Temperature  Difference 

NPRES:  Number  of  Pressures 

NTEMP:  Number  of  Temperatures 

NTR :  Number  of  Temperature  Ratios 

P1BAR:  Initial  Pressure 

PFAC:  Pressure  Factor 

TEMIN:  Initial  Temperature 

XE1:  Initial  guess  of  mole  fraction 

For  running  the  program  AIRTEST  limiting  values  of  N2LIM  and 
02LIM  are  l.e-5. 

D.4.2  Flowcharts: 

The  following  are  the  flow  charts  of  the  demonstration  routines: 

AIRTEST  < - >  AIRPR1  < - >  OUTPUT 

RDITAB  < -  PLASMA1.DAT 

RDECCS  < -  PLASMA2.DAT 

AIRPR2  < - >  APRHLP  < - >  EQUCOM 

TTGPTP  < - >  LENJON 

CRGCRS 

ELAT 

DETERM 

BROKAW 

ATOMTEST  < - >  OUTPUT 

ATOMPL  < - >  RDITAB  < -  PLASMA1.DAT 

RDECCS  < -  PLASMA2.DAT 

PARTI 

SOLEQ 

TTGPTP  < - >  LENJON 

CRGCRS 

ELAT 

DETERM 

BROKAW 

MGDTEST  < - >  OUTPUT 

MOLGAS  < - >  RDITAB  < -  PLASMA1.DAT 

RDECCS  < -  PLASMA 2 . DAT 

PARTI 

SOLEQ 
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TTGPTP  < - >  LENJON 

CRGCRS 

ELAT 

DETERM 

BROKAW 

MICRTEST  < - >  RDSEEN  < - >  RDITAB  < -  PLASMAl.DATt-PLASMA2.DAT 

OUTPUT 

ECMTTP 

TPSSIP  < - >  LENJON 

MONCH 

IONAT 

CRGCRS 

ELATOM 

DETERM 


NOBLE  < - >  OUTPUT 

ATOMPL  < - >  RDITAB  < 

RDECCS  < 
PARTI 
SOLEQ 
TTGPTP  < 


SEENTEST  < - >  RDSEEN  < - >  RDITAB  < 

OUTPUT 

ECMTTP 

TPSSIP  < - >  LENJON 

MONCH 

IONAT 

CRGCRS 

ELATOM 

DETERM 


PLASMA1.DAT 

PLASMA2.DAT 


>  LENJON 
CRGCRS 
ELAT 
DETERM 
BROKAW 

PLASMA1 .DAT+PLASMA2 .DAT 


D.5  Description  of  records  (lines  of  data)  in  data  files  ; 

Table  1:  Description  of  records  in  PLASMA2.DAT.  For  a  qualitative 
description  of  what  types  of  information  are  in  Parts  I,  II,  and 
III,  see  section  D.2.1  on  data  files  .  A  more  complete 
description  of  what  variables  these  values  are  assigned  to  is 
given  in  the  explanations  of  variables  in  subroutine  RDSEEN  in 
section  D.6. 


Record  1:  NCODE,  NCARD ,  NCARDE,  NENER,  IDIM,  NCARDC,  NELATL, 
NCARDP,  NENP 
where 

NCOLE  :  code  number  of  gas 

NCARD  :  number  of  records  (lines  of  data)  for  each  gas 
:  NCARDE + NCARDC +NCARDP+ 2 
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NCARDE:  number  of  records  in  Part  I  for  each  gas 
NENER  :  number  of  ionization  levels 
IDIM  :  energy  dimension  code  (IDIM  =1,  in  J;  -2, 
in  cm**(-l);  =3,  in  °K;  = 4 ,  in  eV) 

NCARDC:  number  of  records  in  Part  II  for  each  gas 
NELATL:  number  of  electron-atom  gas-kinetic  energy 
levels 

NCARDP :  number  of  records  in  Part  III  for  each  gas 
NENF  :  number  of  enthalpy  (MJ/kmole)  and  entropy  (MJ/ 
kmole-K)  values  for  temperatures  upto  15,000K 
in  3  ranges  (15,000  to  6,000K:  temp,  interval  = 
500K;  6,000  to  3000K:  temp.  interval  =  200K; 

3,000  to  400K:  temp.  interval  =  100K;  298. 15K; 
and  OK) 

Record  2:  Text 

Record  NCARDE  (Part  I  for  each  gas): 

total  number  of  ioniz . level ;  no.  of  gi,Ei  in  each  level 
ionization  potential 
(gi,Ei)  values  in  each  level 
Record  NCARDC  (Part  II  for  each  gas): 

number  of  electron-atom  gas-kinetic  energy  ranges 
initial  energy  in  the  range;  energy  step  size 
electron-atom  gas-kinetic  collision  cross-section 
Record  NCARDP  (Part  III  for  each  gas): 

NENP  enthalpy  values 
NENP  entropy  values 
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Table  2: 

List  of 

gases  for 

which  data  is  in 

PLASMA2.DAT 

Code  No. 

Gas 

Part 

Code 

No.  Gas 

Part 

1 

He 

I,  II 

16 

C12 

II 

2 

Ne 

I,  II 

17 

NO 

1 1 ,  II I 

3 

Ar 

I  - 1 1 1 

18 

N20 

II 

4 

Kr 

I,  II 

19 

NH3 

II ,  III 

5 

Xe 

I,  II 

20 

H20 

II,  III 

6 

Li 

I,  II 

21 

CO 

II ,  III 

7 

Na 

I,  II 

22 

HC1 

II 

8 

K 

1 ,  1 1 

23 

CH4 

II ,  III 

9 

Cs 

I,  II 

24 

C02 

II, III 

10 

H 

I  - 1 1 1 

25 

OH 

III 

11 

N 

I-I II 

26 

Ar* 

III 

12 

0 

I  - 1 1 1 

27 

H+ 

III 

13 

H2 

II ,  III 

28 

N+ 

III 

14 

02 

II , III 

29 

0+ 

III 

15 

N2 

II , III 

30 

e- 

III 

***************************************************************** 

D.6  Description  of  Subroutines: 


AILGTB 

Call  for  the  subroutine  is: 

CALL  AI LGTB ( XT , YT , NT , NZ , X , Y , N , I FL , I AL ) 


where 


XT ( NT )  :  Table  of  x-values.  INPUT 

YT( NT , NZ )  :  Two-dimensional  array  of  y-values.  INPUT 
NT:  number  of  rows  in  the  table  of  Y  values  and  number  of  x 
values 

NZ  :  number  of  columns  in  the  table  of  Y  values.  INPUT 
X  :  Value  of  x  at  which  interpol.  values  are  required.  INPUT 
Y ( NZ )  :  interpolated  y-values.  OUTPUT 
N  :  order  of  interpolation.  INPUT 

IFL:  out  of  range  extrapolation  done/error.  (IFL  =  0,  no 
extrapolation;  =  -1,  extrapolation  at  near  end;  =  -2, 

extrapolation  at  far  end;  =  2,  monotonicity  not 
present:  error).  OUTPUT 

IAL :  index  for  interpolation  by  Aitken-Neville/Lagrange 
method  (IAL  =  0,  A . N .  method;  .NE.  0,  Lagrange 
method ) .  OUTPUT 

Internal  Variables: 

A:  This  variable  has  several  different  assignments: 

1)  The  difference  between  the  first  and  last  values  of 
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x 

2)  The  product  of  all  delta  x's  of  adjacent  points  in 

the  immediate  vicinity  of  the  desired  x 

3)  The  difference  between  adjacent  delta  x's 

B:  Product  of  all  possible  differences  between  delta  x's 
of  adjacent  points 

JA:  the  number  of  points  -  half  of  the  order  of 

interpolation 

JB:  the  number  of  points  +  half  of  the  order  of 

interpolation 

K:  Loop  variable  representing  a  column  of  Y  data 

Nl:  Half  of  the  order  of  interpolation 

Sign:  Indicator  as  to  whether  x-values  are  increasing  or 
decreasing( +1  =  increasing,  -1  =  decreasing) 

WA(N-l):  difference  between  the  current  x-value  and  the 
previous  one  examined 

WB(N-1,NZ):  Y  value  at  each  point  surrounding  the 
interpolation  point. 

WC(N1):  I  don't  know.  Ask  Dr.  Bose. 

WD ( N ) :  Ratio  of  product  of  delta  x's  over  product  of  all 
possible  differences  between  delta  x's 
YA:  difference  between  last  and  second  to  last  points 
YB:  difference  between  second  to  last  and  third  to  last 
points 

YC:  difference  between  last  and  third  to  last  points 
YD(NZ):  First  derivative  (Slope)  used  in  far-side 
extrapolation 

YDD(NZ) :  Second  derivative(numerical )  used  in  extrapolation 
AIRPR1 

Call  for  the  subroutine  is: 

CALL  AIRPR1 (NP, NPG, TT, CPRT, XPG) 


where 


NP  :  printing  index  (=0,  print  data  input  values;  ,NE.  0, 
do  not  print).  INPUT 

NPG(8,4)  :  contains  in  each  row:  charge  (for  molecules, 
value  put  iB  -2),  index  for  Part  I-III  data  in  Plasma2.DAT 
to  be  stored.  INPUT 

TT  :  temperature  table.  OUTPUT 
CPRT  :  specific  heat  (MJ/kmole-K)  table.  OUTPUT 
XPG( 8,8)  :  contains  in  each  row:  mole  mess,  Leonard- Jones 
potent ia  parameter, molecular  dia  (in  Ang.).  3  cols* 

INPUT.  On  OUTPUT  in  each  row  (cols.  4-8)  it  contains: 

mole  fraction,  molar  spec,  heat  (MJ/kmole-K) ,  molar  en 
thalpy  (MJ/kmole),  mean  free  path  (m),  and  species  dlf 
fusion  coeff.  (before  mix.  rule  is  applied) (m2/s ) . 


AIRPR2 
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Call  for  the  subroutine  is: 

CALL  AIRPR2 ( NPG, TT , CPRT , XPG , T , PBAR , IRANGE) 


where 


NPG, TT, CPRT, XPG  :  see  AIRPR1 
T  :  temperature  (deg.K).  INPUT 
PBAR  :  pressure  in  bar.  INPUT 

IRANGE  :  an  index  to  indicate  extrapolation  in  the  inter 
polation  subroutine.  INPUT 

APRHLP 

Call  for  the  subroutine  is: 

CALL  APRHLP(NTOT, NPRIM, XNMX , HMX , SV, CPRV, X, XPG, NPG, IND, T , 
PBAR, EPLN ) 


where 


NPG, XPG, T, PBAR  :  see  AIRPR1  and  AIRPR2 

NTOT  :  total  number  of  components.  INPUT 

NPRIM  :  number  of  primary  components.  INPUT 

SV, CPRV  :  specie  enthalpy,  entropy  and  spec.  heat.  INPUT 

X  :  mole  fraction.  INPUT/OUTPUT 

EPLN  :  relative  accuracy  of  calculation  of  mole  fraction. 
IND  :  index  table  to  refer  to  the  row  number  in  NGG,XPG. 
INPUT 


ATOMPL 


Call  for  the  subroutine  is: 

CALL  ATOMPL ( NTOT , NPRIM , NSPEC , NT EM I , I FLAG , NP , NBRKW , NPT , NCHRG , 
I ND , XPRIM, XPRAT , XMOL , POTX , SIGX , HSO , NPG , XPG , HMX , XNMX , PBAR , TE , 
TH , ENTHM , ENTRM , XMOLM , DENS , VI SCOS , XKTH , XKTE , GAMEH , ELCOND , PLFR 
, GAUSS ) 


where 


NTOT, NPRIM, NPG, XPG, PBAR, TE,TH  :  see  APRHLP 

NSPEC  :  number  of  element  specie.  INPUT 

NT EM I  :  direction  of  change  of  temperature  (increasing: 

NT EM I  =  1;  decreasing:  NTEMI  =  -1).  INPUT 
I FLAG  :  an  index  ( I FLAG  =  -1:  read  data  through  RDITAB  and 
RDECCS;  I FLAG  =  0:  initialize  guess  values  and  other 
variables;  and  I FLAG  =  1:  normal  calculation).  INPUT 

NP  : printing  index  (=1:  only  final  results;  =  0:  no  print; 
=  -1:  print  all  values  except  input  data;  and  =  -2: 
print  everything).  INPUT 
NBRKW  :  =  0 


D-13 


NPT  :  table  of  data  input  code  for  different  components  used 
by  RDECCS .  INPUT 

NCHRG  :  initial  charges  of  NSPEC  elements 
IND  :  index  table  used  by  Subr.  TTGPTP 

XPRIM  :  guess  mole  fraction  of  primary  comp.  , incl.  elec.) 
XPRAT  :  ratio  of  mole  fraction  of  elements 
XMOL, POTX, SIGX,HS0  :  for  each  element  mole  mass,  Lennard- 
Jones  potential  function,  molecular  diameter,  and  heat 
of  formation  at  absolute  temperature.  INPUT 
ENTHM  :  mix.  enthalpy  (MJ/kg).  OUTPUT 
ENTRM  :  mix.  entropy  (MJ/kg-K).  OUTPUT 
XMOLM  :  mix.  mole  mass  (kg/kmole).  OUTPUT 
DENS  :  mix.  density  (kg/m3).  OUTPUT 
VISCOS  :  dyn.  viscosity  coeff.  (kg/ms).  OUTPUT 
XKTH  :  total  heavies  heat  cond.  coeff.  (kW/mK).  OUTPUT 
XKTE  :  total  electron  heat  cond.  coeff.  (kW/mK).  OUTPUT 
GAMEH  :  coll,  frequency  (s-1).  OUTPUT 
ELCOND  :  elec.  cond.  (A/Vm).  OUTPUT 
PLFR  :  plasma  freq.  (s-1).  OUTPUT 

GAUSS  :  magnetic  induction  (Gauss).  INPUT.  (If  required, 

ECYCF  -  cyclotron  freq. /coll.  freq.  can  also  be  ob¬ 
tained  by  changing  the  subr.  arguments). 

BROKAW 

Call  for  the  subroutine  is: 

CALL  BR0KAW(NIN,M1,M2,NM0L, IELEC, IND ,XPG,WA,WB, VISCOS, XKCH) 

This  subroutine  Is  an  auxiliary  subroutine  for  TTGPTP  to  use 
Brokaw  mixing  formulas  for  calculation  of  the  transport 
properties,  where  there  is  at  least  one  molecular  component. 

CRGCRS 

Call  for  the  subroutine  is: 

CALL  CRGCRS ( TREF , TH , TE , PBAR , X E , QBAR , Z 1 , Z 2 , NRANGE ) 


where 


TE, TH  :  electron  and  heavy  particles  temperature,  respective 
ly.  INPUT 

TREF  :  reference  temperature  (=TE,  if  at  least  one  collision 
partner  is  electron;  otherwise  =  TH).  INPUT 
XE  :  electron  mole  fraction.  INPUT 
Z1,Z2  :  charge  for  collision  partners. 

QBAR  :  collision  cross-section  (Ang2).  OUTPUT 
PBAR  :  pressure  (bar).  INPUT 

Internal  Variables: 

NRANGE:  Out  of  range  indicator  requested 
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Pi:  3.14159265358979 

OB(ll):  Internal  values  of  collision  cross  section 
TEXA(26,11):  Attractive  coulomb  potential  omega  integrals 
TEXR(26,11):  Repulsive  Coulomb  potential  omega  integrals 
TS:  Dimensionless  Temperature  this  corresponds  to  X  when 
AILGTB  is  called  to  interpolate  values 
TTS3(26):  First  column  of  Repulsive  Coulomb  potential 

tables.  This  corresponds  to  XT  when  AILGTB 
is  called  to  interpolate  values. 

TTS4 ( 26 ) :  First  column  of  attractive  Coulomb  potential 

tables.  This  corresponds  to  XT  when  AILGTB 
is  called. 

XA:  Several  assignments: 

1)  A  portion  of  the  denominator  of  the  debye  length 

2)  Dummy  variable  =  In  of  XL 

3)  Area  of  influence  for  shielded  coulomb  potential 
XB:  Dummy  variable  =  XL**2 

XL:  dummy  variable  =  4*TS 
ZZ:  product  of  charges 

DETERM 


Call  for  the  subroutine  is: 


CALL  DETERM( A, MXD , N , DETL , SIGN, IER) 


where 


A ( MXD , MXD )  :  coefficient  matrix.  INPUT 
N  :  order  of  A.  INPUT 

DETL  :  natural  logarithm  of  absolute(determinant ) .  OUTPUT 

SIGN  :  sign  of  determinant.  OUTPUT 

IER  :  output  error  flag  (=0,  no  error).  OUTPUT 

ECMTTP 

Call  for  the  subroutine  is: 

CALL  ECMTTP ( I , ITC, IER , NP , PBAR , TE , TH, XMASS , XE, XI , XI I , XI I I , 

XN , EIGT , ENTHT , ENTRT , XMOLT , HJ , XKHR, XKER, XEDTE, XEDTH , EPS, NITMX 
, NS AHA, ND I ) 

where 


I  :  minimum  ionization  level  for  input,  but  the  actual 
ionization  level  for  output(=l  for  neutrals). 
INPUT/OUTPUT 

ITC  :  number  of  iteration  done.  OUTPUT 

IER  :  output  error  code  (=0,  no  error).  OUTPUT 

NP  :  printing  index  (>0,  no  print;  =  0,  final  values).  INPUT 

XMASS  :  mole  mass  of  the  heavy  particle.  INPUT 

XE, XI, XI I, XI 1 1  :  mole  fractions  of  electrons  and  heavy 
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species.  INPUT/OUTPUT 

XN  :  number  density  (ra'-1).  OUTPUT 

ENTHT , ENTRT , XMOLT  :  enthalpy  (J/kmole),  entropy ( J/kmole-K ) 
and  mole  mass  of  the  mixture  plasma.  OUTPUT 
HJ  :  individual  component  enthalpy  (kJ/kmole).  OUTPUT 
NITMX  :  max.  number  of  iteration.  INPUT 

NSAHA  :  index  for  using  Ker rebrock/Vei s  equilibrium  formula 
(=0,  use  Veis  formula).  INPUT 
NDI :  index  number  for  using  "lowering  of  ionization  poten 
tial”  (=0,  lowering  enabled).  INPUT 
EPS:  relative  tolerance  for  equilibrium  composition 
calculation  INPUT 

Internal  Variables 

A:  Dummy  variable  =  (l-i)/i 
AA:  Several  assignments 

1)  Dummy  variable  for  upper  energy  limit  (ELIM) 

2)  Dummy  variable  for  ionization  potential  (EIT(NENER)) 

3)  Dummy  variable  for  electron  mole  fraction  (XE) 

4)  Dummy  variable  for  heavy  mole  fraction  (XII) 

AB:  Working  variable  for  In  (i-l/i) 

AE:  Dummy  variable  =  XE+XI+XII 

ALXE :  Natural  log  of  electron  mole  fraction 

ALX2I:  Natural  log  of  XII 

AI :  Dummy  variable  for  Ionization  level  in  REAL  format 
B:  Dummy  variable  for  (l+i)/i 

DBK ( 2 ) :  delta  b  from  equation  (16)  in  reference  [5] 

DFS :  Difference  between  actual  and  desired  values  for  the 
Newton  iteration  method  of  solving  implicit  relations 
EFACT :  Energy  conversion  factor  to  convert  to  *K 
IDIM  =  1:  EFACT  =  *K/J 
IDIM  =  2:  EFACT  =  °K/cra**(-l) 

IDIM  =  3:  EFACT  *  °K/°K 
IDIM  =  4:  EFACT  *  °K/eV 

EIT(NENER) :  Array  of  Ionization  potential  of  each  energy 
level 

EIGT(NENER):  Sum  of  all  lower  ionization  potentials  up  to 
the  current  energy  level . 

ELIM:  Maximum  excitation  Energy  that  has  to  be  considered  to 
evaluate  the .partition  function 
FS:  Implicit  function  for  electron  mole  fraction  used  in 
Newton's  method  solution 

FSD:  Derivative  of  FS  for  Newton's  method  solution 
GKK( 2 , 2 ) :  Matrix  represented  by  equation  (16a)  in  reference 
[5) 

GT(NT):  Statistical  weight  of  E(i);  known  as  g(i)  in 

statistical  thermodynamics;  the  degeneracy  of  an 
energy  level 
IDIM:  Energy  units  code 
J  =  1 

cm**(-l)  a  2 
°K  =3 
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ev  =4 

IJ:  Loop  variable  for  newton  method 
IMAX:  Maximum  ionization  level 
IMIN:  Minimum  ionization  level 
IMN1:  Dummy  Variable  for  I 

K:  Loop  variable  representing  the  current  g;,E.  pair  for  the 
current  ionization  level 

Kl:  Upper  limit  of  DO  loop  corresponding  to  NENER-1 
K2:  Two  assignments 

1 )  Loop  variable 

2)  Flag  for  the  number  of  ionization  levels 

KA:  Loop  variable  representing  the  current  energy  level 
KB:  The  current  energy  level  +  i-1 
KC:  The  number  of  g.  , E;  pairs  in  the  energy  level 
corresponding  to  'KB 

KIT:  Dummy  variable  for  the  number  of  iterations  performed 

NAA:  Flag  for  AA  if  electron  mole  fraction  is  within  limits 

NENER:  Number  of  ionization  levels 

NT:  Number  of  g.  ,  E.  pairs  in  each  ionization  level 

NTEST :  Flag  f or' del tas ( =0 ,  small  delta,  =  1,  large  delta) 

ST(3):  Saha  functions  for  each  heavy  species 

TA:  Working  Variable  for  Temperature  Ratio 

TE:  Electron  Temperature 

TH:  Heavy  Particle  Temperature 

Third:  Decimal  value  for  1/3 

UT ( NT , NENER ) :  E.  value  for  each  energy  level 

XA:  2  assignments 

1)  Dummy  variable  for  UT 

2)  coefficient  in  calculation  of  XEDTE,  XEDTH 
XB:  two  assignments 

1)  Dummy  variable  for  individual  component  enthalpy, HJ 

2)  coefficient  used  in  calculating  XEDTE, XEDTH 
XC:  coefficient  used  in  calculating  XEDTE, XEDTH 

XD:  coefficient  used  in  calculating  XEDTE, XEDTH 
XEDTE:  Partial  derivative  of  XE  w/  respect  to  TE 
XEDTH:  Partial  derivative  of  XE  w/respect  to  TH 
XEN:  New  guess  for  electron  mole  fraction  in  newton 
iteration 

XEQ:  Equilibrium  electron  mole  fraction 
XEX:  Dimensionless  energy  =  -E./KT. 

XEX1:  exp(-E,/KTe) 

XKER:  Intermediate  value  of  reactive  heat  conductivity  for 
electrons . 

XKHR:  Intermediate  value  of  reactive  heat  conductivity  for 
heavy  particles. 

XLD:  Lowered  mole  fraction,  maybe?  (from  lowering  of 
ionization  potential) 

XMASS:  Molar  mass  of  the  heavy  particle 
XMOLT:  Mole  mass  of  the  mixture  plasma 
XM:  Dummy  value  for  temperature  ratio 
XN:  Number  density  in  m**(-3) 

ZLN(3):  Energy  weighted  partition  function 
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ZT(3):  Partition  function 

ELAT  &  EL ATOM 

Call  for  the  subroutines  is: 

CALL  ELAT ( TE , QBAR , KS ) 

CALL  ELATOM ( TE , QBAR ) 

where 


TE  :  electron  temperature  (K).  INPUT 

QBAR  :  electron-neutral  collision  cross-section.  OUTPUT 
KS  :  index  for  the  neutral.  INPUT 

Internal  Variables: 

AA:  REAL  value  of  ICOUN 
AAK:  REAL  value  of  ICOUN/2 
AKB:  Boltzmann’s  constant 

AIN(25):  Integrand  of  equation  (35)  , ref erence  (5) 

AINI(5):  Initial  energy  in  electron-atom  gas-kinetic  range 
AINTE:  Integral  in  equation  (35),  reference  [5] 

AN:  multiplier  (by  eV)  to  get  dimensionless  E/KT 
ANIT:  Dummy  variable  for  initial  energy  in  electron-atom  gas 
kinetic  range 

AQ (25,5):  Electron-Atom  gas-kinetic  collision  cross  sections 
for  given  range  and  level 
DH:  Dummy  value  for  energy  step  size 
DHI(5):  Energy  step  size 
E:  Charge  of  an  electron  in  coulombs 
EV(25):  Value  of  energy  at  a  given  step 
GG:  Dimensionless  energy  E/KT 
I:  Loop  Variable 

ICOUN:  I  don't  see  a  use  for  this  variable — ask  Dr*  Bose 
IJ:  Loop  variable 

IOE:  Error  Flag,  but  for  what,  I  have  no  idea — ask  Dr.  Bose 
J:  Counter  representing  current  electron-atom  gas-kinetic 
level 

JJ:  Multiplicative  counter 
L:  Counter  for  loop  variable  I  +  1 
M:  Column  number  of  collision  matrix 
MF(10):  (1+1)1,  where  i  =.1,2, ..9 

NELATL:  Number  of  electron-atom  gas-kinetic  energy  levels 
NELAT:  Dummy  variable  for  NELAT 
NOL:  Dummy  variable  for  NOV-1 

NOM( 5 ) :  Number  of  electron-atom  gas-kinetic  energy  ranges 
NOV  :  Dummy  value  for  NOM 

OBDMR:  I  Don't  know,  but  a  guess  is  the  denominator  tens  to 
convert  s  values  of  the  Q  matrix  to  array  form. 
OBNMR:  I  Don't  Know,  but  a  guess  is  the  numerator  term  to 
convert  s  values  of  the  Q  matrix  to  the  array 
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form . 

OMF:  Multiplication  f actors ( m+ 1 ) i 
Pi:  3.14159265358979 

Q(25):  Electron-Atom  Gas-Kinetic  Cross  section  for  a  given 
range  and  level--each  element  represents  a 
different  range  in  he  same  level 
S:  Floating  Point  version  of  M 
SEVEN:  Even  element  integrands 
SODD:  Odd  element  integrands 
TE:  Electron  Temperature 
XEX :  exp( -E/KT ) 


EQUCOM 

Call  for  the  subroutine  is: 

CALL  EQUCOM ( NIN , NTOT , NPRIM, DLNMX , X J , XLNX , XLNKX , XNU , H I J ) 


where 

NIN  :  dimension  of  arrays  (must  be  .LE.  NTOT).  INPUT 
NTOT  :  number  of  total  components.  INPUT 
NPRIM  :  number  of  primary  components.  INPUT 
DLNMX  :  relative  tolerance.  INPUT 
XJ  :  mole  fraction.  OUTPUT 

XLNX  :  guess  value  of  logarithm  of  primary  mole  fraction. 
INPUT 

XLNKX  :  logarithm  of  equilibrium  constant.  OUTPUT 
XNU  :  chemical  valency  matrix.  INPUT 
HIJ  :  compatibility  equations  matrix 

FCMTTP 

Call  for  the  subroutine  is: 

CAL L  FCMTTP ( I , ND I , NSAHA , PBAR , TE , TH , XMASS , XE , ENTHT , ENTRT , 
XMOLT , XEQ, XNQ , ENTHTQ , ENTRTQ , XMOLTQ ) 


where  the  arguments  are  equivalent  of  ECMTTP,  except  that  the 
last  letter  Q  refers  to  equilibrium  values. 

INSPPR 

Call  for  the  subroutine  is: 

CALL  INSPPR( INDEX, PBAR, Z,Q) 


where 


INDEX  :  signifies  type  of  collision  (=l,i-i;  =2, ( i  +  1 )-( i+1 ) 


=  3 , e-e;  =4,i-e;  =5,(i+l)-e;  =6, i-( i+1) ) ;  ( i=ionization 

index;  value  of  INDEX  refers  to  row  number  of  Z).  INPUT 
PBAR  :  pressure (bar ) .  INPUT 

Z  :  auxiliary  array,  for  which  columns  in  each  row  denote: 
(col . 1 : ref . temp . , col . 2 : ref .mole  mass, cols. 3  & 

4: charge; 

col . 5 : heat  cond.coef f . , col . 6 : viscos . coef f . , col . 7 :dif f . 
coeff . ) .  INPUT/OUTPUT 

0  :  collision  cross-section  (for  details  of  locations,  see 
I ON AT ) .  INPUT 


IONAT 

Call  for  the  subroutine  is: 

CALL  IONAT (TH, PROPT, QBAR, NRANGE) 

where 


TH  :  heavy  particles  temperature  (K).  INPUT 
NRANGE  :  out  of  range  indicator  requested.  INPUT 
PROPT: species  properties  array  ( col . 1 : heavies  mole  mass; 
col . 2 .neutral  particles  dia. ;col . 3: Lennard-Jones 
potential ; cols. 4  to  ll:data  for  use  of  Monchik 
tables) .INPUT 

QBAR  :  collision  cross-section  Q'*,s'(ang.  )  with  column¬ 
wise  distribution  of  (l,s)  as  follows: 


1=  123456789 


( 1 i s]=  (1.1)  (1,2)  (1,3)  (1,4)  (1,5)  (1,6)  (1,7)  (2,2)  (2,3) 


1=  10  11  12  13  14 


(l,s)=  (2,4)  (2,5)  (2,6)  (3,3)  (4,4) 


In  above,  all  columns  are  not  filled  up  by  all  the  subroutines. 
These  are  given  in  the  following: 


col.  no. 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

11 

12 

13 

14 

Subr . : 

LENJON 

x 

X 

X 

X 

X 

X 

X 

X 

MONCH 

X 

X 

X 

X 

X 

X 

X 

X 

X 

IONAT 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

CRGCRS 

X 

X 

X 

X 

X 

X 

X 

X 

n 

X 

X 

X 

X 

X 

ELATOM 

X 

X 

X 

X 

X 

X 

X 
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Internal  Variables: 

A:  Constant  =  propt(lO) 

B:  Constant  =  Propt(ll) 

Beta(8):  Summation,  determined  from  equation  (36a)  in 
reference  [5] 

Eta(8):  Summation,  determined  from  equation  (36a)  in 
reference  [5] 

Phil:  Potential  value  #1  (propt(6)) 

Phi2:  Potential  value  #2  (propt(8)) 

R01:  Radial  position  #1  (Propt(7)) 

R02:  Radial  Position  #2  (Propt(9)) 

QBA(14):  Monchik  Cross  section  #1 
QBB(14):  Monchik  Cross  Section  #2 
XA:  Dummy  Variable  for  In  (TH/MH) 

LENJON 

Call  for  the  subroutine  is: 

CALL  LENJON (TH, POT, SIGMA, QBAR, NFL, NRANGE ) 


where 


TH,0BAR, NRANGE  :  see  IONAT 

POT  :  Lennard- Jones  potential  parameter  (K).  INPUT 
SIGMA  :  molecular  diameter  (ang.).  INPUT 
NFL  :  out  of  range  flag.  OUTPUT 

MOL GAS 

Call  for  the  subroutine  is: 

CALL  MOLGAS ( NP , NPT , NINI , KHS , NHS , NCP , NTOT , NPRIM, I FLAG, NBRKW, 
NHS 1 , IND , NPG, XPG, HMX , XNMX , TC , PBAR ) 


where 


NP  :  optional  printing  index( =0, -2, -1 , 1  for  no  printing,  all 
printing,  all  printing  except  input  data,  only  final 
results).  INPUT 

NPT ( NINI , 2 )  :  reading  from  PLASMA2.DAT  for  NINI  number  of 
gases;  rowwise  col.l:  species  code  ICODE,  col.  2:  row 
no.  of  (NPG, XPG).  INPUT 

KHS (NHS 1,2)  :  NHS  (=NHS1.LE.5)  values  of  enthalpy-entropy 
data  are  to  be  calculated  and  stored  (in 
addition  to  directly  read  from  PLASMA2.DAT) 
for  atom  and  ion  (referred  to  by  KF  (for  atoms 
only)  and  KT  in  NPG  table,  referred  by  the  two 
column  values  per  row  in  KHS).  INPUT 
NCP  :  total  number  of  spec,  heat  to  be  computed.  INPUT 
NTOT, NPRIM  :  total  number  of  components,  prim. comp.  INPUT 
I FLAG  :  initialize  also,  if  IFLAG=0 .  INPUT 
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NBRKW  :  use  Brokaw  mixing  formula,  if  NBRKW  =  1.  INPUT 

IND(NTOT)  index  for  current  data  (compatible  with  NPG) 

NPG,XPG  :  see  ATOMPL 

TC  :  current  temperature  (K).  INPUT 

PBAR  :  pressure  (bar).  INPUT 

Output  results  are  obtained  with  the  help  of  COMMON  /RESULT/  as 
follows: 

ENTHT  :  mix.  enthalpy  (MJ/kg) 

ENTRT  :  mix.  entropy  (MJ/kg-K) 

XMOLT  :  mix.  mole  mass  (kg/kmole) 

VISCOS  :  dyn.  viscosity  coeff.  (kg/ms) 

XKT  :  total  heat  cond.  coeff.  (kW/mK) 

DENS  :  mix.  mass  density  (kg,/ra3) 

DEN LG  =  LoglO (DENS) 

PLFR  :  plasma  freq.  (s-1) 

GAMEH  :  coll.  freq.  (s-1) 

PRAN  =  Prf:  frozen  Prandtl  number 
ELCOND  :  elec.  cond.  (A/Vm) 

CPEFF  :  effective  spec,  heat  (MJ/kg-K) 

RATIO  =  cycl . freq. /coll . freq.  (calculated  on  the  basis  of 
the  value  of  GAUSS  (magnetic  induction  in  Gauss)  as  INPUT) 

XT  :  mole  fraction  of  different  specie 


MONCH 

Call  for  the  subroutine  is: 

CALL  MONCH ( TH , PH I O , RO , QBAR , NF L , NRANGE ) 


where 


TH, QBAR, NFL, NRANGE  :  see  LENJON 

PHIO, RO  :  species  parameter  for  using  Monchik  tables.  INPUT 
Internal  Variables: 

ALPHA:  Collision  parameter  =  In  4to/kT;  corresponds  to  x  when 
AILGTB  is  called 

FLS(9):  Value  in  denominator  of  equation  (34b)  in  reference 
[5] 

Pi:  =  3.14159265358979 

QB(9):  Internal  values  of  collision  cross  section 
TMONCH(50, 9) :  Exponential  repulsive  potential  omega 
integrals 

TTS2:  First  column  of  Monchik  Tables.  Corresponds  to  xt 
when  AILGTB  is  called 
XA:  effective  cross  sectional  area 

OUTPUT 


Call  for  the  subroutine  is: 
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CALL  OUTPUT 


PARTI 

Call  for  the  subroutine  is: 

CALL  PART I ( NCHRG , KF , TE , TH , XE , PBAR , XMOL , HSPR , SSPR , ZEXC , E ION , 
EIONC ) 

where 


NCHRG  :  charge  index  INPUT 

KF  :  column  number  of  ( stat .weight ,  excit .energy ) 

TE , TH, XE, PBAR, XMOL  :  elec . temp heavies  temp elec . mole  fra 
ction, pressure( bar ), heavies  mole  mass.  INPUT 
HSPR, SSPR  :  (e  thalpy, entropy )/universal  gas  const.  OUTPUT 
ZEXC  :  excitation  partition  function.  OUTPUT 
EION, EIONC  :  ioniz.  energy,  total  ioniz.  energy  (from  ground 
level).  OUTPUT 


RDECCS 

Call  for  the  subroutine  is: 

CALL  RDECCS ( INPUT , NP , ICODE , KF , KS , KT ) 

where 


INPUT  :  input  channel  unit  no.  INPUT 

NP  :  optional  print,  if  NP=0.  INPUT 

ICODE  :  species  code  no.  (for  PLASMA2.DAT).  INPUT 

KF , KS , KT  :  column  numbers  (for  storing  input  data).  INPUT 

RDITAB 

Call  for  the  subroutine  is: 

CALL  RDITAB 

Internal  Variables 

I:  Loop  variable 

Title(72):  Dummy  variable  for  the  title  of  the  type  of 
table 

TEXA(26,11):  Attractive  coulomb  potential  omega 

integrals 

TEXR(26,11):  Repulsive  coulomb  potential  omega 
integrals 

TLEN(82,9):  Lennard  Jones  6-12  Potential  Omega 
Integrals 

TMONCH( 50, 9 ) :  Exponential  repulsive  potential  omega 

integrals 

TTS1(82):  First  column  in  Lennard  Jones  tables 
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TTS2(50):  First  column  in  Monchik  tables  for 
exponential  repulsive  potential 
TTS3{26):  First  column  in  repulsive  coulomb  potential 
tables 

TTS4(26):  First  Column  in  attractive  coulomb  potential 
tables 


RDSEEN 

Call  for  the  subroutine  is: 

CALL  RDSEEN 
Internal  Variables: 

AINI (NELATL) :  Initial  energy  in  electron-atom  gas-kinetic 
range 

AQ( NOV, NELATL) :  electron-atom  gas-kinetic  collision  cross 

section  for  a  given  range  and  level. 

DHI( NELATL):  energy  step  size 

EIT(NENER):  Ionization  potential  in  cm,  (see  the  CRC 

handbook  of  Chemistry  and  Physics,  pp.  10-210  to 
10-211  in  1991  version.) 

GT ( NT ( I ) ,  NENER ) :  degeneracy  of  each  energy  level (g,;  value) 

for  each  ionization  level. 

I:  Loop  variable  (This  subroutine  only) 

IC:  Flag  for  the  gas  used 
=1,  argon 
=2,  Helium 
=  3 ,  Neon 
=4,  Xenon 
=5,  Krypton 

ICT:  Gas  code  read  from  data  statement 
=1,  Helium 
-2,  Neon 
=3,  Argon 
=4,  Krypton 
=5,  xenon 

IDIM:  energy  dimension  code 


J  =  1 
cm**(-l)  =  2 
°  K  =3 
ev  =4 


NCODE:  Code  number  of  gas 
=1,  Helium 
=2,  Neon 
=3,  Argon 
=4,  Krypton 
=  5 ,  xenon 

NCARD:  total  number  of  lines  of  data  for  each  gas 
NCARD  =  NCARDE+ NCARDC  +  NCARD P + 2 
NCARDC :  number  of  lines  of  data  in  Part  II  of  Plasma2.Dat 
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for  each  gas 

NCARDE:  number  of  lines  of  data  in  Part  I  of  Piasma2.Dat 
for  each  gas 

NCARDP:  number  of  lines  of  data  in  Part  III  of  Plasma2.Dat 
for  each  gas 

NELATL:  Number  of  e lect ron-atom  gas  kinetic  energy  levels 

NENER:  Number  of  ionization  levels 

NENP:  Number  of  enthalpy (MJ/mole )  and  Entropy ( MJ/kMol - 0 K ) 

values  for  temperatures  up  to  15000  K  in  3  ranges 
15000-6000  K  :  500  K  intervals 

6000-3000  K  :  200  K  intervals 

3000-400  K  :  100  K  intervals 

also  at  298.15  and  0  K 

NOM( NELATL)  :  Number  of  electron-atom  gas  kinetic  ranges 

NOV:  Number  of  ranges  in  current  electron-atom  gas-kinetic 
energy  level 

NT ( NENER ) :  number  of  g, ,  E.  pairs  in  each  ionization  level 

PRO(ll,5):  Table  of  properties  for  the  5  noble  gases  read 
in  from  a  data  statement 

PROPT(ll):  Gas  properties  for  one  gas  (Determined  based  on 
IC) 

UT( NT ( I ) , NENER ) :  E,  value  for  each  ionization  level 


SOLEQ 

Call  for  the  subroutine  is: 

CALL  SOLEQ ( R,  A, MD, N, DET ) 


where 


R(MD)  :  right  hand  side  vector  (replaced  by  result).  INPUT 
A ( MD , MD )  :  coefficient  matrix  (destroyed).  INPUT 
N  :  number  of  equations  (N.LE.MD).  INPUT 
DET  :  determinant.  OUTPUT 

TPSSIP 

Call  for  the  subroutine  is: 

CALL  TPSSIP ( I , XE, XI , XI I , PBAR, TE, TH, PROPT, ELCOND, XKCE, XKCH, 
XMUT , XKRE , XKRH , XKHR , XKER, DAMB , Q, Z, GAMEH , N RANGE ) 


where 


I  :  ionization  index  (=1  for  neutrals).  INPUT 
XE, XI , XI I  :  mole  fraction  of  electrons  and  ions.  INPUT 
PBAR, TE, TH  :  pressure ( bar ) ,  elec . temp. , heavies  temp.  INPUT 
PROPT, Q, Z,NRANGE  :  see  IONAT  and  CRGCRS.  INPUT 
ELCOND  :  elec,  conductivity  ( amp/volt-m) .  OUTPUT 
XKCE, XKCH  :  heat  cond.  coeff.  (by  pure  conduction)  for  elec, 
and  heavies,  respectively  (kW/mK).  OUTPUT 
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XKRE , XKRH  :  respective  reactive  heat  cond.  coeff.  OUTPUT 
XMUT  :  viscosity  coeff.  (kg/ms).  OUTPUT 
XKHR,XKER  :  intermediate  values  from  ECMTTP. 

DAMB  :  ambi-polar  diff.  coeff.  (m'/s) .  OUTPUT 
GAMEH  :  electron-heavies  coll.  freq.  (s’1).  OUTPUT 

Internal  Variables: 

AI :  REAL  equivalent  of 

1)  i-1 

2)  i 

ASJK:  A-star  (j,k)  used  in  calculating  viscosity  and  thermal 
conductivity 

BSJK:  B-Star  (j,k)  used  in  calculating  thermal  conductivity 
C1,C2:  I  don't  know.  Ask  Dr.  Bose 

DEE:  4th  order  electron-electron  diffusion  coefficient 


DET: 

Thermo- 

diffusion  coefficient 

DTL1 

ln  of 

!QA| 

DTL2 

In  of 

1  QB ! 

DTL3 

In  of 

:qc! 

DTL4 

In  of 

IQDj 

IER1 

Output 

error 

flag 

in 

calculating 

:qai 

IER2 

Output 

error 

flag 

in 

calculating 

!  QB  i 

IER3 

Output 

error 

flag 

in 

calculating 

!QC ! 

IER4 

Output 

error 

flag 

in 

calculating 

!QD| 

IFL :  Don’t  know.  Ask  Dr.  Bose 

I FLAG:  Dummy  value  for  NIN  used  to  determine  type  of  cross- 
section  in  if-then  statements 
IB1:  Dummy  variable  for  row  of  Z 
IBIG:  Loop  variable  for  row  of  Z 
J:  Loop  variable  for  Row  of  QB 
Jl:  Dummy  value  for  row  of  QB 
J2:  Dummy  value  for  row  of  QB 
JA:  Loop  variable  for  row  of  Z 
JB:  Loop  variable  for  column  of  Z,Q 
K:  Loop  variable  for  column  of  QB 
K2:  Dummy  variable  for  column  of  QB 

NIN(6):  Flag  for  what  kind  of  Collision  cross  section  to 
use. 

1  =  >  Lennard  Jones  6-12  potential 
2=>  Monchik  Tables 

3=>  Charge  transfer  cross  section  for  ion-atom 
collisions 

4=>  Attractive  repulsive  coulomb  potential 
5=>  Electron-atom  cross  section 
NRANGE:  Out  of  range  indicator  for  Interpolation  of  tables 
PBAR:  Pressure  in  Bar 

PHI:  Value  used  in  calculating  exponential  repulsive 
potential  cross  sections 

POT:  Dummy  value  used  in  calculating  Lennard  Jones  6-12 
Potential  cross  sections 
PROPT(l):  Molecular  Mass  of  Heavy  particles 
PROPT( 2 ) :  Diameter  of  neutral  particles  in  angstroms 
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PR0PT(3):  Lennard- Jones  6-12  potential 
PROPT(4-12);  Data  for  use  in  Monchik  Tables 
Q:  Collision  cross  section 

QA(4,4):  4th  order  Sonine  Polynomial  expansion  in  matrix 
form 

QB(3,3):  First  minor  of  QA 
QC(3,3):  5th  minor  of  QA 
QD (3,3):  6th  Minor  of  QA 

RO:  Radial  position  for  use  in  Monchik  Tables 
SIGMA:  Dummy  value  for  diameter  of  neutral  particles  in 
angstroms 

SIGN1:  Sign  of  ! QA ! 

SIGN2:  Sign  of  i QB ! 

SIGN3:  Sign  of  ! QC ! 

SIGN4:  Sign  of  ! QD ! 

TR:  Temperature  Ratio 

TREF:  Reference  temperature  used  in  calculation  of  collision 
cross  sections 

U1,UY,UZ:  I  don't  know.  Ask  Dr.  Bose 
XA:  1)  weighted  heavy  particle  cross  section 
2)  Dummy  variable  for  XI, XII 
XE:  Electron  Mole  fraction 
XK:  I  don't  know.  Ask  Dr.  Bose 
XKD:  Correction  to  XKCE 

XI,XII,XIII:  mole  fractions  of  three  heavy  species  being 

examined.  XI  could  be  either  neutral  if  I  = 
1,  or  an  ion  mole  fraction  if  I>1. 

XMLR:  Reference  molecula.r  mass 
XN:  Number  density  in  m'J 
XTQ :  SQRT(TE) 

XXMU, YK, YXMU, ZK, ZXMU:  I  don’t  know,  ask  Dr,  Bose 
Z(6,7):  Auxiliary  matrix 

col  1:  reference  temperature 

col  2:  reference  molecular  mass 

col  3,4:  charge  of  colliding  partners 

col  5:  heat  conductivity  coefficient 

col  6:  viscosity  coefficient 

col  7:  Diffusion  coefficient 

row  1:  heavy  species  1-heavy  species  1  collision 

row  2:  heavy  species  2-heavy  species  2  collision 

row  3:  electron-electron  collision 

row  4:  heav'  species  1  -  electron  collision 

row  5:  heavy  species  2  -  electron  collision 

row  6:  heavy  species  1-heavy  species  2  collision 

For  the  modifications  made  for  microinstabilities,  the  following 
variables  were  also  used: 

CALDOCON:  effective  conductivity  calculated  by  mean  free 
path  formula  (used  by  Caldo  and  Choueiri) 

COND:  classical  conductivity  calculated  by  mean  free  path 
formula 
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DEEAN: effective  electron-electron  diffusion  coefficient 
DETAN:  Effective  thermo-diffusion  coefficient 
effcl:  ratio  of  effective  to  classical  resistivity 
elcondef:  effective  electrical  conductivity 
elcyc:  electron  cyclotron  frequency  in  radians/sec 
IERAN1 :  Output  error  flag  in  calculating  j QAAN ! 

IERAN2:  Output  error  flag  in  calculating  IQBANj 
IERAN3:  Output  error  flag  in  calculating  jQCANj 
IERAN4:  Output  error  flag  in  calculating  ! QDAN 
ji:  loop  variable  for  collision  type 
ki:  loop  variable  for  cross  section  type 
Hall:  electron  Hall  parameter 

GAMANP:  anomalous  momentum  exchange  collision  frequency 
QAAN:  Sonine  polynomial  expansion  for  effective  collision 
cross  sections  in  matrix  form 
QAN:  effective  collision  cross  section 

QANP:  anomalous  momentum  exchange  collision  cross  section 

QBAN (3,3):  First  minor  of  QAAN 

QCAN (3,3):  5th  minor  of  QAAN 

QDAN (3,3):  6th  Minor  of  QAAN 

SIGNAN1 :  Sign  of  ! QAAN  j 

SIGNAN2:  Sign  of  IQBANj 

SIGNAN3:  Sign  of  j  QCAN  j 

SIGNAN4 :  Sign  of  j  QDAN  j 

vthe:  electron  thermal  velocity 

vthi :  ion  thermal  velocity 

TTGPTP 

Call  for  the  subroutine  is: 

CALL  TTGPTP ( NN , Ml , M2 , M3 , I FL , TE  r  TH , PBAR , VISCOS , XKCH , XKCE , 
GAMEH , ELCOND , DEE , DET , NPG , XPG , WA , WB , WC , WD , GAUSS , ECYCF , RATIO, 
PLFR, IDN) 


where 


IDN(NN)  :  index  for  different  species  (to  refer  to  NPG,XPG) 
INPUT 

NPG(M1 , 4 ) , XPG (Ml , 8) :  see  ATOMPL.  INPUT 

(WA, WB) (M2, M2) , (WC, WD) (M3, M3 ) :  working  locations. 

TE,TH, PBAR, XKCH, XKCE, GAMEH, ELCOND  :  see  TTSSIP. 

VISCOS  :  viscosity  coeff.  .(kg/ms).  OUTPUT 

DEE, DET  :  electron  molar  and  thermo-diff.  coeff.  OUTPUT 

GAUSS  :  magnetic  field  (Gauss).  INPUT 

ECYCF  :  cyclotron  frequency.  OUTPUT 

RATIO  :  ratio  of  cyclotron  to  collision  freq.  OUTPUT 

PLFR  :  plasma  frequency.  OUTPUT 

XEGRDT 


Call  for  the  subroutine  is: 
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CALL  XEGRDT ( NCHRG, NSAHA , TE , TH, XE , ENER , XEDTE . XEDTH ) 
where 


NCHRG  :  charge  index.  INPUT 

NSAHA  :  Kerrebrock  (  .  NE .  0  )  ./Veis  (  =0  )  formulation  to  be 
TE, TH, XE  :  elec . temp heavies  temp. , elec .mole  frac. 
ENER:  ioniz.  energy  (K).  INPUT 

XEDTE, XEDTH  :  partial  derivative  of  XE  with  respect 
and  TH.  OUTPUT 


used . 
INPUT 

to  TE 


